結果

問題 No.2602 Real Collider
ユーザー cho435
提出日時 2024-01-12 23:07:48
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,662 bytes
コンパイル時間 4,183 ms
コンパイル使用メモリ 253,844 KB
最終ジャッジ日時 2025-02-18 18:58:42
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 WA * 2
other AC * 65 WA * 13
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using ll = long long;
#define rep(i,n) for(int i=0;i<(int)(n);i++)
using mint = atcoder::modint998244353;
using ld = long double;
int main(){
    int q;
    cin>>q;
    vector<pair<int,int>> vp(3);
    vector<pair<int,int>> vq(q);
    for(auto&[x,y]:vp) cin>>x>>y;
    for(auto&[x,y]:vq) cin>>x>>y;
    ld cx,cy,r2=-1;
    rep(i,3){
        int x1,x2,y1,y2;
        x1=vp.at(1).first-vp.at(0).first;
        x2=vp.at(2).first-vp.at(0).first;
        y1=vp.at(1).second-vp.at(0).second;
        y2=vp.at(2).second-vp.at(0).second;
        if(x1*x2+y1*y2<=0){
            cout<<x1*x2+y1*y2<<endl;
            cx=(vp.at(1).first+vp.at(2).first);
            cy=(vp.at(1).second+vp.at(2).second);
            cx/=2;
            cy/=2;
            r2=(vp.at(1).first-vp.at(2).first)*(vp.at(1).first-vp.at(2).first)
                    +(vp.at(1).second-vp.at(2).second)*(vp.at(1).second-vp.at(2).second);
            r2/=4;
            break;
        }
        swap(vp.at(0),vp.at(1));
        swap(vp.at(1),vp.at(2));
    }
    ld eps=1e-10;
    if(r2<0){
        ld k=0;
        k+=(vp.at(0).first-vp.at(2).first)*(vp.at(1).first-vp.at(2).first);
        k+=(vp.at(0).second-vp.at(2).second)*(vp.at(1).second-vp.at(2).second);
        ld dv=0;
        dv+=(vp.at(0).first-vp.at(2).first)*(vp.at(1).second-vp.at(2).second);
        dv-=(vp.at(0).second-vp.at(2).second)*(vp.at(1).first-vp.at(2).first);
        assert(abs(dv)>eps);
        k/=dv;
        cx=(vp.at(0).first+vp.at(1).first)+k*(vp.at(0).second-vp.at(1).second);
        cx/=2;
        cy=(vp.at(0).second+vp.at(1).second)-k*(vp.at(0).first-vp.at(1).first);
        cy/=2;
        r2=(vp.at(0).first-cx)*(vp.at(0).first-cx)+(vp.at(0).second-cy)*(vp.at(0).second-cy);
    }
    for(auto[x,y]:vq){
        ld d2=(x-cx)*(x-cx)+(y-cy)*(y-cy);
        if(d2<r2+eps) cout<<"Yes\n";
        else cout<<"No\n";
    }
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0