結果
| 問題 |
No.117 組み合わせの数
|
| コンテスト | |
| ユーザー |
takkun_k_t
|
| 提出日時 | 2020-05-11 17:18:41 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
RE
|
| 実行時間 | - |
| コード長 | 1,336 bytes |
| コンパイル時間 | 1,380 ms |
| コンパイル使用メモリ | 168,068 KB |
| 実行使用メモリ | 6,812 KB |
| 最終ジャッジ日時 | 2024-07-18 14:36:50 |
| 合計ジャッジ時間 | 2,075 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | RE * 1 |
ソースコード
#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;
}
takkun_k_t