組合數相關演算法合集

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.普通大數求法（擴充套件歐幾里得）

4.double精確版組合數