sumcomb (组合数)

10.19

找规律题目,在图上推推就发现,一个是起点下面的数,一个是起点右下角的数。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#define LL long long
using namespace std;

const LL mod = 1e9 + 7;

LL mub[1000010];
LL x, y;

void init(){
    mub[0] = 1;//
    for(int i=1; i<=1000010; i++){
        mub[i] = mub[i-1] * i % mod;
    }
}

LL exgcd(LL a, LL b, LL &x, LL &y){
    if(a == 0 && b == 0)
        return -1;
    if(b == 0){
        x = 1; y = 0;
        return a;
    }
    LL d = exgcd(b, a % b, y, x);
        y -= a / b * x;
    return d;
}

LL mod_reverse(LL a, LL n){
    LL d = exgcd(a, n, x, y);
    if(d == 1)
        return ( x % n + n ) % n;
    else
        return -1;
}

LL solve(int a, int b){
    if(a > b) return 0;
    LL nn = mod_reverse((mub[a] * mub[(b + mod - a) % mod]) % mod, mod);
    return mub[b] * nn % mod;
}

int main(){
    freopen ("sumcomb.in", "r", stdin);
    freopen ("sumcomb.out", "w", stdout);
    int T; scanf("%d", &T);
    init();
    while ( T-- ){
        int opt, a, b; scanf("%d%d%d", &opt, &a, &b);
        if(a > b) { printf("0\n"); continue;}
        if(opt == 1) printf("%I64d\n", solve(b+1, a+1));
        else printf("%I64d\n", solve(b, a+1));
    }
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值