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<bits/stdc++.h> #define L long long using 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; }
标签:pq,return,Lucas,定理,ans,mod From: https://www.cnblogs.com/201929-whx/p/17606822.html