组合数取模

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
const int MOD=1e9+7;
LL n,m;

LL quick_mod(LL a, LL b) {
    LL ans=1;
    a%=MOD;
    while(b) {
        if(b&1) {
            ans=ans*a%MOD;
            b--;
        }
        b>>= 1;
        a=a*a%MOD;
    }
    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)%MOD;
        LL b=i%MOD;
        ans=ans*(a*quick_mod(b,MOD-2)%MOD)%MOD;
    }
    return ans;
}

LL Lucas(LL n, LL m) {
    if(m==0) return 1;
    return C(n%MOD,m%MOD)*Lucas(n/MOD,m/MOD)%MOD;
}

int main() {
    cin>>n>>m;
    cout<<Lucas(max(n,m),min(n,m))<<endl;
    return 0;
}

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值