結果
| 問題 |
No.705 ゴミ拾い Hard
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-03-26 01:26:30 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 100 ms / 1,500 ms |
| コード長 | 1,094 bytes |
| コンパイル時間 | 1,544 ms |
| コンパイル使用メモリ | 167,112 KB |
| 実行使用メモリ | 9,984 KB |
| 最終ジャッジ日時 | 2024-09-30 14:12:36 |
| 合計ジャッジ時間 | 5,884 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 40 |
ソースコード
#include <bits/stdc++.h>
typedef long long LL;
const int N = 3e5 + 7;
struct dat {
int l, r, p;
};
dat q[N], *fr, *ba;
int n, a[N], x[N], y[N];
LL f[N];
LL tri(LL x) {
return x * x * x;
}
LL F(int i, int j) {
return f[j] + tri(std::abs(a[i] - x[j + 1])) + tri(y[j + 1]);
}
int main() {
scanf("%d", &n);
for(int i = 1; i <= n; ++i)
scanf("%d", &a[i]);
for(int i = 1; i <= n; ++i)
scanf("%d", &x[i]);
for(int i = 1; i <= n; ++i)
scanf("%d", &y[i]);
*(fr = ba = q + 1) = {1, n, 0};
for(int i = 1; i <= n; ++i) {
f[i] = F(i, fr->p);
if(++fr->l > fr->r)
++fr;
while(fr <= ba && F(ba->l, i) <= F(ba->l, ba->p))
--ba;
if(fr > ba)
*++ba = {i + 1, n, i};
else if(F(ba->r, i) <= F(ba->r, ba->p)) {
int l = ba->l, r = ba->r;
while(l + 1 < r) {
int mid = l + r >> 1;
if(F(mid, i) <= F(mid, ba->p))
r = mid;
else
l = mid;
}
ba->r = r - 1;
*++ba = {r, n, i};
}
else if(ba->r != n) {
ba[1] = {ba->r + 1, n, i};
++ba;
}
}
printf("%lld\n", f[n]);
return 0;
}