Lucas定理:
主要是求$C_{n}^{m}$在模$p$情况下($mod \, p$)(一般$p$较小,而$n,m$较大的情况)
公式:
$ C_{n}^{m} ≡ C_{n \, mod \, p}^{m \, mod \, p} \times C_{n/p}^{m/p} (mod \, p) $
证明以后补吧
就以这题来说明具体解法:
题目
Luogu P3807 【模板】卢卡斯定理/Lucas 定理
Code:
//From:201929#include#define L long longusing namespace std;L pq[100005];L n,m,mod;L quick(L x,L y)//快速幂{ L ans=1; while(y) { if(y%2) ans=(ans*x)%mod; y>>=1; x=(x*x)%mod; } return ans;}L q(L x,L y){ if(y>x) return 0; return (pq[x]*quick(pq[y],mod-2))%mod*quick(pq[x-y],mod-2)%mod; //pq[x]/pq[y]/pq[x-y] //逆元变成pq[x]*(pq[y]^(mod-2))*(pq[x-y]^(mod-2)) //暴力求解C x,y的值}L lucas(L x,L y){ if(y==0) return 1ll; return (lucas(x/mod,y/mod)*q(x%mod,y%mod))%mod;}int main(){ int t; scanf("%d",&t); while(t--) { scanf("%lld%lld%lld",&n,&m,&mod); pq[0]=1; for(int i=1;i<=mod;i++) pq[i]=(pq[i-1]*i)%mod; //预处理小情况(x<mod)的答案 printf("%lld\n",lucas(n+m,m)); } return 0;}