結果

問題 No.703 ゴミ拾い Easy
ユーザー fumiphysfumiphys
提出日時 2019-03-30 18:49:22
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 105 ms / 1,500 ms
コード長 2,729 bytes
コンパイル時間 1,028 ms
コンパイル使用メモリ 119,408 KB
実行使用メモリ 12,672 KB
最終ジャッジ日時 2024-06-10 18:38:45
合計ジャッジ時間 5,124 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 1 ms
6,944 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,944 KB
testcase_05 AC 2 ms
6,944 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 2 ms
6,944 KB
testcase_08 AC 1 ms
6,944 KB
testcase_09 AC 2 ms
6,944 KB
testcase_10 AC 2 ms
6,940 KB
testcase_11 AC 2 ms
6,944 KB
testcase_12 AC 2 ms
6,944 KB
testcase_13 AC 1 ms
6,944 KB
testcase_14 AC 2 ms
6,944 KB
testcase_15 AC 2 ms
6,940 KB
testcase_16 AC 2 ms
6,944 KB
testcase_17 AC 2 ms
6,944 KB
testcase_18 AC 2 ms
6,940 KB
testcase_19 AC 2 ms
6,940 KB
testcase_20 AC 2 ms
6,940 KB
testcase_21 AC 2 ms
6,940 KB
testcase_22 AC 2 ms
6,940 KB
testcase_23 AC 2 ms
6,940 KB
testcase_24 AC 102 ms
12,664 KB
testcase_25 AC 105 ms
12,672 KB
testcase_26 AC 102 ms
12,544 KB
testcase_27 AC 102 ms
12,672 KB
testcase_28 AC 104 ms
12,672 KB
testcase_29 AC 102 ms
12,588 KB
testcase_30 AC 103 ms
12,648 KB
testcase_31 AC 101 ms
12,672 KB
testcase_32 AC 102 ms
12,544 KB
testcase_33 AC 104 ms
12,632 KB
testcase_34 AC 75 ms
12,544 KB
testcase_35 AC 76 ms
12,500 KB
testcase_36 AC 75 ms
12,544 KB
testcase_37 AC 77 ms
12,672 KB
testcase_38 AC 78 ms
12,580 KB
testcase_39 AC 78 ms
12,544 KB
testcase_40 AC 77 ms
12,544 KB
testcase_41 AC 77 ms
12,544 KB
testcase_42 AC 76 ms
12,544 KB
testcase_43 AC 76 ms
12,544 KB
testcase_44 AC 2 ms
6,940 KB
testcase_45 AC 2 ms
6,944 KB
testcase_46 AC 86 ms
12,564 KB
testcase_47 AC 87 ms
12,572 KB
testcase_48 AC 2 ms
6,940 KB
testcase_49 AC 2 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// includes
#include <cstdio>
#include <cstdint>
#include <iostream>
#include <iomanip>
#include <string>
#include <queue>
#include <stack>
#include <vector>
#include <set>
#include <map>
#include <unordered_map>
#include <algorithm>
#include <utility>
#include <functional>
#include <cmath>
#include <climits>
#include <bitset>
#include <list>
#include <random>

// macros
#define ll long long int
#define pb emplace_back
#define mk make_pair
#define pq priority_queue
#define FOR(i, a, b) for(int i=(a);i<(b);++i)
#define rep(i, n) FOR(i, 0, n)
#define rrep(i, n) for(int i=((int)(n)-1);i>=0;i--)
#define all(x) (x).begin(),(x).end()
#define sz(x) ((int)(x).size())
#define UNIQUE(v) v.erase(unique(v.begin(), v.end()), v.end())

using namespace std;

//  types
typedef pair<int, int> P;
typedef pair<ll, int> Pl;
typedef pair<ll, ll> Pll;
typedef pair<double, double> Pd;
 
// constants
const int inf = 1e9;
const ll linf = 1LL << 50;
const double EPS = 1e-10;
const int mod = 1e9 + 7;

// solve
template <class T>bool chmax(T &a, const T &b){if(a < b){a = b; return 1;} return 0;}
template <class T>bool chmin(T &a, const T &b){if(a > b){a = b; return 1;} return 0;}

ll calc(ll a, ll x, ll y){
  return (x - a) * (x - a) + y * y;
}

template <typename T>
struct ConvexHullTrick_ {
  vector<pair<T, T> > lines;
  function<bool(T l, T r)> comp;
  ConvexHullTrick_(function<bool(T l, T r)> comp_=[](T l, T r){return l >= r;}): comp(comp_){}
  bool isnot_require(pair<T, T> l1, pair<T, T> l2, pair<T, T> l3){
    if(l1 < l3)swap(l1, l3);
    return (l3.second - l2.second) * (l2.first - l1.first) >= (l2.second - l1.second) * (l3.first - l2.first);
  }
  void add(pair<T, T> p){
    while(lines.size() >= 2 && isnot_require(*(lines.end()-2), *(lines.end()-1), p))lines.pop_back();
    lines.push_back(p);
  }
  T f(int i, T x){
    return lines[i].first * x + lines[i].second;
  }
  T val(T x){
    int ld = -1;
    int rd = lines.size() - 1;
    while(rd - ld > 1){
      int md = (rd + ld) / 2;
      if(comp(f(md, x), f(md + 1, x)))ld = md;
      else rd = md;
    }
    return f(rd, x);
  }
};

using ConvexHullTrickI = ConvexHullTrick_<int>;
using ConvexHullTrickL = ConvexHullTrick_<long long>;

ll dp[300001];

int main(int argc, char const* argv[])
{
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  int n;
  cin >> n;
  vector<ll> a(n), x(n), y(n);
  rep(i, n)cin >> a[i];
  rep(i, n)cin >> x[i];
  rep(i, n)cin >> y[i];
  ConvexHullTrickL cht;
  dp[0] = calc(a[0], x[0], y[0]);
  cht.add(mk(- 2 * x[0], x[0] * x[0] + y[0] * y[0]));
  for(int i = 1; i < n; i++){
    cht.add(mk(- 2 * x[i], x[i] * x[i] + y[i] * y[i] + dp[i-1]));
    dp[i] = a[i] * a[i] + cht.val(a[i]);
  }
  cout << dp[n - 1] << endl;
	return 0;
}
0