886. 求组合数 II

输入格式
第一行包含整数n。

接下来n行，每行包含一组a和b。

输出格式
共n行，每行输出一个询问的解。

数据范围
1≤n≤10000,
1≤b≤a≤10⁵

输入样例：
3
3 1
5 3
2 2

输出样例：
3
10
1

---

适用于： 数字较大，同时查询多（此时Ⅰ已经不再适用）
公式：
C(a,b) = fact[a] * infact[b] * infact[a-b]

方法：使用快速幂求逆元,预处理阶乘

1. infact：存储逆元的阶乘 O(logn)
   infact[i] = infact[i-1] * qmi(i,mod-2,mod)%mod;
2. fact：存储阶乘 O(1)
3. 所以从前往后枚举，计算阶乘 O(n)
   时间复杂度 1*3 --> O(nlogn)

注意：如果这个算式中出现除法，我们就需要逆元了，将除法运算转换为乘法运算。
具体参考：
https://www.cnblogs.com/czc1999/p/11682068.html

#include<iostream>
using namespace std;
typedef long long LL;

const int N=1e5+10,mod = 1e9+7;

LL fact[N],infact[N];


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

int main(){
    //阶乘预处理
    fact[0]=1;  //1的阶乘为1
    infact[0]=1;  //1的倒数阶乘为1
    for(int i = 1;i<N;i++){
        fact[i] = fact[i-1]*i%mod;
        infact[i]=infact[i-1]*qmi(i,mod-2,mod)%mod;  //(i-1)倒数阶乘 * i的倒数 = i 的倒数的阶乘
    }
    
    int n;
    cin>>n;
    
    int a,b;
    while(n--){
        cin>>a>>b;
        cout<<fact[a]*infact[b]%mod*infact[a-b]%mod<<endl;
    }
    
    return 0;
}