lucas定理 给自己打个板子!

本文介绍了一种利用 Lucas 定理加速计算组合数 C(n, m) 的算法,并通过 C(n, n+m) 的实例演示了如何在 C/C++ 中实现。该方法特别适用于大数模运算,通过取模和递归策略简化了计算过程。
摘要由CSDN通过智能技术生成

定理很有用
C(n,m)=C(n%mod,m%mod)*C(n/mod,m/mod);

#include<iostream>
#include<cstdio>
using namespace std;
#define ll long long
const ll int MAX_N = 1e5 + 10;
#define bug(x) cout<<#x<<" == "<<x<<endl;
ll mod ;
ll a[MAX_N]= {0};
ll epow(ll t)
{
    ll b = mod-2;
    ll ans = 1;
    while(b)
    {
        if(b&1)
            ans *= t;
        t *= t;
        ans %= mod;
        t %= mod;
        b >>= 1;
    }
    return ans;
}
ll C(ll n,ll m)
{
    if(m-n<0)
        return 0;
    return a[m]*epow(a[n])%mod*epow(a[m-n])%mod;
}
ll lucas(ll n,ll m)
{
    if(m<mod)
    {
        return C(n,m);
    }
    else
    {
        return C(n%mod,m%mod)%mod*lucas(n/mod,m/mod)%mod;
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        ll n,m;
        scanf("%lld %lld %lld",&n,&m,&mod);
        a[0] = 1;
        for(ll i=1; i<=100003; i++)
        {
            a[i] = a[i-1]*i;
            a[i] %= mod;
        }
        ll ans=lucas(n,n+m);
        printf("%lld\n",ans%mod);
    }
}

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值