Submission #154873

# Submission time Handle Problem Language Result Execution time Memory
154873 2019-09-25T12:34:09 Z srvlt Split the sequence (APIO14_sequence) C++14
0 / 100
4 ms 380 KB
//#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;
    line(int k, int b) {
        this -> k = k;
        this -> b = b;
    }
    int val(int x) {
        return k * x + b;
    }
    db intrsc(const line & l) {
        return (db)(l.b - b) / (db)(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 = line(x, 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)(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]};
    int mx = -1;
    for (int i = 0; i < pnt.size(); i++) {
        if (i > 0) {
            assert(pnt[i] >= pnt[i - 1]);
        }
        if (mx == -1 || h[i].val(x) > h[mx].val(x)) {
            mx = i;
        }
    }
    return {h[mx].val(x), ind[mx]};
}

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();
    }
}

Compilation message

sequence.cpp: In function 'pii getmax(long long int)':
sequence.cpp:76:23: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
     for (int i = 0; i < pnt.size(); i++) {
                     ~~^~~~~~~~~~~~
sequence.cpp: In function 'int main()':
sequence.cpp:135:12: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)', declared with attribute warn_unused_result [-Wunused-result]
     freopen("input.txt", "r", stdin);
     ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Incorrect 4 ms 376 KB Integer -1 violates the range [1, 6]
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 4 ms 376 KB Integer -1 violates the range [1, 49]
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 4 ms 380 KB Integer -1 violates the range [1, 199]
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 4 ms 376 KB Integer -1 violates the range [1, 999]
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 4 ms 376 KB Integer -1 violates the range [1, 9999]
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 4 ms 380 KB Integer -1 violates the range [1, 99999]
2 Halted 0 ms 0 KB -