結果
問題 | No.1923 Divisor Array |
ユーザー | りあん |
提出日時 | 2022-05-03 23:48:21 |
言語 | C++23 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 100 ms / 2,000 ms |
コード長 | 2,379 bytes |
コンパイル時間 | 3,944 ms |
コンパイル使用メモリ | 260,784 KB |
実行使用メモリ | 67,596 KB |
最終ジャッジ日時 | 2023-09-15 20:56:27 |
合計ジャッジ時間 | 8,485 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 30 ms
36,344 KB |
testcase_01 | AC | 30 ms
36,296 KB |
testcase_02 | AC | 100 ms
67,596 KB |
testcase_03 | AC | 31 ms
36,536 KB |
testcase_04 | AC | 68 ms
57,488 KB |
testcase_05 | AC | 50 ms
46,968 KB |
testcase_06 | AC | 78 ms
63,608 KB |
testcase_07 | AC | 51 ms
47,528 KB |
testcase_08 | AC | 53 ms
48,460 KB |
testcase_09 | AC | 85 ms
65,396 KB |
testcase_10 | AC | 47 ms
44,912 KB |
testcase_11 | AC | 50 ms
46,496 KB |
testcase_12 | AC | 57 ms
50,328 KB |
testcase_13 | AC | 45 ms
44,616 KB |
testcase_14 | AC | 30 ms
36,424 KB |
testcase_15 | AC | 30 ms
36,344 KB |
testcase_16 | AC | 30 ms
36,240 KB |
testcase_17 | AC | 30 ms
36,344 KB |
testcase_18 | AC | 30 ms
36,320 KB |
testcase_19 | AC | 31 ms
36,492 KB |
testcase_20 | AC | 31 ms
36,280 KB |
testcase_21 | AC | 32 ms
36,364 KB |
testcase_22 | AC | 30 ms
36,292 KB |
testcase_23 | AC | 31 ms
36,220 KB |
testcase_24 | AC | 61 ms
50,404 KB |
testcase_25 | AC | 71 ms
55,160 KB |
testcase_26 | AC | 81 ms
59,316 KB |
testcase_27 | AC | 75 ms
56,700 KB |
testcase_28 | AC | 55 ms
47,140 KB |
testcase_29 | AC | 48 ms
44,328 KB |
testcase_30 | AC | 61 ms
49,920 KB |
testcase_31 | AC | 78 ms
58,032 KB |
testcase_32 | AC | 80 ms
59,068 KB |
testcase_33 | AC | 68 ms
54,048 KB |
testcase_34 | AC | 31 ms
36,368 KB |
testcase_35 | AC | 30 ms
36,228 KB |
testcase_36 | AC | 30 ms
36,500 KB |
testcase_37 | AC | 31 ms
36,240 KB |
testcase_38 | AC | 32 ms
36,220 KB |
testcase_39 | AC | 40 ms
41,672 KB |
testcase_40 | AC | 70 ms
60,768 KB |
testcase_41 | AC | 38 ms
40,072 KB |
testcase_42 | AC | 77 ms
64,376 KB |
testcase_43 | AC | 44 ms
44,756 KB |
testcase_44 | AC | 31 ms
36,284 KB |
testcase_45 | AC | 30 ms
36,248 KB |
testcase_46 | AC | 31 ms
36,280 KB |
testcase_47 | AC | 30 ms
36,296 KB |
testcase_48 | AC | 30 ms
36,428 KB |
testcase_49 | AC | 30 ms
36,308 KB |
testcase_50 | AC | 30 ms
36,280 KB |
testcase_51 | AC | 31 ms
36,220 KB |
testcase_52 | AC | 30 ms
36,348 KB |
testcase_53 | AC | 31 ms
36,432 KB |
ソースコード
// #pragma GCC target("avx2") #pragma GCC optimize("O3") // #pragma GCC optimize("unroll-loops") #include<bits/stdc++.h> using namespace std; using P = pair<int, int>; const int M = 998244353; long long powmod(long long a, long long b) { long long res = 1; while (b) { if (b & 1) { res = res * a % M; } b >>= 1; a = a * a % M; } return res; } long long inv(long long a) { return powmod(a, M - 2); } int blim = 2048; vector<vector<long long>> binom_small(blim, vector<long long>(blim)); long long binom(int n, int r) { if (n - r < r) r = n - r; if (n < 0 || r < 0 || n < r) return 0; if (n < blim) return binom_small[n][r]; long long u = 1, d = 1; for (int i = 0; i < r; ++i) { u = u * (n - i) % M; d = d * (i + 1) % M; } return u * inv(d) % M; } int main() { cin.tie(0); ios::sync_with_stdio(0); for (int i = 0; i < blim; ++i) { binom_small[i][0] = binom_small[i][i] = 1; for (int j = 1; j < i; ++j) { binom_small[i][j] = (binom_small[i - 1][j - 1] + binom_small[i - 1][j]) % M; } } int n, m, k; cin >> n >> m >> k; vector<vector<long long>> dp1(m + 3, vector<long long>(m + 3)); dp1[0][1] = 1; for (int i = 1; i < m + 3; ++i) { vector<long long> imos(m + 4); for (int j = 0; j < m + 3; ++j) { imos[j + 1] = (imos[j] + dp1[i - 1][j]) % M; } for (int j = 1; j < m + 3; ++j) { dp1[i][j] = (imos[j] - imos[(j + 1) / 2] + M) % M; } } vector<long long> x(m + 3); for (int i = 0; i < m + 3 && i <= n; ++i) { long long b = binom(n, i); for (int j = 0; j + 1 < m + 3; ++j) { x[j] = (x[j] + dp1[i][j + 1] * b) % M; } } vector<vector<long long>> dp2(13, vector<long long>(m + 3)); dp2[0][1] = 1; for (int i = 0; i + 1 < 13; ++i) { for (int j = 1; j < m + 3; ++j) { for (int l = 2; j * l < m + 3; ++l) { dp2[i + 1][j * l] = (dp2[i + 1][j * l] + dp2[i][j] * x[l - 1]) % M; } } } long long ans = 0; for (int i = 0; i < 13; ++i) { long long p = binom(k, i); for (int j = 0; j <= m; ++j) { ans = (ans + dp2[i][j] * p) % M; } } cout << ans << '\n'; return 0; }