组合数（费马小定理， 快速幂）

给定 n 组询问，每组询问给定两个整数 a，b，请你输出 Cbamod(109+7)的值。

输入格式

第一行包含整数 n。

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

输出格式

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

数据范围

1≤n≤10000,
1≤b≤a≤105

输入样例：

3
3 1
5 3
2 2

输出样例：

3
10
1

使用快速幂求解a！，用费马小定理求解后面的逆元

代码：

#include <iostream>

using namespace std;

typedef long long ll;

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

ll n, a, b;
ll fact[N], ifact[N];

ll qmi(ll a, ll b) {
    ll res = 1;
    while (b) {
        if (b & 1) res = res * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return res;
}

void init() {
    fact[0] = 1;
    ifact[0] = 1;
    for (int i = 1; i < N; i ++ ) {
        fact[i] = fact[i - 1] * i % mod;
        ifact[i] = ifact[i - 1] * qmi(i, mod - 2) % mod;
    }
}

int main() {
    init();
    cin >> n;
    while (n -- ) {
        cin >> a >> b;
        cout << (fact[a] * ifact[b]) % mod * ifact[a - b] % mod << endl;
    }
}