【逆元求法】
1.要求:p 是质数
2.时间复杂度:O(n)
3.求解 的步骤:
1)通过循环,预先算好所有小于 N 的阶乘(%p)的结果,存到数组 fac[] 中 (fac[i] = i!%p)
2)求 的逆元(即求 fac[m] 的逆元),根据费马小定理,x%p 的逆元为 ,通过快速幂,求解 ,记为 M
3)求 的逆元:同上,即求解
4)通过逆元计算组合数,即:
4.实现:
LL powMod(LL x, LL n, LL mod) {//快速幂求x^n%mod
LL res=1;
while(n) {
if(n&1)
res=res*x%mod;
x=x*x%mod;
n>>=1;
}
return res;
}
LL inv(LL x,LL mod) {//求逆元
return powMod(x,mod-2,mod);
}
LL fac[N];
int main() {
LL n,m,mod;//要求mod是质数
scanf("%lld%lld%lld",&n,&m,&mod);
fac[0]=1;
for(int i=1;i<=n;i++)//预处理求fac,fac[i]=i!%mod
fac[i]=fac[i-1]*i%mod;
//C(n,m) = n!*(m!%mod的逆元)*((n-m)!%mod的逆元)%mod
LL res=fac[n]*inv(fac[m],mod)%mod*inv(fac[n-m],mod)%mod;
printf("%lld\n",res);
return 0;
}
有时需要用到 1~n 的所有逆元,这个时候可以利用递推来求
设
那么,有:
即:
两边同除 i*k,有:
即:
const int MOD = 1E9+7;
const int N = 100000+5;
LL fac[N],inv[N];
void init() {
fac[0]=1;
fac[1]=1;
inv[0]=1;
inv[1]=1;
for(LL i=2; i<N; i++) {
fac[i]=fac[i-1]*i%MOD;
inv[i]=(MOD-MOD/i)*inv[MOD%i]%MOD;
}
for(int i=1; i<N; i++)
inv[i]=inv[i]*inv[i-1]%MOD;
}
【递推打表】
1.要求:n、m 不大于 10000
2.时间复杂度:O(n^2)
3.方法:
4.实现
const int mod=1E9+7;
const int N=10000+5;
int comb[N][N];//comb[i][j]内存放的是C(i,j)%mod
void init() {
for(int i=0; i<N; i++) {
comb[i][i]=1;
comb[i][0]=1;
for(int j=1; j<i; j++) {
comb[i][j]=comb[i-1][j]+comb[i-1][j-1];
if(comb[i][j]>=mod)
comb[i][j]-=mod;
}
}
}