結果

問題 No.705 ゴミ拾い Hard
ユーザー Caiiiiiiii
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #

#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;

}
0