結果
問題 | No.117 組み合わせの数 |
ユーザー | takkun_k_t |
提出日時 | 2020-05-11 17:15:49 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 1,335 bytes |
コンパイル時間 | 1,671 ms |
コンパイル使用メモリ | 167,248 KB |
実行使用メモリ | 6,816 KB |
最終ジャッジ日時 | 2024-07-18 14:32:43 |
合計ジャッジ時間 | 2,430 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ソースコード
#include<bits/stdc++.h> #define rep(i,a,b) for(int i=a;i<b;i++) #define rrep(i,a,b) for(int i=a;i>=b;i--) using namespace std; typedef long long ll; ll mod = 1000000007; ll f[101010],rf[101010]; ll inv(ll x){ ll res=1; ll k=mod-2; ll y = x; while(k){ if(k&1) res = (res*y) % mod; y = (y*y) % mod; k /= 2; } return res; } void init(){ f[0]=1; rep(i,1,101010) f[i] = (f[i-1]*i) % mod; rf[101010 - 1] = inv(f[101010 - 1]); rrep(i, 101010 - 2, 0) rf[i] = (rf[i + 1] * (i + 1)) % mod; } ll C(int n, int k){ if (k < 0 || n < k) return 0; ll a = f[n]; // n! ll b = rf[n-k]; // (n-k)! ll c = rf[k]; // k! ll bc = (b*c) % mod; return (a*bc) % mod; } ll P(int n, int k) { if (k < 0 || n < k) return 0; ll a = f[n]; // = n! ll b = rf[n - k]; // = (n-k)! return (a * b) % mod; } ll H(int n, int k) { if (n == 0 && k == 0) return 1; return C(n + k - 1, k); } int main(){ init(); int T; cin >> T; rep(t, 0, T) { string s; cin >> s; char c; int a, b; sscanf(s.c_str(), "%c(%d,%d)", &c, &a, &b); int ans; if (c == 'C') ans = C(a, b); else if (c == 'P') ans = P(a, b); else ans = H(a, b); printf("%d\n", ans); } return 0; }