結果
問題 | No.117 組み合わせの数 |
ユーザー | pianoneko |
提出日時 | 2019-07-12 15:26:48 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 1,655 bytes |
コンパイル時間 | 1,724 ms |
コンパイル使用メモリ | 170,740 KB |
実行使用メモリ | 20,144 KB |
最終ジャッジ日時 | 2024-11-17 15:08:33 |
合計ジャッジ時間 | 3,732 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ソースコード
#include <bits/stdc++.h> #define REP(i, n) for (int i = 0; i < n; i++) #define ALL(obj) obj.begin(), obj.end() const int iINF = 1e9; const long long llINF = 1e18; const int MOD = 1e9 + 7; using namespace std; const int MAX = 1000000; vector<long long> fac(MAX + 1), inverseFac(MAX + 1); long long modPow(long long x, long long n) { long long ans = 1; while (n != 0) { if (n & 1) ans = ans * x % MOD; x = x * x % MOD; n = n >> 1; } return ans; } long long comb(long long a, long long b) { if (a == 0 && b == 0) return 1; if (a < b || a < 0) return 0; long long tmp = inverseFac[a - b] * inverseFac[b] % MOD; return tmp * fac[a] % MOD; } long long perm(long long a, long long b) { if (a < b) return 0; long long tmp = inverseFac[a - b] % MOD; return tmp * fac[a] % MOD; } long long dupc(long long a, long long b) { if (a == 0 && b == 0) return 1; if (a < 0 || b == 0) return 0; long long tmp = inverseFac[a - 1] * inverseFac[b] % MOD; return tmp * fac[a + b - 1] % MOD; } void init() { fac[0] = 1; inverseFac[0] = 1; for (long long i = 0; i < MAX; i++) { fac[i + 1] = fac[i] * (i + 1) % MOD; inverseFac[i + 1] = inverseFac[i] * modPow(i + 1, MOD - 2) % MOD; } } int main() { init(); int t; cin >> t; char com, n1, n2, n3; int n, k; REP(i, t) { cin >> com >> n1 >> n >> n2 >> k >> n3; if (com == 'C') { cout << comb(n, k) << endl; } else if (com == 'P') { cout << perm(n, k) << endl; } else { cout << dupc(n, k) << endl; } } }