結果
問題 | No.117 組み合わせの数 |
ユーザー | もりを |
提出日時 | 2019-09-19 15:26:14 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 92 ms / 5,000 ms |
コード長 | 1,252 bytes |
コンパイル時間 | 1,366 ms |
コンパイル使用メモリ | 168,964 KB |
実行使用メモリ | 50,144 KB |
最終ジャッジ日時 | 2024-07-19 06:46:15 |
合計ジャッジ時間 | 1,902 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ソースコード
#include <bits/stdc++.h> using namespace std; using i64 = int_fast64_t; using ui64 = uint_fast64_t; #define REP(i, n) for (i64 (i) = 0; (i) < (n); ++(i)) #define FOR(i, a, b) for (i64 (i) = (a); (i) < (b); ++(i)) const int MOD = (int)1e9 + 7; vector<uint_fast64_t> fac, finv, inv; void COMinit(int MAX) { fac.assign(MAX, 0); finv.assign(MAX, 0); inv.assign(MAX, 0); fac[0] = fac[1] = finv[0] = finv[1] = inv[1] = 1; for (int i = 2; i < MAX; ++i) { fac[i] = fac[i - 1] * i % MOD; inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD; finv[i] = finv[i - 1] * inv[i] % MOD; } } uint_fast64_t C(int n, int r) { if (n < r or r < 0) return 0; return fac[n] * (finv[r] * finv[n - r] % MOD) % MOD; } uint_fast64_t P(int n, int r) { if (n < r or r < 0) return 0; return (fac[n] * finv[n - r]) % MOD; } uint_fast64_t H(int n, int r) { if (n == 0 and r == 0) return 1; return C(n + r - 1, r); } int main() { COMinit(2001001); int Q; scanf("%d\n" , &Q); while (Q--) { int N, K; char c; scanf("%c(%d,%d)\n", &c, &N, &K); if(c == 'C') cout << C(N ,K) << '\n'; if(c == 'P') cout << P(N, K) << '\n'; if(c == 'H') cout << H(N, K) << '\n'; } }