結果
| 問題 |
No.409 ダイエット
|
| ユーザー |
|
| 提出日時 | 2016-08-06 13:22:48 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,288 bytes |
| コンパイル時間 | 2,048 ms |
| コンパイル使用メモリ | 174,336 KB |
| 実行使用メモリ | 12,416 KB |
| 最終ジャッジ日時 | 2024-11-07 03:10:35 |
| 合計ジャッジ時間 | 6,015 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 74 WA * 18 |
ソースコード
#include "bits/stdc++.h"
using namespace std;
#define FOR(i,j,k) for(int (i)=(j);(i)<(int)(k);++(i))
#define rep(i,j) FOR(i,0,j)
#define sz(x) (int)(x).size()
typedef long long ll;
int N;
ll A, B, W, INF = LLONG_MAX;
int main(){
cin >> N >> A >> B >> W;
vector<ll> D(N + 1), dp(N + 1, INF), a(N + 1), b(N + 1);
rep(i, N)scanf("%lld", &D[i+1]);
dp[0] = b[0] = W;
deque<int> l;
l.push_back(0);
for(int i = 1; i <= N; ++i){
ll a1, a2, a3, b1, b2, b3, n;
while(sz(l) >= 2){
a1 = a[l[0]], a2 = a[l[1]];
b1 = b[l[0]], b2 = b[l[1]];
if(a1*i + b1 >= a2*i + b2)l.pop_front();
else break;
}
dp[i] = D[i] - A*(i - 1) + B*i*(i - 1) / 2;
dp[i] += a[l[0]] * i + b[l[0]];
a3 = a[i] = -B*i;
b3 = b[i] = B*i*(i + 1) / 2 + dp[i] + A*i;
while(sz(l) >= 2){
n = sz(l);
a1 = a[l[n-2]], a2 = a[l[n-1]];
b1 = b[l[n-2]], b2 = b[l[n-1]];
if((a2 - a1)*(b3 - b2) >= (b2 - b1)*(a3 - a2))l.pop_back();
break;
}
l.push_back(i);
}
ll ans = INF;
for(int i = 0; i <= N; ++i){
int n = N - i;
ans = min(ans, dp[i] - A*n + B*n*(n + 1) / 2);
}
cout << ans << endl;
}