結果

問題 No.409 ダイエット
ユーザー 0w1
提出日時 2019-11-26 12:27:36
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 96 ms / 2,000 ms
コード長 1,725 bytes
コンパイル時間 805 ms
コンパイル使用メモリ 84,748 KB
最終ジャッジ日時 2025-01-08 05:38:57
ジャッジサーバーID
(参考情報)
judge3 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 92
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <iomanip>
#include <algorithm>
#include <deque>
using namespace std;

#define Inf 1e15
#define maxn 300010

long long d[maxn];
long long dp[maxn][2]; // dp[i][0]: sleep at i day, dp[i][1]: read at i day

struct line { // y = mx + b
	long long m; // slope
	long long b; // intercept
	line(long long m, long long b): m(m), b(b) {}
};

long long foo(line l, long long i, long long a, long long b) { // return dp[i][1] from line l(j)
	return l.m * i + l.b + b * i * (i + 1) / 2 - a * i;
}

double vertex(line l1, line l2) { // the intersection of l1 and l2
	return 1.0 * (l2.b - l1.b) / (l1.m - l2.m);
}

// reamin: w
// days: n
// study 1 day: -a
// study for x days: + x * b
// sleep in i day: + d[i]

int main() {
	long long n, a, b, w;
	cin >> n >> a >> b >> w;
	for (int i = 1; i <= n; i++)
		cin >> d[i];
	dp[0][0] = dp[0][1] = w;
	deque<line> dq;
	dq.push_front(line(0, w)); // sine dp[j][0] = w
	for (long long i = 1; i <= n; i++) {
		dp[i][0] = min(dp[i - 1][0], dp[i - 1][1]) + d[i]; // sleep
		// consider the value of dp[i][1] by subtituting i into dq[0]
		while (dq.size() > 1 && (foo(dq[1], i, a, b) < foo(dq[0], i, a, b)))
			dq.pop_front(); // we no longer need dq[0], pop out
		dp[i][1] = foo(dq[0], i, a, b);
		// cout << dp[i][0] << " " << dp[i][1] << " " << dq.size() << '\n';
		// add line "i" into dq
		line tmp = line(0, 0);
		tmp.m = -1 * i * b; // since this line will be treat as j in the future, the slope's j = i now
		tmp.b = dp[i][0] + b * i * (i - 1) / 2 + i * a;
		while (dq.size() >= 2 && (vertex(dq[dq.size() - 1], dq[dq.size() - 2]) > vertex(tmp, dq[dq.size() - 2])))
			dq.pop_back();
		dq.push_back(tmp);
	}
	cout << min(dp[n][1], dp[n][0]) << '\n';
	return 0;
}
0