#include <bits/stdc++.h>
using i64 = long long;
struct line {
i64 a, b;
line(i64 _a, i64 _b) : a(_a), b(_b) {}
long double intersect(line rhs) {
return (long double)(b - rhs.b) / (rhs.a - a);
}
i64 eval(i64 x) {
return a * x + b;
}
};
namespace std {
std::string to_string(line l) {
return '(' + l.a + "x + " + l.b + ')';
}
};
#ifdef DEBUG
#include "debug.h"
#else
#define debug(...) void(23)
#endif
constexpr int maxN = 1E5 + 5;
constexpr int maxK = 200 + 5;
int A[maxN], go[maxK][maxN];
i64 dp[maxK][maxN], p[maxN];
signed main() {
std::ios::sync_with_stdio(false);
std::cin.tie(nullptr);
int N, K;
std::cin >> N >> K;
for(int i = 1; i <= N; i++) {
std::cin >> A[i];
}
for(int i = 1; i <= N; i++) {
p[i] = p[i - 1] + A[i];
}
std::deque<std::pair<line, int>> dq;
auto size = [&]() -> int { return (int) dq.size(); };
auto Add = [&](line l, int idx) -> void {
while(size() >= 2 && l.intersect(dq[0].first) >= dq[0].first.intersect(dq[1].first)) {
dq.pop_front();
}
dq.emplace_back(l, idx);
};
auto Query = [&](i64 x) -> std::pair<i64, int> {
while(size() >= 2 && dq[size() - 1].first.eval(x) <= dq[size() - 2].first.eval(x)) {
dq.pop_back();
}
return {dq[size() - 1].first.eval(x), dq[size() - 1].second};
};
for(int k = 1; k <= K; k++) {
dq.clear();
Add({0, 0}, 0);
for(int i = 1; i <= N; i++) {
// for(int j = 0; j < i; j++) {
// if(dp[k - 1][j] + 1LL * (p[i] - p[j]) * (p[N] - p[i]) >= dp[k][i]) {
// go[k][i] = j;
// dp[k][i] = dp[k - 1][j] + 1LL * (p[i] - p[j]) * (p[N] - p[i]);
// }
// }
auto [v, j] = Query(p[N] - p[i]);
dp[k][i] = v + 1LL * p[i] * (p[N] - p[i]);
go[k][i] = j;
Add(line{-p[i], dp[k - 1][i]}, i);
}
}
int ans = 1;
for(int i = 1; i < N; i++) {
if(dp[K][i] >= dp[K][ans]) {
ans = i;
}
}
std::cout << dp[K][ans] << '\n';
std::vector<int> v;
for(int k = K; k >= 1; k--) {
v.emplace_back(ans);
ans = go[k][ans];
}
for(int i = K - 1; i >= 0; i--) {
std::cout << v[i] << " \n"[i == 0];
}
return 0;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
8540 KB |
contestant found the optimal answer: 108 == 108 |
2 |
Correct |
1 ms |
6492 KB |
contestant found the optimal answer: 999 == 999 |
3 |
Correct |
1 ms |
4444 KB |
contestant found the optimal answer: 0 == 0 |
4 |
Incorrect |
1 ms |
10712 KB |
contestant didn't find the optimal answer: 1542522 < 1542524 |
5 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
1 ms |
8536 KB |
contestant didn't find the optimal answer: 1092528 < 1093956 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
8540 KB |
contestant found the optimal answer: 610590000 == 610590000 |
2 |
Correct |
1 ms |
8536 KB |
contestant found the optimal answer: 311760000 == 311760000 |
3 |
Correct |
35 ms |
86096 KB |
contestant found the optimal answer: 1989216017013 == 1989216017013 |
4 |
Incorrect |
1 ms |
8536 KB |
contestant didn't find the optimal answer: 1066891648920 < 1499437552673 |
5 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
1 ms |
8540 KB |
contestant didn't find the optimal answer: 20166072 < 21503404 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
2 ms |
8536 KB |
contestant didn't find the optimal answer: 1794250000 < 1818678304 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
12 ms |
9048 KB |
contestant didn't find the optimal answer: 13198103640 < 19795776960 |
2 |
Halted |
0 ms |
0 KB |
- |