結果
問題 | No.1431 東大文系数学2020第2問改 |
ユーザー | st1vdy |
提出日時 | 2021-07-30 16:42:12 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 819 ms / 5,000 ms |
コード長 | 1,955 bytes |
コンパイル時間 | 1,679 ms |
コンパイル使用メモリ | 167,300 KB |
実行使用メモリ | 73,360 KB |
最終ジャッジ日時 | 2023-10-13 16:11:22 |
合計ジャッジ時間 | 15,397 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 117 ms
73,356 KB |
testcase_01 | AC | 117 ms
73,356 KB |
testcase_02 | AC | 115 ms
73,312 KB |
testcase_03 | AC | 629 ms
73,280 KB |
testcase_04 | AC | 118 ms
73,308 KB |
testcase_05 | AC | 807 ms
73,288 KB |
testcase_06 | AC | 728 ms
73,232 KB |
testcase_07 | AC | 759 ms
73,256 KB |
testcase_08 | AC | 727 ms
73,312 KB |
testcase_09 | AC | 773 ms
73,220 KB |
testcase_10 | AC | 478 ms
73,232 KB |
testcase_11 | AC | 426 ms
73,304 KB |
testcase_12 | AC | 495 ms
73,272 KB |
testcase_13 | AC | 330 ms
73,308 KB |
testcase_14 | AC | 338 ms
73,220 KB |
testcase_15 | AC | 323 ms
73,268 KB |
testcase_16 | AC | 268 ms
73,268 KB |
testcase_17 | AC | 244 ms
73,292 KB |
testcase_18 | AC | 248 ms
73,264 KB |
testcase_19 | AC | 810 ms
73,224 KB |
testcase_20 | AC | 249 ms
73,268 KB |
testcase_21 | AC | 617 ms
73,336 KB |
testcase_22 | AC | 348 ms
73,184 KB |
testcase_23 | AC | 810 ms
73,308 KB |
testcase_24 | AC | 819 ms
73,260 KB |
testcase_25 | AC | 117 ms
73,300 KB |
testcase_26 | AC | 118 ms
73,356 KB |
testcase_27 | AC | 115 ms
73,360 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; namespace SimpleModInt { constexpr int md = (int)998244353; inline int norm(long long a) { return (a % md + md) % md; } inline int add(int a, int b) { a += b; if (a >= md) a -= md; return a; } inline int sub(int a, int b) { a -= b; if (a < 0) a += md; return a; } inline int mul(int a, int b) { return (int)((long long)a * b % md); } inline int powmod(int a, long long b) { int res = 1; while (b > 0) { if (b & 1) res = mul(res, a); a = mul(a, a); b >>= 1; } return res; } inline int inv(int a) { a %= md; if (a < 0) a += md; int b = md, u = 0, v = 1; while (a) { int t = b / a; b -= t * a; swap(a, b); u -= t * v; swap(u, v); } assert(b == 1); if (u < 0) u += md; return u; } } using namespace SimpleModInt; constexpr int top = 3000 * 3000; vector<int> fac(top + 1), ifac(top + 1); int binom(int n, int m) { if (n < m || n < 0 || m < 0) return 0; return mul(fac[n], mul(ifac[m], ifac[n - m])); } void init() { fac[0] = 1; for (int i = 1; i <= top; i++) fac[i] = mul(fac[i - 1], i); ifac[top] = inv(fac[top]); for (int i = top; i; i--) ifac[i - 1] = mul(ifac[i], i); } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr); init(); int n, m, k, res = 0; cin >> n >> m >> k; for (int i = 0; i <= n; i++) { for (int j = 0; j <= n; j++) { int cnt = mul(binom((n - i) * (n - j), m), mul(binom(n, i), mul(binom(n, j), binom(i + j, k)))); if ((i + j - k) & 1) { res = sub(res, cnt); } else { res = add(res, cnt); } } } cout << res; return 0; }