結果
| 問題 |
No.2215 Slide Subset Sum
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-02-29 17:49:02 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 1,513 ms / 3,000 ms |
| コード長 | 1,917 bytes |
| コンパイル時間 | 4,716 ms |
| コンパイル使用メモリ | 277,544 KB |
| 実行使用メモリ | 98,380 KB |
| 最終ジャッジ日時 | 2024-09-29 12:47:26 |
| 合計ジャッジ時間 | 20,744 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 45 |
ソースコード
#pragma GCC optimize("Ofast,unroll-loops")
#pragma GCC target("avx2,popcnt")
#include <bits/stdc++.h>
using namespace std;
const int P = 998244353;
int main() {
ios::sync_with_stdio(false);
cin.tie(nullptr);
int N, M, K;
cin >> N >> M >> K;
vector<int> A(N);
for (int i = 0; i < N; i++) {
cin >> A[i];
}
int Q = N - M + 1;
vector<int> ans(Q);
vector<tuple<int, int, int>> querys;
for (int i = 0; i < Q; i++) {
querys.emplace_back(i, i + M, i);
}
vector<vector<int>> dp(N + 1, vector<int>(K));
auto add = [&](vector<int> &f, int v) {
vector<int> g = f;
for (int i = 0; i < K; i++) {
int j = (i + v < K ? i + v : i + v - K);
g[j] += f[i];
if (g[j] >= P) {
g[j] -= P;
}
}
return g;
};
auto dfs = [&](auto &dfs, int L, int R,
vector<tuple<int, int, int>> &query) {
vector<tuple<int, int, int>> qL, qM, qR;
int M = (L + R) / 2, lb = M, rb = M;
for (auto &[l, r, i] : query) {
if (l <= M && M <= r) {
lb = min(lb, L);
rb = max(rb, R);
qM.emplace_back(l, r, i);
} else if (r < M) {
qL.emplace_back(l, r, i);
} else {
qR.emplace_back(l, r, i);
}
}
for (int j = 0; j < K; j++) {
dp[M][j] = (j == 0);
}
for (int i = M - 1; i >= lb; i--) {
dp[i] = add(dp[i + 1], A[i]);
}
for (int i = M; i < rb; i++) {
dp[i + 1] = add(dp[i], A[i]);
}
for (auto [l, r, i] : qM) {
int res = P - 1;
for (int j = 0; j < K; j++) {
int k = (j == 0 ? 0 : K - j);
res += 1ll * dp[l][j] * dp[r][k] % P;
if (res >= P) {
res -= P;
}
}
ans[i] = res;
}
if (L + 1 == R) {
return;
} else {
dfs(dfs, L, M, qL);
dfs(dfs, M, R, qR);
}
};
dfs(dfs, 0, N, querys);
for (int i = 0; i < Q; i++) {
cout << ans[i] << "\n";
}
}