結果
問題 |
No.2602 Real Collider
|
ユーザー |
|
提出日時 | 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 |
ソースコード
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")