結果

問題 No.913 木の燃やし方
ユーザー square1001
提出日時 2019-10-18 22:57:44
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,770 bytes
コンパイル時間 1,337 ms
コンパイル使用メモリ 97,888 KB
実行使用メモリ 42,428 KB
最終ジャッジ日時 2024-06-25 18:31:08
合計ジャッジ時間 11,215 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 9 WA * 25
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <set>
#include <deque>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
int main() {
    cin.tie(0);
    ios_base::sync_with_stdio(false);
    int N;
    cin >> N;
    vector<long long> A(N + 1);
    for (int i = 0; i < N; ++i) {
        cin >> A[i + 1];
        A[i + 1] += A[i];
    }
    // minimize (A[i] - A[j]) + (i - j)^2 = (A[i] + i * i) + (-A[j] + j * j) - 2 * j * i for each i
    vector<long long> intercept(N + 1), slope(N + 1);
    for (int i = 0; i <= N; ++i) {
        intercept[i] = -A[i] + 1LL * i * i;
        slope[i] = -2 * i;
    }
    vector<int> seq(N + 1);
    deque<int> que;
    for (int i = 0; i <= N; ++i) {
        if (i >= 1) {
            seq[i] = que[0];
        }
        while (que.size() >= 2) {
            int pa = que[que.size() - 2];
            int pb = que[que.size() - 1];
            double cra = -double(intercept[pa] - intercept[i]) / (slope[pa] - slope[i]);
            double crb = -double(intercept[pb] - intercept[i]) / (slope[pb] - slope[i]);
            if (cra > crb) que.pop_back();
            else break;
        }
        que.push_back(i);
        while (que.size() >= 2) {
            int pa = que[0];
            int pb = que[1];
            long long va = 1LL * slope[pa] * (i + 1) + intercept[pa];
            long long vb = 1LL * slope[pb] * (i + 1) + intercept[pb];
            if (va > vb) que.pop_front();
            else break;
        }
    }
    vector<vector<long long> > qs;
    for (int i = 1; i <= N; ++i) {
        for (int j = seq[i - 1]; j <= seq[i]; ++j) {
            qs.push_back({ A[i] - A[j] + 1LL * (i - j) * (i - j), j, i });
        }
    }
    sort(qs.begin(), qs.end());
    set<int> st;
    for (int i = 0; i < N; ++i) {
        st.insert(i);
    }
    vector<long long> ans(N);
    for (vector<long long> i : qs) {
        set<int>::iterator it = st.lower_bound(i[1]);
        while (it != st.end() && *it < i[2]) {
            ans[*it] = i[0];
            it = st.erase(it);
        }
    }
    for (int i = 0; i < N; ++i) {
        cout << ans[i] << '\n';
    }
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0