可能是求逆元太慢了,
请把你的求逆元代码改成:(N=2e7)
fac[0] = 1;
for (int i = 1; i <= N; i++) fac[i] = fac[i - 1] * i % mod;
inv[0] = 1, inv[N] = qpow(fac[N], mod - 2);
for (int i = N; i >= 1; i--) inv[i] = inv[i + 1] * (i + 1) % mod;
求C(n,m):fac[n] * inv[m] % mod * inv[n - m] % mod。