結果

問題 No.2602 Real Collider
ユーザー anonymouslyanonymously
提出日時 2024-01-13 06:49:29
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
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
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 27 ms
10,752 KB
testcase_01 AC 31 ms
10,752 KB
testcase_02 AC 30 ms
10,624 KB
testcase_03 AC 31 ms
10,752 KB
testcase_04 AC 27 ms
10,752 KB
testcase_05 AC 30 ms
10,752 KB
testcase_06 AC 29 ms
10,752 KB
testcase_07 AC 27 ms
10,752 KB
testcase_08 AC 27 ms
10,752 KB
testcase_09 AC 29 ms
10,752 KB
testcase_10 AC 1,376 ms
10,880 KB
testcase_11 AC 512 ms
10,880 KB
testcase_12 AC 641 ms
10,752 KB
testcase_13 AC 332 ms
10,752 KB
testcase_14 AC 689 ms
10,752 KB
testcase_15 AC 390 ms
10,752 KB
testcase_16 AC 577 ms
10,752 KB
testcase_17 AC 656 ms
10,880 KB
testcase_18 AC 470 ms
10,624 KB
testcase_19 AC 558 ms
10,752 KB
testcase_20 AC 773 ms
10,880 KB
testcase_21 AC 470 ms
10,624 KB
testcase_22 AC 518 ms
10,752 KB
testcase_23 AC 363 ms
10,880 KB
testcase_24 AC 483 ms
10,752 KB
testcase_25 AC 535 ms
10,752 KB
testcase_26 AC 381 ms
10,880 KB
testcase_27 AC 596 ms
10,880 KB
testcase_28 AC 667 ms
10,752 KB
testcase_29 AC 530 ms
10,880 KB
testcase_30 AC 579 ms
10,752 KB
testcase_31 AC 575 ms
10,752 KB
testcase_32 AC 505 ms
10,752 KB
testcase_33 AC 587 ms
10,752 KB
testcase_34 AC 599 ms
10,752 KB
testcase_35 AC 374 ms
10,880 KB
testcase_36 AC 367 ms
10,624 KB
testcase_37 AC 634 ms
10,880 KB
testcase_38 AC 654 ms
10,752 KB
testcase_39 AC 638 ms
10,752 KB
testcase_40 AC 321 ms
10,880 KB
testcase_41 AC 747 ms
10,752 KB
testcase_42 AC 575 ms
10,880 KB
testcase_43 AC 563 ms
10,880 KB
testcase_44 AC 723 ms
10,880 KB
testcase_45 AC 469 ms
10,880 KB
testcase_46 AC 451 ms
10,752 KB
testcase_47 AC 636 ms
10,624 KB
testcase_48 AC 511 ms
10,752 KB
testcase_49 AC 411 ms
10,752 KB
testcase_50 AC 338 ms
10,752 KB
testcase_51 AC 355 ms
10,880 KB
testcase_52 AC 259 ms
10,752 KB
testcase_53 AC 613 ms
10,880 KB
testcase_54 AC 477 ms
10,624 KB
testcase_55 AC 539 ms
10,880 KB
testcase_56 AC 518 ms
10,752 KB
testcase_57 AC 470 ms
10,880 KB
testcase_58 AC 197 ms
10,752 KB
testcase_59 AC 548 ms
10,624 KB
testcase_60 AC 528 ms
10,752 KB
testcase_61 AC 401 ms
10,752 KB
testcase_62 AC 617 ms
10,752 KB
testcase_63 AC 705 ms
10,880 KB
testcase_64 AC 771 ms
10,752 KB
testcase_65 AC 408 ms
10,624 KB
testcase_66 AC 641 ms
10,880 KB
testcase_67 AC 317 ms
10,752 KB
testcase_68 AC 351 ms
10,752 KB
testcase_69 AC 246 ms
10,752 KB
testcase_70 AC 317 ms
10,880 KB
testcase_71 AC 383 ms
10,880 KB
testcase_72 AC 580 ms
10,880 KB
testcase_73 AC 479 ms
10,752 KB
testcase_74 AC 600 ms
10,752 KB
testcase_75 AC 631 ms
10,880 KB
testcase_76 AC 530 ms
10,880 KB
testcase_77 AC 572 ms
10,752 KB
testcase_78 AC 693 ms
10,880 KB
testcase_79 AC 604 ms
10,880 KB
testcase_80 AC 699 ms
10,880 KB
権限があれば一括ダウンロードができます

ソースコード

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