组合数的几种求法

本文介绍在不同条件下计算组合数的方法,包括直接计算、使用卢卡斯定理等,并提供了多种情况下的具体实现代码。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Cnm%pC_n^m\%pCnm%p

  1. 0≤m≤n≤10000\leq m\leq n\leq 10000mn10001≤p≤1e91\leq p \leq 1e91p1e9,直接求
void Com(int n,int p){
	C[0][0]=1;
	for(int i=1;i<=1000;++i){
		for(int j=0;j<=i;++i){
			if(j==0||j==i) C[i][j]=1%p;
			else C[i][j]=(C[i-1][j-1]+C[i-1][j])%p;
		}	
	}
}
  1. 0≤m≤n≤1e180\leq m\leq n\leq 1e180mn1e181≤p≤1e61\leq p \leq 1e61p1e6,用卢卡斯定理
LL f[N];  //N为组合数的底数 的范围  
void init(int p){  
    f[0] = 1;  
    for(int i = 1; i <= p; ++i)  
        f[i]  = f[i-1] * i % p;  
}  
LL pow_mod(LL a, LL x, int p){  
    LL ret = 1;  
    a %= p;  
    while(x){  
        if(x & 1){  
            ret = ret * a % p;  
            --x;  
        }  
        else{  
            a = a * a % p;  
            x >>= 1;    
        }  
    }  
    return ret;  
}  
LL Lucas(LL n, LL k, int p){  
    LL ret = 1;  
    while(n && k){  
        LL nn = n % p, kk = k % p;  
        if(nn < kk) return 0;  
        ret = ret * f[nn] * pow_mod(f[kk] * f[nn - kk] % p, p - 2, p) % p;  
        n /= p;  
        k /= p;  
    }  
    return ret;  
}
  1. 0≤n≤1e180\leq n\leq 1e180n1e180≤m≤1e60\leq m\leq 1e60m1e61≤p≤1e91\leq p \leq 1e91p1e9,用卢卡斯定理
LL quick_mod(LL a, LL b)  
{  
    LL ans = 1;  
    a %= p;  
    while(b)  
    {  
        if(b & 1)  
        {  
            ans = ans * a % p;  
            b--;  
        }  
        b >>= 1;  
        a = a * a % p;  
    }  
    return ans;  
}  
  
LL C(LL n, LL m)  
{  
    if(m > n) return 0;  
    LL ans = 1;  
    for(int i=1; i<=m; i++)  
    {  
        LL a = (n + i - m) % p;  
        LL b = i % p;  
        ans = ans * (a * quick_mod(b, p-2) % p) % p;  
    }  
    return ans;  
}  
  
LL Lucas(LL n, LL m)  
{  
    if(m == 0) return 1;  
    return C(n % p, m % p) * Lucas(n / p, m / p) % p;  
} 
  1. 对于底数固定的,递推求所有组合数
C[0]=1;
for(int i=1;i<=n;++i)  
            C[i]=C[i-1]*(n-i+1)%mod*inv[i]%mod;
  1. 0≤m≤n≤1e60\leq m\leq n\leq 1e60mn1e61≤p≤1e51\leq p \leq 1e51p1e5,p可能为合数。因为p为合数好多东西不能用,所以直接对结果的每个质因子取模,可以避免除法。
    nefu628
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值