//#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;
line(int k, int b) {
this -> k = k;
this -> b = b;
}
int val(int x) {
return k * x + b;
}
db intrsc(const line & l) {
return (db)(l.b - b) / (db)(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 = line(x, 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]};
int mx = -1;
for (int i = 0; i < pnt.size(); i++) {
if (i > 0) {
assert(pnt[i] >= pnt[i - 1]);
}
if (mx == -1 || h[i].val(x) > h[mx].val(x)) {
mx = i;
}
}
return {h[mx].val(x), ind[mx]};
}
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();
}
}
Compilation message
sequence.cpp: In function 'pii getmax(long long int)':
sequence.cpp:76:23: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (int i = 0; i < pnt.size(); i++) {
~~^~~~~~~~~~~~
sequence.cpp: In function 'int main()':
sequence.cpp:135:12: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)', declared with attribute warn_unused_result [-Wunused-result]
freopen("input.txt", "r", stdin);
~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
4 ms |
376 KB |
Integer -1 violates the range [1, 6] |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
4 ms |
376 KB |
Integer -1 violates the range [1, 49] |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
4 ms |
380 KB |
Integer -1 violates the range [1, 199] |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
4 ms |
376 KB |
Integer -1 violates the range [1, 999] |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
4 ms |
376 KB |
Integer -1 violates the range [1, 9999] |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
4 ms |
380 KB |
Integer -1 violates the range [1, 99999] |
2 |
Halted |
0 ms |
0 KB |
- |