[组合数]lucas定理 排列组合取模

博客介绍了lucas计算组合数C的两种处理方法。当n、m范围大而mod范围小时,可将n、m缩减至mod范围内计算;当n、m范围小而mod范围大时,可预处理阶乘逆元进行O(1)计算。

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

lucas: 计算 C 时 当 n, m范围大 mod范围小时可以将 n m 缩减至 mod 范围内计算

预处理: 当 n m 范围小 mod 范围大时可以预处理阶乘逆元 O(1)计算 

const int FN = 1e5 + 10;
const int MOD = 1e9 + 7;
ll fac[FN] = {1, 1}, inv[FN] = {1,1}, f[FN] = {1,1};
ll C(ll a,ll b)
{
	if(b > a) return 0;
	return fac[a] * inv[b] % MOD * inv[a - b] % MOD;
}
void init() //预处理逆元和阶乘
{
	for(int i = 2; i < FN; i++)
	{
		fac[i] = fac[i - 1] * i % MOD;
		f[i] = (MOD - MOD / i) * f[MOD % i] % MOD;
		inv[i] = inv[i - 1] * f[i] % MOD;
	}
}

lucas

 C  m 
    n  = lucas(a, b)

ll mulit(ll a, ll b, ll m)
{
    ll ans = 0;
    while(b)
    {
        if(b & 1)
            ans = (ans + a) % m;
        a = (a << 1) % m;
        b >>= 1;
    }
    return ans;
}
 
ll quick_mod(ll a, ll b, ll m)
{
    ll ans = 1;
    while(b)
    {
        if(b&1)
            ans = mulit(ans, a, m);
        a = mulit(a, a, m);
        b >>= 1;
    }
    return ans;
}
 
ll comp(ll a, ll b, ll m)
{
    if(a < b)
        return 0;
    if(a == b)
        return 1;
    if(b > a - b)
        b = a - b;
    ll ans = 1, ca = 1, cb = 1;
 
    for(int i = 0; i < b; i++)
    {
        ca = ca * (a - i) % m;
        cb = cb * (b - i) % m;
    }
 
    ans = ca * quick_mod(cb, m - 2, m) % m;
    return ans;
}
 
ll lucas(ll a,ll b,ll m)
{
    ll ans = 1;
    while(a && b)
    {
        ans = (ans * comp(a % m, b % m, m)) % m;
        a /= m;
        b /= m;
    }
    return ans;
}
 

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值