結果

問題 No.2331 Maximum Quadrilateral
ユーザー srjywrdnprktsrjywrdnprkt
提出日時 2023-07-25 00:12:00
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 216 ms / 2,000 ms
コード長 3,172 bytes
コンパイル時間 1,812 ms
コンパイル使用メモリ 197,400 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-04-09 19:18:40
合計ジャッジ時間 6,552 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 185 ms
6,812 KB
testcase_01 AC 118 ms
6,940 KB
testcase_02 AC 129 ms
6,944 KB
testcase_03 AC 131 ms
6,940 KB
testcase_04 AC 168 ms
6,940 KB
testcase_05 AC 178 ms
6,944 KB
testcase_06 AC 167 ms
6,944 KB
testcase_07 AC 151 ms
6,944 KB
testcase_08 AC 105 ms
6,944 KB
testcase_09 AC 96 ms
6,944 KB
testcase_10 AC 87 ms
6,944 KB
testcase_11 AC 95 ms
6,940 KB
testcase_12 AC 87 ms
6,940 KB
testcase_13 AC 92 ms
6,940 KB
testcase_14 AC 182 ms
6,944 KB
testcase_15 AC 16 ms
6,944 KB
testcase_16 AC 62 ms
6,944 KB
testcase_17 AC 5 ms
6,944 KB
testcase_18 AC 19 ms
6,944 KB
testcase_19 AC 6 ms
6,940 KB
testcase_20 AC 203 ms
6,940 KB
testcase_21 AC 201 ms
6,940 KB
testcase_22 AC 205 ms
6,940 KB
testcase_23 AC 216 ms
6,940 KB
testcase_24 AC 205 ms
6,944 KB
testcase_25 AC 2 ms
6,944 KB
testcase_26 AC 2 ms
6,944 KB
testcase_27 AC 2 ms
6,940 KB
testcase_28 AC 2 ms
6,940 KB
testcase_29 AC 2 ms
6,940 KB
testcase_30 AC 2 ms
6,944 KB
testcase_31 AC 2 ms
6,940 KB
testcase_32 AC 1 ms
6,944 KB
testcase_33 AC 2 ms
6,944 KB
testcase_34 AC 2 ms
6,944 KB
testcase_35 AC 2 ms
6,940 KB
testcase_36 AC 2 ms
6,944 KB
testcase_37 AC 2 ms
6,940 KB
testcase_38 AC 2 ms
6,944 KB
testcase_39 AC 2 ms
6,944 KB
testcase_40 AC 1 ms
6,944 KB
testcase_41 AC 1 ms
6,944 KB
testcase_42 AC 1 ms
6,940 KB
testcase_43 AC 2 ms
6,944 KB
testcase_44 AC 2 ms
6,940 KB
testcase_45 AC 1 ms
6,944 KB
testcase_46 AC 1 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
using ll = long long;

// vector OA = (x, y)
template <typename T>
struct vec{
    T x, y;

    vec (T xx=0, T yy=0) : x(xx), y(yy) {};
    vec operator-() const {
      return vec(-x, -y);
    }
    vec& operator+=(const vec &w) {
        x += w.x; y += w.y;
        return *this;
    }
    vec& operator-=(const vec &w) {
        x -= w.x; y -= w.y;
        return *this;
    }
    vec operator+(const vec &w) const {
        vec res(*this);
        return res += w;
    }
    vec operator-(const vec &w) const {
        vec res(*this);
        return res-=w;
    }
};

//Inner product of vectors v and w
template <typename T>
T dot(vec<T> v, vec<T> w){
    return v.x * w.x + v.y * w.y;
}

//Outer product of vector v and w
template <typename T>
T outer(vec<T> v, vec<T> w){
    return v.x * w.y - w.x * v.y;
}

//size of triangle ABC
template <typename T>
T heron(vec<T> &a, vec<T> &b, vec<T> &c){
    return abs(outer(b-a, c-a));
}

//Convex Hull(Smallest Convex set containing all given vecs)
//Grahum Scan (O(NlogN))
template <typename T>
vector<vec<T>> convex_hull(vector<vec<T>> P){

    sort(P.begin(), P.end(), [](vec<T> &p1, vec<T> &p2) {
        if (p1.x != p2.x) return p1.x < p2.x;
        return p1.y < p2.y;
    });

    int N=P.size(), k=0, t;
    vector<vec<T>> res(N*2);

    for (int i=0; i<N; i++){
        while(k > 1 && outer(res[k-1]-res[k-2], P[i]-res[k-1]) <= 0) k--;
        res[k] = P[i];
        k++;
    }
    t = k;

    for (int i=N-2; i>=0; i--){
        while(k > t && outer(res[k-1]-res[k-2], P[i]-res[k-1]) <= 0) k--;
        res[k] = P[i];
        k++;
    }
    res.resize(k-1);

    return res;
}

//distance between line PQ and point R
template <typename T>
T dist_line_point(vec<T> &p, vec<T> &q, vec<T> &r){
    // ax+by+c=0
    T a = (q.y-p.y), b = (p.x-q.x), c = -p.x*q.y + p.y*q.x;
    cout << a << " " << b << " " << c << endl;
    return abs(a*r.x+b*r.y+c) / sqrt(a*a+b*b);
};

//judge if line AB intersects line CD
template <typename T>
bool intersect(vec<T> &a, vec<T> &b, vec<T> &c, vec<T> &d){
    T s, t;
    s = outer(b-a, c-a);
    t = outer(b-a, d-a);
    if (s == 0 && t == 0){
        if (max(a.x, b.x)<min(c.x, d.x) || max(a.y, b.y)<min(c.y, d.y) || min(a.x, b.x)>max(c.x, d.x) || min(a.y, b.y)>max(c.y, d.y)) return 0;
        else return 1;
    }
    if ((s > 0 && t > 0) || (s < 0 && t < 0)) return 0;
    s = outer(d-c, a-c);
    t = outer(d-c, b-c);
    if ((s > 0 && t > 0) || (s < 0 && t < 0)) return 0;

    return 1;
}

int main(){

    ll N, x, y, u, l, ans=0;
    cin >> N;
    vector<vec<ll>> p(N);
    for (int i=0; i<N; i++){
        cin >> x >> y;
        p[i] = vec(x, y);
    }

    for (int i=0; i<N; i++){
        for (int j=i+1; j<N; j++){
            l = -1e18, u = -1e18;
            for (int k=0; k<N; k++){
                if (i == k || j == k) continue;
                //upper side
                if (outer(p[j]-p[i], p[k]-p[i]) > 0) u = max(u, heron(p[i], p[j], p[k]));
                else l = max(l, heron(p[i], p[j], p[k]));
                ans = max(ans, l+u);
            }
        }
    }

    cout << ans << endl;

    return 0;
}
0