#include using namespace std; #define overload2(a, b, c, ...) c #define overload3(a, b, c, d, ...) d #define overload4(a, b, c, d, e ...) e #define overload5(a, b, c, d, e, f ...) f #define overload6(a, b, c, d, e, f, g ...) g #define fast_io ios::sync_with_stdio(false); cin.tie(nullptr); #pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math") typedef long double ld; #define chmin(a,b) a = min(a,b); #define chmax(a,b) a = max(a,b); #define bit_count(x) __builtin_popcountll(x) #define leading_zero_count(x) __builtin_clz(x) #define trailing_zero_count(x) __builtin_ctz(x) #define gcd(a,b) __gcd(a,b) #define lcm(a,b) a / gcd(a,b) * b #define rep(...) overload3(__VA_ARGS__, rrep, rep1)(__VA_ARGS__) #define rep1(i,n) for(int i = 0 ; i < n ; i++) #define rrep(i,a,b) for(int i = a ; i < b ; i++) #define repi(it,S) for(auto it = S.begin() ; it != S.end() ; it++) #define pt(a) cout << a << endl; #define print(...) printall(__VA_ARGS__); #define debug(a) cout << #a << " " << a << endl; #define all(a) a.begin(), a.end() #define endl "\n"; #define v1(T,n,a) vector(n,a) #define v2(T,n,m,a) vector>(n,v1(T,m,a)) #define v3(T,n,m,k,a) vector>>(n,v2(T,m,k,a)) #define v4(T,n,m,k,l,a) vector>>>(n,v3(T,m,k,l,a)) templateistream &operator>>(istream&is,pair&p){is>>p.first>>p.second;return is;} templateostream &operator<<(ostream&os,const pair&p){os<istream &operator>>(istream&is,vector&v){for(T &in:v){is>>in;}return is;} templateostream &operator<<(ostream&os,const vector&v){for(auto it=v.begin();it!=v.end();){os<<*it<<((++it)!=v.end()?" ":"");}return os;} templateistream &operator>>(istream&is,vector>&v){for(T &in:v){is>>in;}return is;} templateostream &operator<<(ostream&os,const vector>&v){for(auto it=v.begin();it!=v.end();){os<<*it<<((++it)!=v.end()?"\n":"");}return os;} templateostream &operator<<(ostream&os,const set&v){for(auto it=v.begin();it!=v.end();){os<<*it<<((++it)!=v.end()?" ":"");}return os;} templateostream &operator<<(ostream&os,const multiset&v){for(auto it=v.begin();it!=v.end();){os<<*it<<((++it)!=v.end()?" ":"");}return os;} template void printall(Args... args){for(auto i:initializer_list>{args...}) cout<(const Int128t &a) const { return x > a.x; } bool operator<=(const Int128t &a) const { return x <= a.x; } bool operator>=(const Int128t &a) const { return x >= a.x; } bool operator==(const Int128t &a) const { return x == a.x; } bool operator!=(const Int128t &a) const { return x != a.x; } Int128t operator<<(const Int128t &a) const { return Int128t(x << a.x); } Int128t operator>>(const Int128t &a) const { return Int128t(x >> a.x); } Int128t& operator<<=(const Int128t &a) { x <<= a.x; return *this; } Int128t& operator>>=(const Int128t &a) { x >>= a.x; return *this; } friend istream &operator>>(istream &is, Int128t &x) { string s; is >> s; x.parse(s); return is; } friend ostream &operator<<(ostream &os, const Int128t &x) { if ( !ostream::sentry(os) ) return os; char buf[64]; char *d = end(buf); __uint128_t tmp = (x.x < 0 ? -x.x : x.x); do { d--; *d = char(tmp % 10 + '0'); tmp /= 10; } while ( tmp ); if ( x.x < 0 ) { d--; *d = '-'; } int len = end(buf) - d; if ( os.rdbuf()->sputn(d, len) != len ) { os.setstate(ios_base::badbit); } return os; } }; using ll = Int128t; int n; ll d_square(ll xg, ll yg, ll gx, ll gy) { return (xg - gx) * (xg - gx) + (yg - gy) * (yg - gy); } tuple center(vectorx, vectory){ ll xp = (y[0] - y[1]) * (x[2] * x[2] - x[0] * x[0] + y[2] * y[2] - y[0] * y[0]) - (y[0] - y[2]) * (x[1] * x[1] - x[0] * x[0] + y[1] * y[1] - y[0] * y[0]); ll yp = (x[0] - x[2]) * (x[1] * x[1] - x[0] * x[0] + y[1] * y[1] - y[0] * y[0]) - (x[0] - x[1]) * (x[2] * x[2] - x[0] * x[0] + y[2] * y[2] - y[0] * y[0]); ll cp = (x[0] - x[1]) * (y[0] - y[2]) - (x[0] - x[2]) * (y[0] - y[1]); cp *= 2; return {xp,yp,cp}; } void solve(){ int q; cin >> q; vector> P(3); cin >> P; rep(i,3) swap(P[i].first,P[i].second); rep(i,3) P[i].second *= -1; sort(all(P),greater>()); vector X(3), Y(3); rep(i,3){ auto[y,x] = P[i]; X[i] = -x; Y[i] = y; } auto[xg,yg,c] = center(X,Y); sort(all(X)); sort(all(Y)); if(X[0] + X[2] == X[1] * 2 && Y[0] + Y[2] == Y[1] * 2) { xg = X[0] + X[2]; yg = Y[0] + Y[2]; c = 2; } ll r_square = d_square(xg,yg,X[2]*c,Y[2]*c); rep(i,q){ ll x, y; cin >> x >> y; ll dist = d_square(xg,yg,x*c,y*c); if(dist <= r_square) { pt("Yes") } else { pt("No") } } } int main(){ fast_io int t = 1; // cin >> t; rep(i,t) solve(); }