题目大意:发上来就过不了审核了……总之大意就是求C(n,m) mod 10007 m,n∈[1,2*10^8]
卢卡斯定理:C(n,m)=C(n%p,m%p)*C(n/p,m/p) mod p 要求p是质数
当中n%p可能会小于m%p 这样的情况下直接返回0就可以
证明去问卢卡斯 我不知道
#include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define p 10007 using namespace std; int fac[p],inv[p]; void Linear_Shaker() { int i; fac[0]=1; for(i=1;i<p;i++) fac[i]=fac[i-1]*i%p; inv[1]=1; for(i=2;i<p;i++) inv[i]=(p-p/i)*inv[p%i]%p; inv[0]=1; for(i=1;i<p;i++) inv[i]=inv[i]*inv[i-1]%p; } int C(int n,int m) { if(n<m) return 0; if(n<p&&m<p) return fac[n]*inv[m]%p*inv[n-m]%p; return C(n%p,m%p)*C(n/p,m/p)%p; } int main() { int T,n,m; Linear_Shaker(); for(cin>>T;T;T--) { scanf("%d%d",&n,&m); printf("%d ",C(n,m)); } }