This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include <bits/stdc++.h>
#define pb push_back
#define mp make_pair
using namespace std;
typedef long long ll;
typedef long double ld;
const ll INF = numeric_limits<ll>::max();
struct Line {
ll k, m;
int id;
ll f(const ll &x) { return k * x + m; }
bool operator< (const Line &oth) {
return mp(k, m) < mp(oth.k, oth.m);
}
};
struct cvh {
int ptr = 0;
vector<Line> hull;
void reset() { ptr = 0; hull.clear(); }
inline int sz() { return hull.size(); }
bool bad(const Line &a, const Line &b, const Line &c) {
if(b.k == c.k && b.m <= c.m) return true;
return (ld) (a.m - b.m) / (ld) (b.k - a.k) >= (ld) (b.m - c.m) / (ld) (c.k - b.k);
}
void add(ll k, ll m, int id) {
Line newline = {k, m, id};
while(sz() >= 2 && bad(hull[sz() - 2], hull[sz() - 1], newline)) hull.pop_back();
hull.pb(newline);
}
Line query(ll x) {
if(!sz()) {
Line amogus = {-INF, -INF, -1};
return amogus;
}
if(ptr >= sz())
ptr = sz() - 1;
while(ptr + 1 < sz() && hull[ptr].f(x) <= hull[ptr + 1].f(x)) ptr++;
return hull[ptr];
}
} Cvh;
const int maxN = 1e5 + 5, maxK = 205;
int n, k;
int trace[maxK][maxN];
ll a[maxN], dp[maxN];
signed main() {
cin >> n >> k;
for(int i = 1; i <= n; i++) {
cin >> a[i];
a[i] += a[i - 1];
}
dp[0] = -INF;
for(int i = 2; i <= k + 1; i++) {
Cvh.reset();
for(int j = 1; j <= n; j++) {
Line best = Cvh.query(a[j]);
if(dp[j] != -INF)
Cvh.add(a[j], dp[j] - a[j] * a[j], j);
trace[i][j] = best.id;
if(best.id != -1)
dp[j] = best.f(a[j]);
else
dp[j] = -INF;
}
cout << i << '\n';
for(int j = 1; j <= n; j++)
cout << dp[j] << ' ';
cout << '\n';
}
cout << dp[n] << '\n';
vector<int> ans;
for(int i = k + 1, cur = n; i; i--) {
ans.pb(trace[i][cur]);
cur = trace[i][cur];
}
for(int i = (int)ans.size() - 1; i >= 0; i--)
if(ans[i]) cout << ans[i] << ' ';
}
# | 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... |