結果

問題 No.913 木の燃やし方
ユーザー square1001square1001
提出日時 2020-02-13 22:31:58
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 315 ms / 3,000 ms
コード長 2,681 bytes
コンパイル時間 1,182 ms
コンパイル使用メモリ 88,280 KB
実行使用メモリ 18,116 KB
最終ジャッジ日時 2024-04-15 23:26:56
合計ジャッジ時間 13,304 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 5 ms
5,376 KB
testcase_04 AC 5 ms
5,376 KB
testcase_05 AC 6 ms
5,376 KB
testcase_06 AC 4 ms
5,376 KB
testcase_07 AC 3 ms
5,376 KB
testcase_08 AC 305 ms
17,492 KB
testcase_09 AC 289 ms
15,924 KB
testcase_10 AC 290 ms
16,808 KB
testcase_11 AC 295 ms
17,212 KB
testcase_12 AC 315 ms
17,996 KB
testcase_13 AC 303 ms
17,488 KB
testcase_14 AC 289 ms
15,796 KB
testcase_15 AC 290 ms
16,808 KB
testcase_16 AC 234 ms
17,616 KB
testcase_17 AC 224 ms
15,768 KB
testcase_18 AC 225 ms
16,936 KB
testcase_19 AC 250 ms
17,240 KB
testcase_20 AC 243 ms
18,004 KB
testcase_21 AC 301 ms
18,008 KB
testcase_22 AC 297 ms
17,992 KB
testcase_23 AC 284 ms
18,116 KB
testcase_24 AC 253 ms
16,320 KB
testcase_25 AC 224 ms
18,016 KB
testcase_26 AC 313 ms
16,484 KB
testcase_27 AC 268 ms
17,260 KB
testcase_28 AC 287 ms
17,232 KB
testcase_29 AC 287 ms
17,860 KB
testcase_30 AC 295 ms
17,880 KB
testcase_31 AC 267 ms
17,192 KB
testcase_32 AC 268 ms
17,260 KB
testcase_33 AC 270 ms
18,024 KB
testcase_34 AC 260 ms
17,100 KB
testcase_35 AC 251 ms
17,996 KB
testcase_36 AC 231 ms
18,008 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <set>
#include <deque>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
const long long inf = 1LL << 62;
vector<long long> subsolve(int N, vector<long long> A) {
	A.insert(A.begin(), 0);
	for (int i = 0; i < N; ++i) {
		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;
	}
	int mid = N / 2 + 1;
	vector<long long> seq(N, inf);
	deque<int> que;
	for (int i = 0; i < mid; ++i) {
		while (que.size() >= 2) {
			int pa = que[que.size() - 2];
			int pb = que[que.size() - 1];
			__int128 cra = __int128(intercept[pa] - intercept[i]) * (i - pb);
			__int128 crb = __int128(intercept[pb] - intercept[i]) * (i - pa);
			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;
		}
	}
	for (int i = mid; i <= N; ++i) {
		seq[i - 1] = que[0];
		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;
		}
	}
	for (int i = 0; i < N; ++i) {
		if (seq[i] != inf) {
			seq[i] = A[i + 1] - A[seq[i]] + 1LL * (i + 1 - seq[i]) * (i + 1 - seq[i]);
		}
	}
	for (int i = N - 2; i >= mid - 1; --i) {
		seq[i] = min(seq[i], seq[i + 1]);
	}
	return seq;
}
vector<long long> solve(int N, vector<long long> A) {
	if (N == 1) return vector<long long>({ A[0] + 1 });
	int mid = N / 2;
	vector<long long> res1 = subsolve(N, A);
	reverse(A.begin(), A.end());
	vector<long long> res2 = subsolve(N, A);
	reverse(res2.begin(), res2.end());
	for (int i = 0; i < N; ++i) {
		res1[i] = min(res1[i], res2[i]);
	}
	reverse(A.begin(), A.end());
	vector<long long> res3 = solve(mid, vector<long long>(A.begin(), A.begin() + mid));
	vector<long long> res4 = solve(N - mid, vector<long long>(A.begin() + mid, A.end()));
	for (int i = 0; i < mid; ++i) {
		res1[i] = min(res1[i], res3[i]);
	}
	for (int i = mid; i < N; ++i) {
		res1[i] = min(res1[i], res4[i - mid]);
	}
	return res1;
}
int main() {
	cin.tie(0);
	ios_base::sync_with_stdio(false);
	int N;
	cin >> N;
	vector<long long> A(N);
	for (int i = 0; i < N; ++i) cin >> A[i];
	vector<long long> ans = solve(N, A);
	for (int i = 0; i < N; ++i) {
		cout << ans[i] << '\n';
	}
	return 0;
}
0