結果

問題 No.2602 Real Collider
コンテスト
ユーザー Jeroen Op de Beek
提出日時 2024-01-12 22:21:59
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 45 ms / 2,000 ms
コード長 2,119 bytes
コンパイル時間 2,161 ms
コンパイル使用メモリ 195,220 KB
最終ジャッジ日時 2025-02-18 18:32:00
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 78
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
#define all(x) begin(x),end(x)
template<typename A, typename B> ostream& operator<<(ostream &os, const pair<A, B> &p) { return os << '(' << p.first << ", " << p.second << ')'; }
template<typename T_container, typename T = typename enable_if<!is_same<T_container, string>::value, typename T_container::value_type>::type> ostream& operator<<(ostream &os, const T_container &v) { string sep; for (const T &x : v) os << sep << x, sep = " "; return os; }
#define debug(a) cerr << "(" << #a << ": " << a << ")\n";
typedef long long ll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef pair<int,int> pi;
const int mxN = 1e5+1, oo = 1e9;
typedef complex<int> pt;
#define X real()
#define Y imag()
auto cross(pt u, pt v) {return (ll)u.X*v.Y-(ll)u.Y*v.X;}
auto sgn(ll a) {return a==0?0:(a>0?1:-1);}
auto ccw(pt p1, pt p2, pt p3) {auto u = p2-p1, v = p3-p2;return cross(u,v);}
auto in(pt p1, pt p2) {return (ll)p1.X*p2.X+(ll)p1.Y*p2.Y;}
auto norm2(pt p) {return (ll)p.X*p.X+(ll)p.Y*p.Y;}
bool comp(const pt& a, const pt& b) { return a.X<b.X or (a.X==b.X and a.Y < b.Y);}
void read(pt& p) {
    int a,b; cin >> a >> b;
    p = {a,b};
}
typedef __int128_t lll;
typedef pt P;
lll circV(P p, P a, P b, P c) { // is p in the circumcircle?
	lll p2 = norm2(p), A = norm2(a)-p2,
	    B = norm2(b)-p2, C = norm2(c)-p2;
	return ccw(p,a,b)*C + ccw(p,b,c)*A + ccw(p,c,a)*B;
}
bool incircle(pt p, pt a, pt b) {
    auto mid = (a+b)/2;
    return norm2(p-mid)*4 <= norm2(b-a);
}
int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
    int q; cin >> q;
    array<pt,3> cs;
    for(auto& i : cs) read(i);
    bool twop=0;
    for(int i=0;i<3;++i) {
        rotate(cs.begin(),cs.begin()+1,cs.end());
        if(incircle(cs[2]*2,cs[0]*2,cs[1]*2)) {
            twop=1;
            break;
        }
    }
    if(!twop) {
        if(ccw(cs[0],cs[1],cs[2])>0) swap(cs[0],cs[1]);
    }
    

    while(q--) {
        pt p; read(p);
        if(twop?incircle(p*2,cs[0]*2,cs[1]*2):circV(p,cs[0],cs[1],cs[2])<=0) {
            cout << "Yes\n";
        } else cout << "No\n";
    }
}
0