이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <iostream>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <algorithm>
#include <string>
#include <cmath>
#include <cstdio>
#include <iomanip>
#include <fstream>
#include <cassert>
#include <cstring>
#include <unordered_set>
#include <unordered_map>
#include <numeric>
#include <ctime>
#include <bitset>
#include <complex>
#include <chrono>
#include <random>
#include <functional>
using namespace std;
#define int long long
const int INF = 1e18 + 239;
const int N = 1e5 + 7;
const int K = 203;
int dp[K][N];
int wr[K][N];
int a[N];
int ps[N];
vector<int> out;
void rec(int k, int n) {
if (k == 0) return;
rec(k - 1, wr[k][n]);
out.push_back(n - wr[k][n]);
}
namespace CHT {
typedef long double ld;
struct Line {
int k, b;
ld start;
int id;
Line(int k_, int b_, ld start_, int id_) {
k = k_;
b = b_;
start = start_;
id = id_;
}
ld eval_ld(ld x) {
return k * x + b;
}
};
int pnt = 0;
vector<Line> st;
void init() {
pnt = 0;
st.clear();
}
void add(int k, int b, int id) {
Line cur(k, b, -1, id);
while (!st.empty()) {
Line bc = st.back();
if (cur.eval_ld(bc.start) < bc.eval_ld(bc.start)) {
st.pop_back();
} else {
break;
}
}
if (st.empty()) {
cur.start = 0;
st.push_back(cur);
} else {
if (st.back().k == cur.k) return;
ld start = (st.back().b - cur.b) / (cur.k - st.back().k);
cur.start = start;
st.push_back(cur);
}
}
int get(int x) {
if (st.empty()) return -1;
pnt = min(pnt, (int)st.size() - 1);
while (pnt + 1 < (int)st.size() && st[pnt + 1].start < x) {
pnt++;
}
return st[pnt].id;
}
}
signed main() {
ios_base::sync_with_stdio(false);
cin.tie(0);
int n, k;
cin >> n >> k;
k++;
int s = 0;
for (int i = 0; i < n; i++) {
cin >> a[i];
s += a[i] * a[i];
ps[i + 1] = ps[i] + a[i];
}
fill(dp[0], dp[0] + N, INF);
dp[0][0] = 0;
for (int i = 1; i <= k; i++) {
dp[i][0] = INF;
CHT::init();
for (int j = 1; j <= n; j++) {
dp[i][j] = INF;
int add = j - 1;
if (dp[i - 1][add] != INF) {
CHT::add(-ps[add], dp[i - 1][add] + ps[add] * ps[add], add);
}
int p = CHT::get(2 * ps[j]);
if (p != -1) {
dp[i][j] = dp[i - 1][p] + ps[p] * ps[p] - 2 * ps[p] * ps[j] + ps[j] * ps[j];
wr[i][j] = p;
}
}
}
cout << (ps[n] * ps[n] - dp[k][n]) / 2 << endl;
rec(k, n);
{
int c = 0;
for (int i = 0; i + 1 < k; i++) {
c += out[i];
cout << c << ' ';
}
cout << endl;
}
}
# | 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... |