求组合数

凡是对大数取余的 一般都是用卢卡斯逆元


#include <iostream>
#include<stdio.h>
using namespace std;
const int p=1e9+7;
typedef long long ll;
ll qmi(ll a,ll k)//用到快速幂/逆元/组合公式
{
    ll res=1;
    while(k)
    {
        if(k&1)
            res=res*a%p;
        a=a*a%p;
        k>>=1;
    }
    return res;
}
ll C(int a,int b)
{
    ll res=1;
    for(ll i=1,j=a; i<=b; i++,j--)
    {
        res=res*j%p;
        res=res*qmi(i,p-2)%p;
    }
    return res;
}
ll lucas(ll a,ll b)
{
    if(a<p&&b<p)
        return C(a,b);
    return
        C(a%p,b%p)*lucas(a/p,b/p)%p;
}
int main()
{
    ll m,n,i;
    ll ans=0,ans1=0;
    scanf("%lld%lld",&m,&n);
    
    printf("%lld\n",lucas(m,n));
    return 0;
}

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值