結果

問題 No.2602 Real Collider
ユーザー anonymously
提出日時 2024-01-13 06:49:29
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 1,376 ms / 2,000 ms
コード長 1,403 bytes
コンパイル時間 83 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 10,880 KB
最終ジャッジ日時 2024-09-28 01:12:38
合計ジャッジ時間 41,712 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 78
権限があれば一括ダウンロードができます

ソースコード

diff #

Q=int(input())
Xa,Ya,Xb,Yb,Xc,Yc=map(int,input().split())
#各頂点の「内積」
inA=(Xb-Xa)*(Xc-Xa)+(Yb-Ya)*(Yc-Ya)
inB=(Xc-Xb)*(Xa-Xb)+(Yc-Yb)*(Ya-Yb)
inC=(Xa-Xc)*(Xb-Xc)+(Ya-Yc)*(Yb-Yc)

#中心
X0=None
Y0=None
#det=1
det=2
if inA<0:
    #X0=(Xb+Xc)/2
    #Y0=(Yb+Yc)/2
    X0=(Xb+Xc)
    Y0=(Yb+Yc)
elif inB<0:
    X0=(Xc+Xa)
    Y0=(Yc+Ya)
elif inC<0:
    X0=(Xa+Xb)
    Y0=(Ya+Yb)
else:
    #(X0-Xa)**2+(Y0-Ya)**2==(X0-Xb)**2+(Y0-Yb)**2==(X0-Xc)**2+(Y0-Yc)**2
    #-2*Xa*X0+Xa**2-2*Ya*Y0+Ya**2==-2*Xb*X0+Xb**2-2*Yb*Y0+Yb**2==-2*Xc*X0+Xc**2-2*Yc*Y0+Yc**2
    #2*(Xb-Xa)*X0+2*(Yb-Ya)*Y0==Xb**2-Xa**2+Yb**2-Ya**2
    #2*(Xc-Xa)*X0+2*(Yc-Ya)*Y0==Xc**2-Xa**2+Yc**2-Ya**2
    det=4*((Xb-Xa)*(Yc-Ya)-(Yb-Ya)*(Xc-Xa))
    #X0=(2*(Yc-Ya)(Xb**2-Xa**2+Yb**2-Ya**2)-2*(Yb-Ya)(Xc**2-Xa**2+Yc**2-Ya**2))/det
    X0=(2*(Yc-Ya)*(Xb**2-Xa**2+Yb**2-Ya**2)-2*(Yb-Ya)*(Xc**2-Xa**2+Yc**2-Ya**2))
    Y0=(-2*(Xc-Xa)*(Xb**2-Xa**2+Yb**2-Ya**2)+2*(Xb-Xa)*(Xc**2-Xa**2+Yc**2-Ya**2))
#print((X0,Y0),det)
#print((det*Xa-X0)**2+(det*Ya-Y0)**2,(det*Xb-X0)**2+(det*Yb-Y0)**2,(det*Xc-X0)**2+(det*Yc-Y0)**2)
R2=(det*Xa-X0)**2+(det*Ya-Y0)**2 if inA>=0 else (det*Xb-X0)**2+(det*Yb-Y0)**2

for i in range(Q):
    x,y=map(int,input().split())
    if (det*x-X0)**2+(det*y-Y0)**2>R2:
        #print((det*x-X0)**2+(det*y-Y0)**2,"No")
        print("No")
    else:
        #print((det*x-X0)**2+(det*y-Y0)**2,"Yes")
        print("Yes")
0