제출 #915870

#제출 시각아이디문제언어결과실행 시간메모리
915870Pekiban수열 (APIO14_sequence)C++17
100 / 100
455 ms82028 KiB
#include <bits/stdc++.h>

using namespace std;

const int mxK = 202, mxN = 1e5+2;
bool vis[mxN];
int to[mxK][mxN], dq[mxN], l = 0, r = 0;
long long dpb[mxN], dp[mxN], pref[mxN];
long long sq(long long a) {
    return a*a;
}
double cord(int x, int y) {
    return (double)(dpb[x] - sq(pref[x]) - dpb[y] + sq(pref[y]))/(pref[y] - pref[x]);// (d-b)/(a-c)
}
int main() {
    ios::sync_with_stdio(0);
    cin.tie(0);
    int n, k;
    cin >> n >> k;
    for (int i = 1; i <= n; ++i) {
        int x;
        cin >> x;
        pref[i] = pref[i-1] + x;
    }
    for (int i = 1; i <= k; ++i) {
        dq[0] = r = l = 0;
        for (int j = 1; j <= n; ++j) {
            while (r - l > 0 && pref[j] != pref[j-1] && cord(dq[l], dq[l+1]) <= pref[j])   l++;
            dp[j] = dpb[dq[l]] - sq(pref[dq[l]]) + pref[dq[l]] * pref[j];
            to[i][j] = dq[l];
            if (pref[j] != pref[dq[r-1]]) {
                while (r - l > 0 && pref[j] != pref[dq[r-1]] && cord(dq[r-1], j) <= cord(dq[r], dq[r-1]))  r--;
                dq[++r] = j;
            } 
        }
        for (int j = 1; j <= n; ++j) {
            dpb[j] = dp[j];
        }
    }
    cout << dp[n] << '\n';
    int idx = n;
    for (int i = k; i >= 1; --i) {
        if (!to[i][idx]) {
            for (int j = n-1; j >= 1 && i >= 1; --j){
                if (!vis[j]) {
                    --i;
                    cout << j << ' ';
                }
            }
            break;
        }
        cout << to[i][idx] << ' ';
        idx = to[i][idx];
        vis[idx] = 1;
    }
}
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...