一.暴力法
复杂度O(n)
适用:小数据,n=60可用
代码:
f[0][0]=1;
for(int i=1;i<=n;i++)
for(int j=0;j<=i;j++)
if(j) f[i][j]=f[i-1][j]+f[i-1][j-1];
else f[i][j]=f[i-1][j];
二.逆元法
复杂度O(nlog(n))
适用:m小于p;
思路:把除以一个数变成乘一个数的逆元。
逆元求法
【对于正整数a和p,如果有ax≡1(mod p),那么把这个同余方程中x的小正整数解叫做a的逆元。由费马小定理知,当a与p互质时,a^(p-1) %p=1。所以当p是质数时,a的逆元为a^(p-2)。】
代码:
1)只需求一次大组合数时可使用
ll power(ll a,ll b){
ll ans=1;
while(b){
if(b&1) ans=ans*a%mod;
a=a*a%mod;
b>>=1;
}
return ans%mod;
}
ll C(ll n,ll m){//(n<=m)
if(n>m) {ans=0;return 0;}
ll ans=1;
n=min(n,m-n);
for(ll i=1;i<=n;i++){
ll res=power(i,mod-2);
ans=ans*res%mod;
}
for(ll i=m-n+1;i<=m;i++){
ans=ans*i%mod;
}
return ans;
}
2)求多次,优化一下,预处理阶乘和阶乘的逆元
ll f[N],inv[N];
ll pow_mod(ll a,ll b,ll c){
ll ans = 1;
ll base = a%c;
while(b){
if(b & 1) ans = (ans*base)%c;
base = (base*base)%c;
b >>= 1;
}
return ans%c;
}
inline ll C(int k, int n){
if(k > n) return 0;
return f[n]*inv[k]%mod*inv[n-k]%mod;
}
void init(){
//rev2 = pow_mod(2, mod-2,mod); // 2的逆元
f[0] = f[1] = 1;
for(ll i = 2; i < N; i++){ //预处理阶乘
f[i] = f[i-1]*i%mod;
}
inv[N-1] = pow_mod(f[N-1],mod-2,mod); //逆推预处理阶乘的逆元
for(int i = N-2; i >= 0; i--){
inv[i] = inv[i+1]*(i+1)%mod;
}
return ;
}
三.Lucas定理
复杂度:O(n*log(n)*log(n/mod))
适用:m>p,m,n很大时
卢卡斯(Lucas)定理的定义
- 设 P为素数,a,b∈N∗a,b∈N∗ ,并且
a=akp^k+ ak−1p^(k−1)+⋯+a1p+a0
b=bkp^k+ bk−1p(k−1)+⋯+b1p+b0
这里 0≤ai,bi≤p−1都是整数,i=0,1,2,⋯,k.
则有 C(a,b)≡C(ak,bk)∗C(bk−1ak−1)∗⋯∗C(a0,b0)(mod P)
再运用秦九韶算法可求得ai,bi;
代码:
ll Lucas(ll a,ll b){
if(a<mod&&b<mod)
return C(a,b);
return Lucas(a/mod,b/mod)*C(a%mod,b%mod);
}
配合逆元法使用即可