/*************************************
* author: marvinthang *
* created: 27.11.2023 15:25:09 *
*************************************/
#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define left ___left
#define right ___right
#define TIME (1.0 * clock() / CLOCKS_PER_SEC)
#define MASK(i) (1LL << (i))
#define BIT(x, i) ((x) >> (i) & 1)
#define __builtin_popcount __builtin_popcountll
#define ALL(v) (v).begin(), (v).end()
#define REP(i, n) for (int i = 0, _n = (n); i < _n; ++i)
#define REPD(i, n) for (int i = (n); i-- > 0; )
#define FOR(i, a, b) for (int i = (a), _b = (b); i < _b; ++i)
#define FORD(i, b, a) for (int i = (b), _a = (a); --i >= _a; )
#define FORE(i, a, b) for (int i = (a), _b = (b); i <= _b; ++i)
#define FORDE(i, b, a) for (int i = (b), _a = (a); i >= _a; --i)
#define scan_op(...) istream & operator >> (istream &in, __VA_ARGS__ &u)
#define print_op(...) ostream & operator << (ostream &out, const __VA_ARGS__ &u)
#ifdef LOCAL
#include "debug.h"
#else
#define file(name) if (fopen(name".inp", "r")) { freopen(name".inp", "r", stdin); freopen(name".out", "w", stdout); }
#define DB(...) 23
#define db(...) 23
#define debug(...) 23
#endif
template <class U, class V> scan_op(pair <U, V>) { return in >> u.first >> u.second; }
template <class T> scan_op(vector <T>) { for (size_t i = 0; i < u.size(); ++i) in >> u[i]; return in; }
template <class U, class V> print_op(pair <U, V>) { return out << '(' << u.first << ", " << u.second << ')'; }
template <size_t i, class T> ostream & print_tuple_utils(ostream &out, const T &tup) { if constexpr(i == tuple_size<T>::value) return out << ")"; else return print_tuple_utils<i + 1, T>(out << (i ? ", " : "(") << get<i>(tup), tup); }
template <class ...U> print_op(tuple<U...>) { return print_tuple_utils<0, tuple<U...>>(out, u); }
template <class Con, class = decltype(begin(declval<Con>()))> typename enable_if <!is_same<Con, string>::value, ostream&>::type operator << (ostream &out, const Con &con) { out << '{'; for (__typeof(con.begin()) it = con.begin(); it != con.end(); ++it) out << (it == con.begin() ? "" : ", ") << *it; return out << '}'; }
template <class A, class B> bool minimize(A &a, B b) { if (a > b) { a = b; return true; } return false; }
template <class A, class B> bool maximize(A &a, B b) { if (a < b) { a = b; return true; } return false; }
const long long INF = 1e18;
struct Line {
long long a, b;
int pos;
Line(long long a = 0, long long b = INF, int pos = -1): a(a), b(b), pos(pos) {}
long long eval(long long x) { return a * x + b; }
};
struct ConvexHullTrick {
vector <Line> lines; int ptr;
ConvexHullTrick(): ptr(0) {}
bool bad(Line a, Line b, Line c) { return (long double) (c.b - a.b) / (a.a - c.a) < (long double) (b.b - a.b) / (a.a - b.a); }
void addLine(long long a, long long b, int c) {
Line l(a, b, c);
if (!lines.empty() && lines.back().a == a) {
if (lines.back().b <= b) return;
lines.pop_back();
}
while (lines.size() >= 2 && bad(lines.end()[-2], lines.back(), l)) lines.pop_back();
lines.push_back(l);
}
Line getMin(long long x) {
if (ptr >= (int) lines.size()) ptr = (int) lines.size() - 1;
while (ptr < (int) lines.size() - 1 && lines[ptr].eval(x) > lines[ptr + 1].eval(x)) ++ptr;
return lines[ptr];
}
};
using CHT = ConvexHullTrick;
// end of template
void process(void) {
int n, k; cin >> n >> k; ++k;
vector <int> a(n), pref(n + 1); cin >> a;
partial_sum(ALL(a), pref.begin() + 1);
vector dp(2, vector(n + 1, INF));
vector trace(n + 1, vector(k + 1, 0));
dp[0][0] = 0;
CHT cht;
REP(c, k) {
cht.lines.clear(); cht.ptr = 0;
int cur = c & 1;
FORE(i, c, n) {
if (i > c) {
Line l = cht.getMin(pref[i]);
dp[!cur][i] = l.eval(pref[i]) + 1LL * pref[i] * pref[i];
trace[i][c + 1] = l.pos;
}
if (dp[cur][i] < INF) cht.addLine(-2 * pref[i], dp[cur][i] + 1LL * pref[i] * pref[i], i);
}
}
cout << (1LL * pref[n] * pref[n] - dp[k & 1][n]) / 2 << '\n';
vector <int> res;
int i = n;
FORDE(c, k, 2) {
i = trace[i][c];
res.push_back(i);
}
reverse(ALL(res));
for (int x: res) cout << x << ' ';
}
int main(void) {
ios_base::sync_with_stdio(false); cin.tie(nullptr); // cout.tie(nullptr);
file("sequence");
// int t; cin >> t; while (t--)
process();
// cerr << "Time elapsed: " << TIME << " s.\n";
return (0^0);
}
Compilation message
sequence.cpp: In function 'int main()':
sequence.cpp:30:61: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
30 | #define file(name) if (fopen(name".inp", "r")) { freopen(name".inp", "r", stdin); freopen(name".out", "w", stdout); }
| ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~
sequence.cpp:110:2: note: in expansion of macro 'file'
110 | file("sequence");
| ^~~~
sequence.cpp:30:94: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
30 | #define file(name) if (fopen(name".inp", "r")) { freopen(name".inp", "r", stdin); freopen(name".out", "w", stdout); }
| ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~
sequence.cpp:110:2: note: in expansion of macro 'file'
110 | file("sequence");
| ^~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 108 == 108 |
2 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 999 == 999 |
3 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 0 == 0 |
4 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 1542524 == 1542524 |
5 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 4500000000 == 4500000000 |
6 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 1 == 1 |
7 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 1 == 1 |
8 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 1 == 1 |
9 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 100400096 == 100400096 |
10 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 900320000 == 900320000 |
11 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 3698080248 == 3698080248 |
12 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 3200320000 == 3200320000 |
13 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 140072 == 140072 |
14 |
Correct |
1 ms |
344 KB |
contestant found the optimal answer: 376041456 == 376041456 |
15 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 805 == 805 |
16 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 900189994 == 900189994 |
17 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 999919994 == 999919994 |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
344 KB |
contestant found the optimal answer: 1093956 == 1093956 |
2 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 302460000 == 302460000 |
3 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 122453454361 == 122453454361 |
4 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 93663683509 == 93663683509 |
5 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 1005304678 == 1005304678 |
6 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 933702 == 933702 |
7 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 25082842857 == 25082842857 |
8 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 687136 == 687136 |
9 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 27295930079 == 27295930079 |
10 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 29000419931 == 29000419931 |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
344 KB |
contestant found the optimal answer: 610590000 == 610590000 |
2 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 311760000 == 311760000 |
3 |
Correct |
1 ms |
604 KB |
contestant found the optimal answer: 1989216017013 == 1989216017013 |
4 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 1499437552673 == 1499437552673 |
5 |
Correct |
1 ms |
348 KB |
contestant found the optimal answer: 1019625819 == 1019625819 |
6 |
Correct |
1 ms |
600 KB |
contestant found the optimal answer: 107630884 == 107630884 |
7 |
Correct |
1 ms |
604 KB |
contestant found the optimal answer: 475357671774 == 475357671774 |
8 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 193556962 == 193556962 |
9 |
Correct |
1 ms |
348 KB |
contestant found the optimal answer: 482389919803 == 482389919803 |
10 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 490686959791 == 490686959791 |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 21503404 == 21503404 |
2 |
Correct |
0 ms |
348 KB |
contestant found the optimal answer: 140412195 == 140412195 |
3 |
Correct |
5 ms |
1116 KB |
contestant found the optimal answer: 49729674225461 == 49729674225461 |
4 |
Correct |
1 ms |
348 KB |
contestant found the optimal answer: 37485571387523 == 37485571387523 |
5 |
Correct |
4 ms |
1120 KB |
contestant found the optimal answer: 679388326 == 679388326 |
6 |
Correct |
4 ms |
1316 KB |
contestant found the optimal answer: 4699030287 == 4699030287 |
7 |
Correct |
4 ms |
1348 KB |
contestant found the optimal answer: 12418819758185 == 12418819758185 |
8 |
Correct |
3 ms |
1116 KB |
contestant found the optimal answer: 31093317350 == 31093317350 |
9 |
Correct |
1 ms |
600 KB |
contestant found the optimal answer: 12194625429236 == 12194625429236 |
10 |
Correct |
2 ms |
860 KB |
contestant found the optimal answer: 12345131038664 == 12345131038664 |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
2 ms |
1560 KB |
contestant found the optimal answer: 1818678304 == 1818678304 |
2 |
Correct |
2 ms |
1564 KB |
contestant found the optimal answer: 1326260195 == 1326260195 |
3 |
Correct |
31 ms |
9332 KB |
contestant found the optimal answer: 4973126687469639 == 4973126687469639 |
4 |
Correct |
2 ms |
1560 KB |
contestant found the optimal answer: 3748491676694116 == 3748491676694116 |
5 |
Correct |
26 ms |
5980 KB |
contestant found the optimal answer: 1085432199 == 1085432199 |
6 |
Correct |
30 ms |
6748 KB |
contestant found the optimal answer: 514790755404 == 514790755404 |
7 |
Correct |
32 ms |
7452 KB |
contestant found the optimal answer: 1256105310476641 == 1256105310476641 |
8 |
Correct |
23 ms |
6428 KB |
contestant found the optimal answer: 3099592898816 == 3099592898816 |
9 |
Correct |
32 ms |
6940 KB |
contestant found the optimal answer: 1241131419367412 == 1241131419367412 |
10 |
Correct |
31 ms |
8476 KB |
contestant found the optimal answer: 1243084101967798 == 1243084101967798 |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
16 ms |
12496 KB |
contestant found the optimal answer: 19795776960 == 19795776960 |
2 |
Correct |
17 ms |
11472 KB |
contestant found the optimal answer: 19874432173 == 19874432173 |
3 |
Correct |
446 ms |
88656 KB |
contestant found the optimal answer: 497313449256899208 == 497313449256899208 |
4 |
Correct |
17 ms |
13008 KB |
contestant found the optimal answer: 374850090734572421 == 374850090734572421 |
5 |
Correct |
485 ms |
88352 KB |
contestant found the optimal answer: 36183271951 == 36183271951 |
6 |
Correct |
326 ms |
66448 KB |
contestant found the optimal answer: 51629847150471 == 51629847150471 |
7 |
Correct |
378 ms |
71100 KB |
contestant found the optimal answer: 124074747024496432 == 124074747024496432 |
8 |
Correct |
333 ms |
60116 KB |
contestant found the optimal answer: 309959349080800 == 309959349080800 |
9 |
Correct |
312 ms |
65984 KB |
contestant found the optimal answer: 124113525649823701 == 124113525649823701 |
10 |
Correct |
408 ms |
80988 KB |
contestant found the optimal answer: 124309619349406845 == 124309619349406845 |