Submission #154864

#	Submission time	Handle	Problem	Language	Result	Execution time	Memory
154864	2019-09-25T12:04:07 Z	srvlt	Split the sequence (APIO14_sequence)	C++14	0 / 100	78 ms	8404 KB

sequence

//#pragma GCC optimize("Ofast")
//#pragma GCC target("sse2,avx")
#include <bits/stdc++.h>
#define ll long long
#define db long double
#define pb push_back
#define pf push_front
#define ppb pop_back
#define ppf pop_front
#define fi first
#define se second
#define mp make_pair
#define endl "\n"
#define int long long
using namespace std;

void dout() {
    cerr << endl;
}

template <typename Head, typename... Tail>
void dout(Head H, Tail... T) {
    cerr << H << ' ';
    dout(T...);
}

mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
typedef pair <int, int> pii;
const int N = 1e5 + 3, K = 203, inf = 1e18;
int n, k, a[N], dp[K][N], p[N], last[K][N];

struct line {
    int k = 0, b = 0;
    int val(int x) {
        return k * x + b;
    }
    db intrsc(const line & l) {
        return (db)(l.b - b) / (k - l.k);
    }
    bool del(line & l1, line & l2) {
        if (k == l1.k) {
            return b > l1.b;
        }
        return intrsc(l2) > l1.intrsc(l2);
    }
};

vector <line> h;
vector <int> ind, ans, pnt;

void upd(int x, int y, int pos) {
    line tmp;
    tmp.k = x, tmp.b = y;
    while (h.size() > 1 && tmp.del(h[h.size() - 1], h[h.size() - 2])) {
        h.ppb();
        pnt.ppb();
        ind.ppb();
    }
    if (h.size() > 0) {
        pnt.pb((int)ceil(tmp.intrsc(h.back())));
    }
    h.pb(tmp);
    ind.pb(pos);
}

pii getmax(int x) {
    if (h.empty()) {
        return {0, 0};
    }
    int tmp = upper_bound(pnt.begin(), pnt.end(), x) - pnt.begin() - 1;
    return {h[tmp].val(x), ind[tmp]};
}

void getans(int x, int y) {
    if (!y) {
        return;
    }
    ans.pb(y);
    getans(x - 1, last[x][y]);
}

void solve(int tc) {
    cin >> n >> k;
    for (int i = 1; i <= n; i++) {
        cin >> a[i];
        p[i] = p[i - 1] + a[i];
    }
    for (int i = 1; i <= k; i++) {
        h = {}, pnt = {-inf}, ind = {};
        for (int j = 1; j < i; j++) {
            if (i > 1) {
                upd(p[j], dp[i - 1][j], j);
            }
        }
        for (int j = i; j < n; j++) {
            int x = p[n] - p[j];
            pii z = getmax(-x);
            dp[i][j] = p[j] * x + z.fi;
            last[i][j] = z.se;
            if (i > 1) {
                upd(p[j], dp[i - 1][j], j);
            }
        }
    }
    int res = -1;
    for (int i = 1; i < n; i++) {
        if (res == -1 || dp[k][res] < dp[k][i]) {
            res = i;
        }
    }
    getans(k, res);
    cout << dp[k][res] << endl;
    reverse(ans.begin(), ans.end());
    for (auto i : ans) {
        cout << i << " ";
    }
}

signed main() {
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
//    freopen("input.txt", "r", stdin);
//    freopen("output.txt", "w", stdout);
    int tc = 1;
//    cin >> tc;
    for (int i = 0; i < tc; i++) {
        solve(i);
//        cleanup();
    }
}

Subtask #1
0 / 11.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	2 ms	376 KB	contestant found the optimal answer: 108 == 108
2	Incorrect	2 ms	376 KB	contestant didn't find the optimal answer: 908 < 999
3	Halted	0 ms	0 KB	-

Subtask #2
0 / 11.0

#	Verdict	Execution time	Memory	Grader output
1	Incorrect	2 ms	376 KB	contestant didn't find the optimal answer: 1093127 < 1093956
2	Halted	0 ms	0 KB	-

Subtask #3
0 / 11.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	2 ms	376 KB	contestant found the optimal answer: 610590000 == 610590000
2	Correct	2 ms	376 KB	contestant found the optimal answer: 311760000 == 311760000
3	Incorrect	5 ms	2812 KB	Unexpected end of file - int32 expected
4	Halted	0 ms	0 KB	-