| # | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
|---|---|---|---|---|---|---|---|
| 1087501 | Thunnus | 수열 (APIO14_sequence) | C++17 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include<bits/stdc++.h>
using namespace std;
using i64 = long long;
//#define int i64
#define vi vector<int>
#define vvi vector<vi>
#define vb vector<bool>
#define pii pair<int, int>
#define fi first
#define se second
#define sz(x) (int)(x).size()
struct Line{
int m, c, i;
Line(i64 _m, i64 _c, int _i) : m(_m), c(_c), i(_i) {}
inline int calc(int x){
return m * x + c;
}
inline long double x_int(Line &other){
return (long double)(c - other.c) / (other.m - m);
}
bool operator ==(Line &other){
return (m == other.m && c == other.c);
}
};
signed main(){
ios_base::sync_with_stdio(false); cin.tie(0);
int n, k;
cin >> n >> k;
vi a(n + 1);
for(int i = 1; i <= n; i++){
cin >> a[i];
a[i] += a[i - 1];
}
vector<vector<i64>> dp(2, vector<i64>(n + 1));
vvi from(n + 1, vi(k + 2));
for(int j = 1; j <= k + 1; j++){
deque<Line> dq;
dq.emplace_front(0, 0, 0);
for(int i = 1; i <= n; i++){
int ps = a[n] - a[i];
while(sz(dq) >= 2 && dq.back().x_int(dq[sz(dq) - 2]) > ps){
dq.pop_back();
}
dp[1][i] = dq.back().calc(ps) + a[i] * ps;
from[i][j] = dq.back().i;
Line cur = {-a[i], dp[0][i], i};
while(sz(dq) >= 2 && cur.x_int(dq[0]) >< dq[1].x_int(dq[0])){
dq.pop_front();
}
if(!dq.empty() && cur == dq[0]){
dq.pop_front();
}
dq.emplace_front(cur);
}
dp[0] = dp[1];
}
cout << dp[0][n] << "\n";
vi ind;
for(int i = k + 1, idx = n; i >= 1; i--){
ind.emplace_back(idx);
idx = from[idx][i];
}
for(int i = sz(ind) - 1; i >= 1; i--){
cout << ind[i] << " ";
}
cout << "\n";
return 0;
}
// dp[i][j] = max(dp[z - 1][j - 1] + (a[i] - a[z - 1]) * (a[n] - a[i]))
// dp[i][j] = max(-a[z - 1] * (a[n] - a[i]) + dp[z - 1][j - 1]) + a[i] * (a[n] - a[i])