이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
using namespace std;
#define int long long
typedef pair<int, int32_t> pii;
#define N 100001
#define K 201
int ps[N];
int32_t dp[N][K];
struct Frac {
int n, d;
Frac(int num, int den) : n{num}, d{den} { if (d < 0) n = -n, d = -d; }
bool operator < (const Frac &o) const { return __int128_t(n)*__int128_t(o.d) < __int128_t(o.n)*__int128_t(d); }
};
struct Line {
int m, c, i;
Frac isect(const Line &o) const { return {c - o.c, o.m - m}; }
int operator () (const int x) const { return m*x + c; }
};
struct CHT {
deque<Line> hull;
void add(Line l){
while (hull.size() >= 2 && hull[1].isect(l) < hull[1].isect(hull[0]))
hull.pop_front();
hull.push_front(l);
}
pii query(int x){
while (hull.size() >= 2 && hull.back()(x) <= hull[hull.size()-2](x))
hull.pop_back();
return {hull.back()(x), hull.back().i};
}
};
int32_t main(){
ios_base::sync_with_stdio(false);
cin.tie(NULL);
int n, k; cin >> n >> k;
for (int i = 1; i <= n; i++){
cin >> ps[i];
ps[i] += ps[i-1];
}
CHT prv;
for (int i = 1; i <= n; i++)
prv.add({ps[i], -ps[i]*ps[i], i});
int ans = 0;
for (int l = 1; l <= k; l++){
CHT curr;
for (int i = 1; i <= n; i++){
auto [v, j] = prv.query(ps[i]);
dp[i][l] = j;
ans = v;
curr.add({ps[i], v - ps[i]*ps[i], i});
}
prv = curr;
}
cout << ans << '\n';
int curr = n;
for (int i = k; i > 0; i--){
cout << dp[curr][i] << ' ';
curr = dp[curr][i];
}
cout << '\n';
return 0;
}
# | 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... |