-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathBinomialCoefficients.cpp
More file actions
62 lines (51 loc) · 1.08 KB
/
BinomialCoefficients.cpp
File metadata and controls
62 lines (51 loc) · 1.08 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
//Binomial Coefficients - https://cses.fi/problemset/task/1079
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MOD = 1e9 + 7;
const ll INF = 1e18;
const int N = 1e6 + 17;
ll fact[N];
ll invfact[N];
ll expo(ll a, ll b, ll m) {
ll res = 1;
while (b > 0) {
if (b & 1) {
res *= a;
res %= m;
}
a *= a;
a %= m;
b >>= 1;
}
return res;
}
ll inverse(ll a, ll mod){
return expo(a , mod-2 , mod);
}
void precalc_factorials() {
fact[0] = fact[1] = 1;
invfact[0] = invfact[1] = 1;
for (int i = 2; i < N; ++i) {
fact[i] = (i * fact[i - 1]) % MOD;
invfact[i] = inverse(fact[i], MOD);
}
}
void solve() {
int a, b;
cin >> a >> b;
ll numerator = fact[a];
ll denominator = (invfact[b] * invfact[a - b]) % MOD;
cout << (numerator * denominator) % MOD << '\n';
}
int main() {
ios_base::sync_with_stdio(false);
cin.tie(nullptr);
precalc_factorials();
int T = 1;
cin >> T;
while (T--) {
solve();
}
return 0;
}