結果
問題 | No.117 組み合わせの数 |
ユーザー | もりを |
提出日時 | 2019-09-19 15:37:02 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 74 ms / 5,000 ms |
コード長 | 1,336 bytes |
コンパイル時間 | 1,360 ms |
コンパイル使用メモリ | 168,536 KB |
実行使用メモリ | 50,688 KB |
最終ジャッジ日時 | 2024-07-19 06:58:11 |
合計ジャッジ時間 | 1,947 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ソースコード
#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)) constexpr bool EXPAND_TO_MINUS = true; constexpr int MOD = (int)1e9 + 7; constexpr int MAX_N = 2020202; uint_fast64_t fac[MAX_N], finv[MAX_N], inv[MAX_N]; void COMinit() { fac[0] = fac[1] = finv[0] = finv[1] = inv[1] = 1; for (int i = 2; i < MAX_N; ++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; if (n < 0) return EXPAND_TO_MINUS ? ((r % 2 ? -1 : 1) * C(-n + r - 1, r)) : 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(); 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'; } }