oj.uz

Submission #966315

# Submission time Handle Problem Language Result Execution time Memory
966315 2024-04-19T17:13:56 Z MisterReaper Split the sequence (APIO14_sequence) C++17
11 / 100
 75 ms 18780 KB 
#include <bits/stdc++.h>

using i64 = long long;
#define int i64

struct line {
    i64 a, b;
    line(i64 _a, i64 _b) : a(_a), b(_b) {}
    long double intersect(line rhs) {
        return (long double)(b - rhs.b) / (rhs.a - a);
    }
    i64 eval(i64 x) {
        return a * x + b;
    }
};
namespace std {
    std::string to_string(line l) {
        return '(' + to_string(l.a) + "x + " + to_string(l.b) + ')';
    }
};

#ifdef DEBUG 
    #include "debug.h"
#else
    #define debug(...) void(23)
#endif

constexpr int maxN = 1E5 + 5;
constexpr int maxK = 200 + 5;

int A[maxN], go[maxK][maxN];
i64 dp[maxK][maxN], p[maxN];

signed main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);

    int N, K;
    std::cin >> N >> K;

    for(int i = 1; i <= N; i++) {
        std::cin >> A[i];
    }

    for(int i = 1; i <= N; i++) {
        p[i] = p[i - 1] + A[i];
    }

    std::deque<std::pair<line, int>> dq;
    auto size = [&]() -> int { return (int) dq.size(); };
    auto Add = [&](line l, int idx) -> void {
        while(size() >= 2 && l.intersect(dq[0].first) >= dq[0].first.intersect(dq[1].first)) {
            dq.pop_front();
        }
        dq.emplace_front(l, idx);
    };
    auto Query = [&](i64 x) -> std::pair<i64, int> {
        int l = 0, r = size() - 1;
        while(l < r) {
            int m = (l + r) >> 1;
            if(dq[m].first.eval(x) < dq[m + 1].first.eval(x)) {
                l = m + 1;
            } else {
                r = m;
            }
        }
        return {dq[l].first.eval(x), dq[l].second};
    };

    for(int k = 1; k <= K; k++) {
        dq.clear();
        Add({0, 0}, 0);
        for(int i = 1; i <= N; i++) {
            // for(int j = 0; j < i; j++) {
            //     if(dp[k - 1][j] + 1LL * (p[i] - p[j]) * (p[N] - p[i]) >= dp[k][i]) {
            //         go[k][i] = j;
            //         dp[k][i] = dp[k - 1][j] + 1LL * (p[i] - p[j]) * (p[N] - p[i]);
            //     }
            // }

            auto [v, j] = Query(p[N] - p[i]);
            dp[k][i] = v + 1LL * p[i] * (p[N] - p[i]);
            go[k][i] = j;
            Add(line{-p[i], dp[k - 1][i]}, i);
            debug(dq);

            debug(k, i, dp[k][i], go[k][i]);
        }
    }

    int ans = 1;
    for(int i = 1; i < N; i++) {
        if(dp[K][i] >= dp[K][ans]) {
            ans = i;
        }
    }

    std::cout << dp[K][ans] << '\n';

    std::vector<int> v;
    for(int k = K; k >= 1; k--) {
        v.emplace_back(ans);
        ans = go[k][ans];
    }

    for(int i = K - 1; i >= 0; i--) {
        std::cout << v[i] << " \n"[i == 0];
    }

    return 0;
}

Subtask #1
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 2 ms 10584 KB contestant found the optimal answer: 108 == 108
2 Correct 2 ms 10588 KB contestant found the optimal answer: 999 == 999
3 Correct 1 ms 6492 KB contestant found the optimal answer: 0 == 0
4 Correct 2 ms 14684 KB contestant found the optimal answer: 1542524 == 1542524
5 Correct 4 ms 18780 KB contestant found the optimal answer: 4500000000 == 4500000000
6 Correct 1 ms 6500 KB contestant found the optimal answer: 1 == 1
7 Correct 1 ms 6608 KB contestant found the optimal answer: 1 == 1
8 Correct 1 ms 6488 KB contestant found the optimal answer: 1 == 1
9 Correct 2 ms 10588 KB contestant found the optimal answer: 100400096 == 100400096
10 Correct 1 ms 10588 KB contestant found the optimal answer: 900320000 == 900320000
11 Correct 1 ms 10588 KB contestant found the optimal answer: 3698080248 == 3698080248
12 Correct 2 ms 12636 KB contestant found the optimal answer: 3200320000 == 3200320000
13 Correct 2 ms 12636 KB contestant found the optimal answer: 140072 == 140072
14 Correct 1 ms 10588 KB contestant found the optimal answer: 376041456 == 376041456
15 Correct 2 ms 12636 KB contestant found the optimal answer: 805 == 805
16 Correct 1 ms 10584 KB contestant found the optimal answer: 900189994 == 900189994
17 Correct 2 ms 12636 KB contestant found the optimal answer: 999919994 == 999919994

Subtask #2
0 / 11.0

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

Subtask #3
0 / 11.0

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

Subtask #4
0 / 17.0

# Verdict Execution time Memory Grader output
1 Incorrect 2 ms 10588 KB declared answer doesn't correspond to the split scheme: declared = 21306261, real = 21394949
2 Halted 0 ms 0 KB -

Subtask #5
0 / 21.0

# Verdict Execution time Memory Grader output
1 Incorrect 6 ms 10892 KB declared answer doesn't correspond to the split scheme: declared = 1606278379, real = 1624050854
2 Halted 0 ms 0 KB -

Subtask #6
0 / 29.0

# Verdict Execution time Memory Grader output
1 Incorrect 75 ms 14544 KB declared answer doesn't correspond to the split scheme: declared = 19737033504, real = 19779800813
2 Halted 0 ms 0 KB -