AcWing 886. 求组合数 II(组合计数+逆元++费马小定理+快速幂)

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const int mod=1e9+7,N=100010;
int fact[N],infact[N];

int qmi(int a,int p,int k){
	int res=1;
	while(p){
		if(p&1){
			res=(ll)res*a%k;
		}
		p=p>>1;
		a=(ll)a*a%k;
	}
	return res;
}

int main(){
	int n;
	cin>>n;
	
	fact[0]=1;
	infact[0]=1;
	for(int i=1;i<N;i++){
		fact[i]=(ll)fact[i-1]*i%mod;
		infact[i]=(ll)infact[i-1]*qmi(i,mod-2,mod)%mod;
	}
	
	for(int i=1;i<=n;i++){
	    int a,b;
		cin>>a>>b;
		cout<<(ll)fact[a]*infact[a-b]%mod*infact[b]%mod<<endl;
	}
	return 0;
}