結果
問題 | No.117 組み合わせの数 |
ユーザー | muripoyotaro |
提出日時 | 2018-10-13 19:14:56 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 997 bytes |
コンパイル時間 | 1,874 ms |
コンパイル使用メモリ | 203,228 KB |
実行使用メモリ | 26,668 KB |
最終ジャッジ日時 | 2024-10-12 18:09:57 |
合計ジャッジ時間 | 3,495 ms |
ジャッジサーバーID (参考情報) |
judge / judge3 |
(要ログイン)
ソースコード
#include<bits/stdc++.h> typedef long long ll; using namespace std; const ll MOD = 1e9+7; const ll MAX= 1000001; ll fac[MAX], finv[MAX], inv[MAX]; void COMinit() { fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for (ll 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; } } //二項計数計算 ll COM(ll n, ll k) { if (n < k) return 0; if (n < 0 || k < 0)return 0; return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD; } ll permutation(ll n,ll k){ if (n < k) return 0; if (n < 0 || k < 0)return 0; return fac[n] * (finv[n-k]%MOD) % MOD; } ll H(ll n, ll k){ if(n==0&&k==0)return 1; return COM(n+k-1,k); } int main(){ cin.tie(0); ios::sync_with_stdio(false); int n; char a,b,c,d; int x,y; cin>>n; COMinit(); while(cin>>a>>b>>x>>c>>y>>d){ if(a=='C'){ cout<<COM(x,y)<<endl; } else if(a=='H'){ cout<<H(x,y)<<endl; } else{ cout<<permutation(x,y)<<endl; } } }