本人没想清为什么

fac[n]*(qpow(fac[n-k],mod-2)*qpow(fac[k],mod-2))%mod;

有些地方会输出负数而如果改成：

fac[n]*(qpow((fac[n-k]*fac[k])%mod,mod-2)%mod;

是精度问题还是一些奇怪的越界问题？

顺便警示一下后人，如果你 20pts 并且只过了第一个点和第四个点可以考虑一下这里。