結果
問題 | No.117 組み合わせの数 |
ユーザー | akusyounin2412 |
提出日時 | 2018-04-03 00:48:33 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 326 ms / 5,000 ms |
コード長 | 2,828 bytes |
コンパイル時間 | 1,297 ms |
コンパイル使用メモリ | 120,744 KB |
実行使用メモリ | 34,432 KB |
最終ジャッジ日時 | 2024-06-26 06:55:06 |
合計ジャッジ時間 | 2,234 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ソースコード
# include <iostream> # include <algorithm> #include <array> # include <cassert> #include <cctype> #include <climits> #include <numeric> # include <vector> # include <string> # include <set> # include <map> # include <cmath> # include <iomanip> # include <functional> # include <tuple> # include <utility> # include <stack> # include <queue> # include <list> # include <bitset> # include <complex> # include <chrono> # include <random> # include <limits.h> # include <unordered_map> # include <unordered_set> # include <deque> # include <cstdio> # include <cstring> #include <stdio.h> #include <stdlib.h> #include <cstdint> using namespace std; using LL = long long; using ULL = unsigned long long; constexpr long long MOD = 1000000000 + 7; constexpr long long INF = std::numeric_limits<long long>::max(); const double PI = acos(-1); #define fir first #define sec second #define debug(x) cerr<<#x<<": "<<x<<'\n' typedef pair<LL, LL> Pll; typedef pair<LL, pair<LL, LL>> Ppll; typedef pair<LL, pair<LL, bitset<100001>>> Pbll; typedef pair<LL, pair<LL, vector<LL>>> Pvll; typedef pair<LL, LL> Vec2; struct Tll { LL first, second, third; }; typedef pair<LL, Tll> Ptll; #define rep(i,rept) for(LL i=0;i<rept;i++) #define Mfor(i,mf) for(LL i=mf-1;i>=0;i--) LL h, w, n, m, k, s, t, q, sum, last, cnt, ans,a; string str; struct Edge { LL to, cost; }; vector<Edge>vec; bool f; char c; void YN(bool f) { if (f) cout << "YES" << endl; else cout << "NO" << endl; } void yn(bool f) { if (f) cout << "Yes" << endl; else cout << "No" << endl; } struct Combination { int mod; vector< int64_t > mfact, rfact; Combination(int sz, int mod) : mfact(sz + 1), rfact(sz + 1), mod(mod) { mfact[0] = 1; for (int i = 1; i < mfact.size(); i++) { mfact[i] = mfact[i - 1] * i % mod; } rfact[sz] = inv(mfact[sz]); for (int i = sz - 1; i >= 0; i--) { rfact[i] = rfact[i + 1] * (i + 1) % mod; } } int64_t fact(int k) const { return (mfact[k]); } int64_t pow(int64_t x, int64_t n) const { int64_t ret = 1; while (n > 0) { if (n & 1) (ret *= x) %= mod; (x *= x) %= mod; n >>= 1; } return (ret); } int64_t inv(int64_t x) const { return (pow(x, mod - 2)); } int64_t P(int n, int r) const { if (r < 0 || n < r) return (0); return (mfact[n] * rfact[n - r] % mod); } int64_t C(int p, int q) const { if (q < 0 || p < q) return (0); return (mfact[p] * rfact[q] % mod * rfact[p - q] % mod); } int64_t H(int n, int r) const { if (n < 0 || r < 0) return (0); return (r == 0 ? 1 : C(n + r - 1, r)); } }; int main() { cin >> n; Combination comb(2000010, MOD); rep(i, n) { char ch; LL x, y; cin >> c >> ch >> x >> ch >> y >> ch; if (c == 'P') cout << comb.P(x, y) << endl; else if (c == 'C')cout << comb.C(x, y) << endl; else cout << comb.H(x, y) << endl; } return 0; }