blob: 02228e215487365cd08441b2a5c9bcfaee0a735c (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
constexpr ll MOD = 1e9+7;
inline ll pw(ll base, ll exp) {
ll res = 1;
while (exp) {
if (exp & 1) (res *= base) %= MOD;
exp >>= 1, (base *= base) %= MOD;
}
return res;
}
inline ll inv(ll x) { return pw(x, MOD-2); }
ll fact[MN] = { 1 }, ifact[MN] = { 1 };
inline ll nCr(int n, int k) { return fact[n]*ifact[k]%MOD*ifact[n-k]%MOD; }
for (int i = 0; i < N; ++i) fact[i+1] = (i+1ll)*fact[i]%MOD, ifact[i+1] = inv(fact[i+1]);
|
