组合数各种算法

组合数相关算法合集

1.大数取模P，P为素数（费马小定理）大佬代码1

const int mod = 1e9+7;  
const int maxn = 1e5+7;  
typedef long long ll;  
ll fac[maxn];  
ll qpow(ll a,ll b)  
{  
    ll ans=1;a%=mod;  
    for(ll i=b;i;i>>=1,a=a*a%mod)  
        if(i&1)ans=ans*a%mod;  
    return ans;  
}  
ll C(ll n,ll m)  
{  
    if(m>n||m<0)return 0;  
    ll s1=fac[n],s2=fac[n-m]*fac[m]%mod;  
    return s1*qpow(s2,mod-2)%mod;  
}  
  
  
//主函数里：  
fac[0]=1;  
    for(int i=1;i<maxn;i++)  
        fac[i]=fac[i-1]*i%mod;  

2.暴力打表版

const int mod = 1e9+7;  
const int maxn = 1e5+7;  
typedef long long ll;  
ll fac[maxn];  
ll qpow(ll a,ll b)  
{  
    ll ans=1;a%=mod;  
    for(ll i=b;i;i>>=1,a=a*a%mod)  
        if(i&1)ans=ans*a%mod;  
    return ans;  
}  
ll C(ll n,ll m)  
{  
    if(m>n||m<0)return 0;  
    ll s1=fac[n],s2=fac[n-m]*fac[m]%mod;  
    return s1*qpow(s2,mod-2)%mod;  
}  
  
  
//主函数里：  
fac[0]=1;  
    for(int i=1;i<maxn;i++)  
        fac[i]=fac[i-1]*i%mod;  

3.普通大数求法（扩展欧几里得）

大佬代码Download

4.double精确版组合数

大佬代码^_^