이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <vector>
#include <algorithm>
#include <iostream>
#include <set>
#include <cmath>
#include <map>
#include <random>
#include <cassert>
#include <ctime>
#include <cstdlib>
#include <queue>
#include <limits.h>
#pragma GCC target ("avx2")
#pragma GCC optimization ("O3")
#pragma GCC optimization ("unroll-loops")
using namespace std;
inline namespace _LineContainer {
bool _Line_Comp_State;
struct Line {
mutable long long k, m, p;
int id;
bool operator<(const Line& o) const {
if (o.k == this->k && o.m == this->m && o.p == this->p) return (o.id < id);
return _Line_Comp_State ? p < o.p : k < o.k;
}
};
struct LineContainer : multiset<Line> {
long long div(long long a, long long b) { return a / b - ((a ^ b) < 0 && a % b); }
bool isect(iterator x, iterator y) {
if (y == end()) {
x->p = LLONG_MAX;
return false;
}
if (x->k == y->k)
x->p = x->m > y->m ? LLONG_MAX : -LLONG_MAX;
else
x->p = div(y->m - x->m, x->k - y->k);
return x->p >= y->p;
}
void add(long long k, long long m, int id) {
auto z = insert({k, m, 0, id}), y = z++, x = y;
while (isect(y, z))
z = erase(z);
if (x != begin() && isect(--x, y))
isect(x, y = erase(y));
while ((y = x) != begin() && (--x)->p >= y->p)
isect(x, erase(y));
}
pair<long long, long long> query(long long x) {
assert(!empty());
_Line_Comp_State = true;
auto l = *lower_bound({0,0,x, -1});
_Line_Comp_State = false;
return {l.k * x + l.m, l.id};
}
};
}
int main () {
ios_base::sync_with_stdio(false);
cin.tie(NULL);
int n, k;
cin >> n >> k;
int64_t arr[n]; int64_t pref[n + 1]; pref[0] = 0;
for (int i = 0; i < n; i++) {
cin >> arr[i];
pref[i + 1] = pref[i] + arr[i];
}
vector<int64_t> dp_cur(n - 1), dp_prev(n - 1);
int prev[n - 1][k + 1];
for (int i = 0; i < n - 1; i++) {
dp_prev[i] = (pref[i + 1]) * (pref[n] - pref[i + 1]);
prev[i][1] = -1;
}
for (int j = 2; j <= k; j++) {
LineContainer cht;
dp_cur[0] = 0;
for (int i = 1; i < n - 1; i++) {
cht.add(pref[i], dp_prev[i - 1], i);
dp_cur[i] = max(cht.query(pref[i + 1] - pref[n]).first, 0ll) + pref[i + 1] * pref[n] - pref[i + 1] * pref[i + 1];
prev[i][j] = cht.query(pref[i + 1] - pref[n]).second - 1;
}
dp_prev = dp_cur;
}
int64_t myMax = 0;
for (int i = 0; i < n - 1; i++) {
myMax = max(myMax, dp_prev[i]);
}
cout << myMax << '\n';
for (int i = 0; i < n - 1; i++) {
if (myMax == dp_prev[i]) {
int mid = i;
int cntr = 0;
while (cntr != k) {
cout << mid + 1 << ' ';
mid = prev[mid][k - cntr];
cntr++;
}
return 0;
}
}
}
컴파일 시 표준 에러 (stderr) 메시지
sequence.cpp:14: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
14 | #pragma GCC optimization ("O3")
|
sequence.cpp:15: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
15 | #pragma GCC optimization ("unroll-loops")
|
# | 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... |