結果
問題 | No.1431 東大文系数学2020第2問改 |
ユーザー | st1vdy |
提出日時 | 2021-07-30 16:42:12 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 767 ms / 5,000 ms |
コード長 | 1,955 bytes |
コンパイル時間 | 1,695 ms |
コンパイル使用メモリ | 169,104 KB |
実行使用メモリ | 73,600 KB |
最終ジャッジ日時 | 2024-09-15 12:11:00 |
合計ジャッジ時間 | 14,378 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge6 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 116 ms
73,600 KB |
testcase_01 | AC | 117 ms
73,556 KB |
testcase_02 | AC | 116 ms
73,516 KB |
testcase_03 | AC | 564 ms
73,472 KB |
testcase_04 | AC | 113 ms
73,600 KB |
testcase_05 | AC | 730 ms
73,600 KB |
testcase_06 | AC | 702 ms
73,364 KB |
testcase_07 | AC | 687 ms
73,432 KB |
testcase_08 | AC | 691 ms
73,364 KB |
testcase_09 | AC | 767 ms
73,536 KB |
testcase_10 | AC | 459 ms
73,544 KB |
testcase_11 | AC | 443 ms
73,596 KB |
testcase_12 | AC | 472 ms
73,428 KB |
testcase_13 | AC | 326 ms
73,428 KB |
testcase_14 | AC | 323 ms
73,524 KB |
testcase_15 | AC | 312 ms
73,536 KB |
testcase_16 | AC | 265 ms
73,360 KB |
testcase_17 | AC | 249 ms
73,432 KB |
testcase_18 | AC | 245 ms
73,580 KB |
testcase_19 | AC | 756 ms
73,472 KB |
testcase_20 | AC | 249 ms
73,424 KB |
testcase_21 | AC | 582 ms
73,600 KB |
testcase_22 | AC | 339 ms
73,396 KB |
testcase_23 | AC | 759 ms
73,576 KB |
testcase_24 | AC | 757 ms
73,560 KB |
testcase_25 | AC | 116 ms
73,600 KB |
testcase_26 | AC | 115 ms
73,472 KB |
testcase_27 | AC | 116 ms
73,600 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; }