結果

問題 No.913 木の燃やし方
ユーザー square1001square1001
提出日時 2019-10-18 23:13:28
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 3,424 bytes
コンパイル時間 1,539 ms
コンパイル使用メモリ 119,324 KB
実行使用メモリ 120,788 KB
最終ジャッジ日時 2023-09-08 02:16:19
合計ジャッジ時間 21,789 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 1 ms
4,376 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 WA -
testcase_04 WA -
testcase_05 AC 7 ms
4,380 KB
testcase_06 WA -
testcase_07 AC 3 ms
4,376 KB
testcase_08 AC 535 ms
58,232 KB
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 AC 499 ms
49,408 KB
testcase_13 AC 232 ms
36,128 KB
testcase_14 AC 223 ms
35,340 KB
testcase_15 AC 222 ms
34,400 KB
testcase_16 AC 2,268 ms
120,788 KB
testcase_17 AC 2,154 ms
118,736 KB
testcase_18 AC 2,267 ms
106,908 KB
testcase_19 AC 1,216 ms
111,216 KB
testcase_20 AC 2,194 ms
113,008 KB
testcase_21 WA -
testcase_22 WA -
testcase_23 AC 2,332 ms
112,888 KB
testcase_24 AC 1,765 ms
111,400 KB
testcase_25 AC 1,087 ms
112,052 KB
testcase_26 AC 224 ms
37,032 KB
testcase_27 AC 219 ms
38,252 KB
testcase_28 AC 220 ms
37,964 KB
testcase_29 AC 232 ms
37,804 KB
testcase_30 AC 597 ms
83,076 KB
testcase_31 AC 260 ms
45,792 KB
testcase_32 AC 221 ms
37,292 KB
testcase_33 AC 226 ms
38,024 KB
testcase_34 AC 914 ms
63,656 KB
testcase_35 AC 945 ms
67,916 KB
testcase_36 AC 2,078 ms
112,780 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <set>
#include <deque>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
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;
	}
	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[max(i - 5, 0)]; j <= seq[i]; ++j) {
			if (j != i) {
				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);
		}
	}
	return ans;
}
vector<long long> solve(int N, vector<long long> A) {
	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]);
	}
	return res1;
}
vector<long long> solve_easy(int N, vector<long long> A) {
	vector<long long> ans(N, 1LL << 62);
	for (int i = 0; i < N; ++i) {
		for (int j = i + 1; j <= N; ++j) {
			long long sum = 1LL * (j - i) * (j - i);
			for (int k = i; k < j; ++k) {
				sum += A[k];
			}
			for (int k = i; k < j; ++k) {
				ans[k] = min(ans[k], sum);
			}
		}
	}
	return ans;
}
#include <random>
#include <string>
using namespace std;
mt19937 mt(1910182301);
int rand_rng(int l, int r) {
	uniform_int_distribution<int> p(l, r - 1);
	return p(mt);
}
string stringify(vector<long long> A) {
	string res = "[";
	for (int i = 0; i < A.size(); ++i) {
		if (i) res += ", ";
		res += to_string(A[i]);
	}
	res += "]";
	return res;
}
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';
	}
	/*
	int N = 3;
	for (int i = 1; i <= 10000; ++i) {
		vector<long long> A(N);
		for (int j = 0; j < N; ++j) {
			A[j] = rand_rng(-6, 7);
		}
		vector<long long> res1 = solve(N, A);
		vector<long long> res2 = solve_easy(N, A);
		if (res1 != res2) {
			cout << "Case #" << i << ":" << endl;
			cout << "N = " << N << ", A = " << stringify(A) << endl;
			cout << "Returns: " << stringify(res1) << endl;
			cout << "Answer: " << stringify(res2) << endl;
			cout << endl;
		}
	}
	*/
	return 0;
}
0