結果

問題 No.703 ゴミ拾い Easy
コンテスト
ユーザー vwxyz
提出日時 2024-12-04 16:23:56
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 347 ms / 1,500 ms
コード長 5,020 bytes
コンパイル時間 4,572 ms
コンパイル使用メモリ 280,068 KB
実行使用メモリ 10,368 KB
最終ジャッジ日時 2024-12-04 16:24:11
合計ジャッジ時間 13,798 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 46
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using uint = unsigned int;
using ll = long long;
using ull = unsigned long long;
template<class T> using V = vector<T>;
template<class T> using VV = V<V<T>>;
constexpr ll TEN(int n) { return (n == 0) ? 1 : 10 * TEN(n-1); }
#define FOR(i, a, b) for(int i=(int)(a);i<(int)(b);i++)
#define rep(i,N) for(int i=0;i<(int)(N);i++)
#define rep1(i,N) for(int i=1;i<=(int)(N);i++)
#define fs first
#define sc second
#define eb emplace_back
#define pb eb
#define all(x) x.begin(),x.end()
template<class T, class U> void chmin(T& t, const U& u) { if (t > u) t = u; }
template<class T, class U> void chmax(T& t, const U& u) { if (t < u) t = u; }
// bit op
int popcnt(uint x) { return __builtin_popcount(x); }
int popcnt(ull x) { return __builtin_popcountll(x); }
int bsr(uint x) { return 31 - __builtin_clz(x); }
int bsr(ull x) { return 63 - __builtin_clzll(x); }
int bsf(uint x) { return __builtin_ctz(x); }
int bsf(ull x) { return __builtin_ctzll(x); }
template<class T>

//que_incrはmax_yに入るものが単調増加ならtrue,単調減少ならfalse
struct ConvexHull {
  using L = array<T, 2>;
  bool que_incr;
  ConvexHull(bool _que_incr) : que_incr(_que_incr) {}
  deque<L> lines;
  // can remove mid?
  static bool is_need(L mid, L left, L right) {
    assert(left[0] <= mid[0] && mid[0] <= right[0]);
    return (right[0]-mid[0])*(left[1]-mid[1]) < (mid[0]-left[0])*(mid[1]-right[1]);
  }
  //work with 2^(60 + 64)
  /*static bool is_need(L mid, L left, L right) {
  assert(left[0] <= mid[0] && mid[0] <= right[0]);
  ll a = (right[0]-mid[0]), b = (left[1]-mid[1]), c = (mid[0]-left[0]), d
  = (mid[1]-right[1]);
  long double x = (long double)(a) * b - (long double)(c) * d;
  if (abs(x) > (1LL << 60)) return x < 0;
  int fl = b < 0, fr = d < 0;
  if (fl != fr) return fl == 1;
  ull z = ull(a) * ull(abs(b)) - ull(c) * ull(abs(d));
  if (fl == 0) return (1ULL << 63) < z;
  return z < (1ULL << 63);
  }*/
  void insert_front(L l) {
    if (lines.empty()) {
      lines.push_front(l);
      return;
    }
    assert(l[0] <= lines[0][0]);
    if (l[0] == lines[0][0]) {
      if (l[1] <= lines[0][1]) return;
      lines.pop_front();
    }
    while (lines.size() >= 2 && !is_need(lines.front(), l, lines[1])) {
      lines.pop_front();
    }
    lines.push_front(l);
  }
  void insert_back(L l) {
    if (lines.empty()) {
      lines.push_back(l);
      return;
    }
    assert(lines.back()[0] <= l[0]);
    if (lines.back()[0] == l[0]) {
      if (l[1] <= lines.back()[1]) return;
      lines.pop_back();
    }
    while (lines.size() >= 2 && !is_need(lines.back(),lines[lines.size()-2], l)) {
      lines.pop_back();
    }
    lines.push_back(l);
  }
  /**
  Insert line
  line's degree must be minimum or maximum
  */
  void insert_line(L line) {
    if (lines.empty()) {
      lines.push_back(line);
      return;
    }
    if (line[0] <= lines[0][0]) insert_front(line);
    else if (lines.back()[0] <= line[0]) insert_back(line);
    else assert(false); //line's degree must be minimum or maximum
  }
  /// get maximum y
  T b_x;
  T first = true;
  T max_y(T x) {
    assert(lines.size());
    auto value = [&](L l) { return l[0] * x + l[1]; };
    if (que_incr) {
      assert(first || b_x <= x);
      first = false; b_x = x;
      while (lines.size() >= 2 &&value(lines[0]) <= value(lines[1])) {
        lines.pop_front();
      }
      return value(lines.front());
    } else {
      assert(first || x <= b_x);
      first = false; b_x = x;
      while (lines.size() >= 2 && value(lines[lines.size()-2]) >= value(lines.back())) {
        lines.pop_back();
      }
      return value(lines.back());
    }
  }
};


struct Line {
  mutable ll k, m, p;
  bool operator<(const Line& o) const { return k < o.k; }
  bool operator<(ll x) const { return p < x; }
};

struct LineContainer : multiset<Line, less<>> {
  // (for doubles, use inf = 1/.0, div(a,b) = a/b)
  const ll inf = LLONG_MAX;
  ll div(ll a, ll b) { // floored division
    return a / b - ((a ^ b) < 0 && a % b); }
  bool isect(iterator x, iterator y) {
    if (y == end()) { x->p = inf; return false; }
    if (x->k == y->k) x->p = x->m > y->m ? inf : -inf;
    else x->p = div(y->m - x->m, x->k - y->k);
    return x->p >= y->p;
  }
  void add(ll k, ll m) {
    auto z = insert({k, m, 0}), y = z++, x = y;
    while (isect(y, z)) z = erase(z);
    if (x != begin() && isect(--x, y)) isect(x, y = erase(y));
    while ((y = x) != begin() && (--x)->p >= y->p)
    isect(x, erase(y));
  }
  ll query(ll x) {//gives max value
    assert(!empty());
    auto l = *lower_bound(x);
    return l.k * x + l.m;
  }
};

int main(){
  int N;cin>>N;
  vector<ll> A(N);
  vector<ll> X(N);
  vector<ll> Y(N);
  for(int i=0;i<N;i++){
    cin>>A[i];
  }
  for(int i=0;i<N;i++){
    cin>>X[i];
  }
  for(int i=0;i<N;i++){
    cin>>Y[i];
  }
  LineContainer dp;
  ll cost=0;
  for(int i=0;i<N;i++){
    dp.add(2*X[i],-cost-X[i]*X[i]-Y[i]*Y[i]);
    cost=-dp.query(A[i])+A[i]*A[i];
  }
  ll ans=cost;
  cout<<ans<<"\n";
}
0