結果
問題 |
No.2602 Real Collider
|
ユーザー |
![]() |
提出日時 | 2024-01-13 16:57:05 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,406 ms / 2,000 ms |
コード長 | 1,212 bytes |
コンパイル時間 | 236 ms |
コンパイル使用メモリ | 82,600 KB |
実行使用メモリ | 91,360 KB |
最終ジャッジ日時 | 2024-09-28 01:35:03 |
合計ジャッジ時間 | 52,902 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 78 |
ソースコード
Q=int(input()) ax,ay,bx,by,cx,cy=map(int, input().split()) A=[] from fractions import Fraction x=Fraction(bx+cx,2) y=Fraction(by+cy,2) r2=(bx-cx)**2+(by-cy)**2 if (ax-x)**2+(ay-y)**2<=Fraction(r2,4): A.append((r2,x,y)) ax,ay,bx,by,cx,cy=bx,by,cx,cy,ax,ay x=Fraction(bx+cx,2) y=Fraction(by+cy,2) r2=(bx-cx)**2+(by-cy)**2 if (ax-x)**2+(ay-y)**2<=Fraction(r2,4): A.append((r2,x,y)) ax,ay,bx,by,cx,cy=bx,by,cx,cy,ax,ay x=Fraction(bx+cx,2) y=Fraction(by+cy,2) r2=(bx-cx)**2+(by-cy)**2 if (ax-x)**2+(ay-y)**2<=Fraction(r2,4): A.append((r2,x,y)) ax,ay,bx,by,cx,cy=bx,by,cx,cy,ax,ay def c(x1, y1, x2, y2, x3, y3): """ 3点を通る円の中心と半径を取得 """ d = 2 * ((y1 - y3) * (x1 - x2) - (y1 - y2) * (x1 - x3)) if d==0: return 0,0,-1 x=Fraction(((y1-y3)*(y1**2-y2**2+x1**2-x2**2)-(y1-y2)*(y1**2-y3**2+x1**2-x3**2)),d) y=Fraction(((x1-x3)*(x1**2-x2**2+y1**2-y2**2)-(x1-x2)*(x1**2-x3**2+y1**2-y3**2)),-d) r=((x - x1) ** 2 + (y - y1) ** 2) return x,y,r x,y,r=c(ax,ay,bx,by,cx,cy) if r!=-1: A.append((4*r,x,y)) A=sorted(A) r,x,y=A[0] for i in range(Q): xx,yy=map(int, input().split()) if ((x-xx)**2+(y-yy)**2)*4>r: print('No') else: print('Yes')