求组合数C(n,m)的三种方法

一.暴力法
复杂度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);
}

配合逆元法使用即可