This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
/**
* author: hashman
* created:
**/
#include <bits/stdc++.h>
using namespace std;
string to_string(string s) {
return '"' + s + '"';
}
string to_string(const char* s) {
return to_string((string) s);
}
string to_string(bool b) {
return (b ? "true" : "false");
}
template <typename A, typename B>
string to_string(pair<A, B> p) {
return "(" + to_string(p.first) + ", " + to_string(p.second) + ")";
}
template <typename A>
string to_string(A v) {
bool first = true;
string res = "{";
for (const auto &x : v) {
if (!first) {
res += ", ";
}
first = false;
res += to_string(x);
}
res += "}";
return res;
}
void debug_out() {
cerr << endl;
}
template <typename Head, typename... Tail>
void debug_out(Head H, Tail... T) {
cerr << " " << to_string(H);
debug_out(T...);
}
#define debug(...) cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__);
int64_t inf = 1e17;
int main() {
ios::sync_with_stdio(false);
cin.tie(0);
int n, k;
cin >> n >> k;
vector<int> a(n + 1);
for (int i = 1; i <= n; i++) {
cin >> a[i];
}
vector<int> p(n + 1);
for (int i = 1; i <= n; i++) {
p[i] = p[i - 1] + a[i];
}
auto cost = [&](int l, int r) {
return ((int64_t) (p[r] - p[l - 1]) * (p[r] - p[l - 1]));
};
k++;
vector<int64_t> dp_cur(n + 1, inf), dp_before(n + 1, inf);
dp_before[0] = 0;
vector<vector<int>> pr(n + 1, vector<int>(k + 1));
auto compute = [&](auto&& self, int at, int l, int r, int bl, int br) -> void {
if (l > r) {
return;
}
int m = (l + r) / 2;
pair<int64_t, int> best = make_pair(inf, -1);
for (int i = max(bl - 100, 1); i <= min(m, br + 100); i++) {
best = min(best, make_pair(dp_before[i - 1] + cost(i, m), i));
}
if (best.second != -1) pr[m][at] = best.second - 1;
dp_cur[m] = best.first;
if (best.second != -1) self(self, at, l, m - 1, bl, best.second);
else self(self, at, l, m - 1, bl, br);
if (best.second != -1) self(self, at, m + 1, r, best.second, br);
else self(self, at, m + 1, r, bl, br);
};
for (int i = 1; i <= k; i++) {
compute(compute, i, 1, n, 1, n);
dp_before = dp_cur;
}
cout << ((int64_t) p[n] * p[n] - dp_before[n]) / 2 << '\n';
int cur = n;
vector<int> splits;
while (cur != 0) {
cur = pr[cur][k];
k--;
if (cur != 0) {
splits.push_back(cur);
}
}
reverse(splits.begin(), splits.end());
for (int split : splits) {
cout << split << " ";
}
cout << '\n';
}
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |