結果

問題 No.703 ゴミ拾い Easy
ユーザー kopricky
提出日時 2018-06-30 17:33:31
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,152 bytes
コンパイル時間 1,565 ms
コンパイル使用メモリ 174,812 KB
実行使用メモリ 10,512 KB
最終ジャッジ日時 2024-07-01 00:58:39
合計ジャッジ時間 4,718 ms
ジャッジサーバーID
(参考情報)
judge4 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample WA * 4
other AC * 5 WA * 41
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#pragma GCC optimize("O3")
#pragma GCC target("avx")
#define rep(i,n) for(int i = 0; i < (int)(n); i++)
#define ll long long
#define fi first
#define se second
#define push_back pb
#define show(x) cout << #x << " = " << x << "\n"

using namespace std;

typedef pair<int,int> P;

#define T long long

#define getchar getchar_unlocked

inline int in() {
    int n, c;
    while ((c = getchar()) < '0') if (c == EOF) return -1;
    n = c - '0';
    while ((c = getchar()) >= '0') n = n * 10 + c - '0';
    return n;
}

//f[j](x) = a[j]x + b[j]でaがjの増加に伴い単調減少する場合
class CHT
{
private:
    using ptt = pair<T, T>;
    bool check(ptt l1, ptt l2, ptt l3){
        return (l2.first-l1.first)*(l3.second-l2.second)>=(l2.second-l1.second)*(l3.first-l2.first);
    }
    T f(int i, T x){
        return lines[i].first * x + lines[i].second;
    }
public:
    vector<ptt> lines;
    CHT(){};
    void add(T a, T b){
        ptt line(a, b);
        while((int)lines.size() >= 2 && check(*(lines.end()-2), lines.back(), line)){
            lines.pop_back();
        }
        lines.emplace_back(line);
    }
    T get(T x){
        static int head = 0;
        while((int)lines.size() - head >= 2 && f(head, x) > f(head + 1, x)){
            head++;
        }
        return f(head, x);
        // int low = -1, high = lines.size() - 1;
		// while (high - low > 1) {
		// 	int mid = (high + low) / 2;
        //     if(f(mid, x) >= f(mid+1, x)){
        //         low = mid;
        //     }else{
        //         high = mid;
        //     }
		// }
		// return f(high, x);
    }
};

int main()
{
    cin.tie(0);
    ios::sync_with_stdio(false);
    int n;
    cin >> n;
    vector<ll> a(n),x(n),y(n);
    rep(i,n){
        a[i] = in();
    }
    rep(i,n){
        x[i] = in();
    }
    rep(i,n){
        y[i] = in();
    }
    CHT cht;
    cht.add(-2*x[0],x[0]*x[0]+y[0]*y[0]);
    rep(i,n){
        if(i < n-1){
            cht.add(-2*x[i+1],cht.get(a[i])+x[i+1]*x[i+1]+y[i+1]*y[i+1]+a[i]*a[i]);
        }else{
            cout << cht.get(a[i])+a[i]*a[i] << "\n";
            return 0;
        }
    }
}
0