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

给定 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;}
}