oj.uz

Submission #966121

# Submission time Handle Problem Language Result Execution time Memory
966121 2024-04-19T12:06:22 Z MisterReaper Split the sequence (APIO14_sequence) C++17
0 / 100
 105 ms 11644 KB 
#include <bits/stdc++.h>

using i64 = long long;

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 '(' + l.a + "x + " + 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_back(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.intersect(dq[m + 1].first) >= 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);
        }
    }

    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
0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 8792 KB contestant found the optimal answer: 108 == 108
2 Incorrect 1 ms 6488 KB declared answer doesn't correspond to the split scheme: declared = 839, real = 888
3 Halted 0 ms 0 KB -

Subtask #2
0 / 11.0

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

Subtask #3
0 / 11.0

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

Subtask #4
0 / 17.0

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

Subtask #5
0 / 21.0

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

Subtask #6
0 / 29.0

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