求组合数(nlogn)

#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const int MAX=1e6+10;
using ll=long long;
const int mod=1e9+7;
long long fact[MAX],invfact[MAX];
ll ksm(ll a, ll b, int mod) {
    ll res = 1 % mod;
    while (b) {
        if (b & 1) res = res * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return res;
}

ll inv(ll x) {
    return ksm(x, mod - 2, mod) % mod;
}

void init(){
    fact[0]=invfact[0]=1;
    for(int i=1;i<=MAX;i++)
    {
        fact[i]=(fact[i-1]*i%mod);
        invfact[i]=(invfact[i-1]*inv(i))%mod;
    }
}
int main()
{
    std::ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    init();
    int n;
    cin>>n;
    while(n--)
    {
        int a,b;
        cin>>a>>b;
        cout<<(fact[a]*invfact[b]%mod*invfact[a-b]%mod)%mod<<"\n";//注意%mod 不然会爆
    }
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值