結果
問題 | No.117 組み合わせの数 |
ユーザー | Shibuyap |
提出日時 | 2020-08-03 05:13:31 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 268 ms / 5,000 ms |
コード長 | 1,680 bytes |
コンパイル時間 | 2,442 ms |
コンパイル使用メモリ | 202,668 KB |
実行使用メモリ | 52,664 KB |
最終ジャッジ日時 | 2024-07-22 07:38:25 |
合計ジャッジ時間 | 2,807 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ソースコード
#include <bits/stdc++.h> #define rep(i,n) for(int i = 0; i < (n); ++i) #define srep(i,s,t) for (int i = s; i < t; ++i) #define drep(i,n) for(int i = (n)-1; i >= 0; --i) using namespace std; typedef long long int ll; typedef pair<int,int> P; #define yn {puts("Yes");}else{puts("No");} #define MAX_N 200005 const int MAX = 2100000; const int MOD = 1000000007; long long fac[MAX], finv[MAX], inv[MAX]; // テーブルを作る前処理 void COMinit() { fac[0] = fac[1] = 1; finv[0] = finv[1] = 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; } } // 二項係数計算 long long COM(int n, int k){ if (n < k) return 0; if (n < 0 || k < 0) return 0; return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD; } long long FINV(int n){ if (n < 0) return 0; return finv[n]; } int main() { COMinit(); int T; cin >> T; rep(i,T){ string s; cin >> s; int m = s.size(); int n = 0, k = 0; int mode = 0; srep(i,2,m-1){ if(s[i] == ','){ mode++; }else if(mode == 0){ n *= 10; n += s[i] - '0'; }else{ k *= 10; k += s[i] - '0'; } } ll ans = 0; if(s[0] == 'C'){ ans = COM(n,k); }else if(s[0] == 'P'){ ans = COM(n,k) * fac[k] % MOD; }else{ ans = COM(n-1+k,k); } if(n==0&&k==0) ans = 1; cout << ans << endl; } return 0; }