//#pragma GCC optimize("Ofast")
//#pragma GCC target("sse2,avx")
#include <bits/stdc++.h>
#define ll long long
#define db long double
#define pb push_back
#define pf push_front
#define ppb pop_back
#define ppf pop_front
#define fi first
#define se second
#define mp make_pair
#define endl "\n"
#define int long long
using namespace std;
void dout() {
cerr << endl;
}
template <typename Head, typename... Tail>
void dout(Head H, Tail... T) {
cerr << H << ' ';
dout(T...);
}
mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
typedef pair <int, int> pii;
const int N = 1e5 + 3, K = 203, inf = 1e18;
int n, k, a[N], dp[K][N], p[N], last[K][N];
struct line {
int k = 0, b = 0;
int val(int x) {
return k * x + b;
}
db intrsc(const line & l) {
return (db)(l.b - b) / (k - l.k);
}
bool del(line & l1, line & l2) {
if (k == l1.k) {
return b > l1.b;
}
return intrsc(l2) > l1.intrsc(l2);
}
};
vector <line> h;
vector <int> ind, ans, pnt;
void upd(int x, int y, int pos) {
line tmp;
tmp.k = x, tmp.b = y;
while (h.size() > 1 && tmp.del(h[h.size() - 1], h[h.size() - 2])) {
h.ppb();
pnt.ppb();
ind.ppb();
}
if (h.size() > 0) {
pnt.pb((int)(tmp.intrsc(h.back())));
}
h.pb(tmp);
ind.pb(pos);
}
pii getmax(int x) {
if (h.empty()) {
return {0, 0};
}
int tmp = upper_bound(pnt.begin(), pnt.end(), x) - pnt.begin() - 1;
return {h[tmp].val(x), ind[tmp]};
}
void getans(int x, int y) {
if (!y) {
return;
}
ans.pb(y);
getans(x - 1, last[x][y]);
}
void solve(int tc) {
cin >> n >> k;
for (int i = 1; i <= n; i++) {
cin >> a[i];
p[i] = p[i - 1] + a[i];
}
for (int i = 1; i <= k; i++) {
h = {}, pnt = {-inf}, ind = {};
for (int j = 1; j < i; j++) {
if (i > 1) {
upd(p[j], dp[i - 1][j], j);
}
}
for (int j = i; j < n; j++) {
int x = p[n] - p[j];
pii z = getmax(-x);
dp[i][j] = p[j] * x + z.fi;
last[i][j] = z.se;
if (i > 1) {
upd(p[j], dp[i - 1][j], j);
}
}
}
int res = -1;
for (int i = 1; i < n; i++) {
if (res == -1 || dp[k][res] < dp[k][i]) {
res = i;
}
}
getans(k, res);
cout << dp[k][res] << endl;
reverse(ans.begin(), ans.end());
for (auto i : ans) {
cout << i << " ";
}
}
signed main() {
ios_base::sync_with_stdio(false);
cin.tie(NULL);
// freopen("input.txt", "r", stdin);
// freopen("output.txt", "w", stdout);
int tc = 1;
// cin >> tc;
for (int i = 0; i < tc; i++) {
solve(i);
// cleanup();
}
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
380 KB |
contestant found the optimal answer: 108 == 108 |
2 |
Incorrect |
3 ms |
376 KB |
contestant didn't find the optimal answer: 908 < 999 |
3 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
2 ms |
376 KB |
contestant didn't find the optimal answer: 1093127 < 1093956 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
376 KB |
contestant found the optimal answer: 610590000 == 610590000 |
2 |
Correct |
2 ms |
376 KB |
contestant found the optimal answer: 311760000 == 311760000 |
3 |
Incorrect |
6 ms |
2936 KB |
Unexpected end of file - int32 expected |
4 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
3 ms |
504 KB |
contestant didn't find the optimal answer: 21343000 < 21503404 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
8 ms |
1144 KB |
contestant didn't find the optimal answer: 1478399266 < 1818678304 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
70 ms |
8316 KB |
contestant didn't find the optimal answer: 19795232584 < 19795776960 |
2 |
Halted |
0 ms |
0 KB |
- |