結果

問題 No.2376 障害物競プロ
ユーザー srjywrdnprktsrjywrdnprkt
提出日時 2023-07-25 18:14:12
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 901 ms / 4,000 ms
コード長 3,533 bytes
コンパイル時間 2,438 ms
コンパイル使用メモリ 202,620 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-04-10 09:58:28
合計ジャッジ時間 63,910 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 458 ms
6,944 KB
testcase_05 AC 624 ms
6,940 KB
testcase_06 AC 273 ms
6,940 KB
testcase_07 AC 901 ms
6,940 KB
testcase_08 AC 865 ms
6,940 KB
testcase_09 AC 893 ms
6,940 KB
testcase_10 AC 863 ms
6,940 KB
testcase_11 AC 813 ms
6,944 KB
testcase_12 AC 769 ms
6,940 KB
testcase_13 AC 837 ms
6,940 KB
testcase_14 AC 848 ms
6,940 KB
testcase_15 AC 798 ms
6,944 KB
testcase_16 AC 851 ms
6,944 KB
testcase_17 AC 778 ms
6,940 KB
testcase_18 AC 784 ms
6,944 KB
testcase_19 AC 767 ms
6,948 KB
testcase_20 AC 796 ms
6,944 KB
testcase_21 AC 799 ms
6,944 KB
testcase_22 AC 677 ms
6,944 KB
testcase_23 AC 474 ms
6,944 KB
testcase_24 AC 477 ms
6,940 KB
testcase_25 AC 229 ms
6,940 KB
testcase_26 AC 519 ms
6,940 KB
testcase_27 AC 460 ms
6,944 KB
testcase_28 AC 283 ms
6,944 KB
testcase_29 AC 226 ms
6,940 KB
testcase_30 AC 305 ms
6,940 KB
testcase_31 AC 284 ms
6,940 KB
testcase_32 AC 36 ms
6,944 KB
testcase_33 AC 117 ms
6,944 KB
testcase_34 AC 148 ms
6,940 KB
testcase_35 AC 93 ms
6,944 KB
testcase_36 AC 524 ms
6,940 KB
testcase_37 AC 609 ms
6,940 KB
testcase_38 AC 223 ms
6,940 KB
testcase_39 AC 675 ms
6,944 KB
testcase_40 AC 255 ms
6,940 KB
testcase_41 AC 202 ms
6,944 KB
testcase_42 AC 884 ms
6,944 KB
testcase_43 AC 891 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)) / 2;
}

//Distance between point a and b
template <typename T>
T distance(vec<T> &a, vec<T> &b){
    return sqrt(dot(a-b, a-b));
}

//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 s1, t1, s2, t2;
    s1 = outer(b-a, c-a);
    t1 = outer(b-a, d-a);
    s2 = outer(d-c, a-c);
    t2 = outer(d-c, b-c);
    if (s1 * t1 < 0 && s2 * t2 < 0) return 1;

    return 0;
}

int main(){

    using ld = long double;
    ll N, Q, a, b, c, d;
    cin >> N >> Q;
    ld x, y;
    vector<vec<ld>> q(N*2);
    vector<vector<ld>> dist(N*2, vector<ld>(N*2, 1e18));
    for (int i=0; i<N*2; i++){
        cin >> x >> y;
        q[i] = vec(x, y);
    }

    for (int i=0; i<N*2; i++){
        dist[i][i] = 0;
        for (int j=i+1; j<N*2; j++){
            bool f=1;
            for (int k=0; k<N; k++){
                if (intersect(q[i], q[j], q[2*k], q[2*k+1])) f=0;
            }
            if (f){
                dist[i][j] = dist[j][i] = distance(q[i], q[j]);
            }
        }
    }

    /*
    for (int i=0; i<N*2; i++){
        for (int j=0; j<N*2; j++) cout << dist[i][j] << " ";
        cout << endl;
    }
    */

    for (int k=0; k<N*2; k++){
        for (int i=0; i<N*2; i++){
            for (int j=0; j<N*2; j++){
                dist[i][j] = min(dist[i][j], dist[i][k]+dist[k][j]);
            }
        }
    }

    while(Q){
        Q--;
        cin >> a >> b >> c >> d; a--; b--; c--; d--;
        cout << setprecision(18) << dist[a*2+b][c*2+d] << endl;
    }

    return 0;
}
0