結果

問題 No.2602 Real Collider
ユーザー 👑 rin204rin204
提出日時 2024-01-12 22:28:52
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 895 ms / 2,000 ms
コード長 1,489 bytes
コンパイル時間 319 ms
コンパイル使用メモリ 82,044 KB
実行使用メモリ 77,008 KB
最終ジャッジ日時 2024-09-27 23:03:08
合計ジャッジ時間 30,238 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 36 ms
52,624 KB
testcase_01 AC 35 ms
53,580 KB
testcase_02 AC 34 ms
52,808 KB
testcase_03 AC 37 ms
54,744 KB
testcase_04 AC 36 ms
53,964 KB
testcase_05 AC 35 ms
52,460 KB
testcase_06 AC 35 ms
53,808 KB
testcase_07 AC 34 ms
53,380 KB
testcase_08 AC 35 ms
52,408 KB
testcase_09 AC 36 ms
52,468 KB
testcase_10 AC 895 ms
76,588 KB
testcase_11 AC 383 ms
76,268 KB
testcase_12 AC 263 ms
75,980 KB
testcase_13 AC 158 ms
76,308 KB
testcase_14 AC 280 ms
76,324 KB
testcase_15 AC 302 ms
76,456 KB
testcase_16 AC 434 ms
76,408 KB
testcase_17 AC 495 ms
76,456 KB
testcase_18 AC 195 ms
76,204 KB
testcase_19 AC 231 ms
75,980 KB
testcase_20 AC 533 ms
76,896 KB
testcase_21 AC 210 ms
76,100 KB
testcase_22 AC 239 ms
76,128 KB
testcase_23 AC 199 ms
76,460 KB
testcase_24 AC 220 ms
76,020 KB
testcase_25 AC 221 ms
76,096 KB
testcase_26 AC 313 ms
76,204 KB
testcase_27 AC 236 ms
76,308 KB
testcase_28 AC 459 ms
76,976 KB
testcase_29 AC 225 ms
76,300 KB
testcase_30 AC 232 ms
75,856 KB
testcase_31 AC 432 ms
76,520 KB
testcase_32 AC 400 ms
76,508 KB
testcase_33 AC 450 ms
76,872 KB
testcase_34 AC 448 ms
76,384 KB
testcase_35 AC 300 ms
76,204 KB
testcase_36 AC 300 ms
76,524 KB
testcase_37 AC 463 ms
76,680 KB
testcase_38 AC 480 ms
76,724 KB
testcase_39 AC 459 ms
76,936 KB
testcase_40 AC 262 ms
76,456 KB
testcase_41 AC 523 ms
76,460 KB
testcase_42 AC 420 ms
76,472 KB
testcase_43 AC 425 ms
76,748 KB
testcase_44 AC 526 ms
76,796 KB
testcase_45 AC 354 ms
76,808 KB
testcase_46 AC 342 ms
76,720 KB
testcase_47 AC 486 ms
76,392 KB
testcase_48 AC 386 ms
76,536 KB
testcase_49 AC 333 ms
76,492 KB
testcase_50 AC 274 ms
76,596 KB
testcase_51 AC 292 ms
76,332 KB
testcase_52 AC 228 ms
76,324 KB
testcase_53 AC 446 ms
76,524 KB
testcase_54 AC 356 ms
76,660 KB
testcase_55 AC 392 ms
76,400 KB
testcase_56 AC 395 ms
76,260 KB
testcase_57 AC 365 ms
76,668 KB
testcase_58 AC 190 ms
76,492 KB
testcase_59 AC 428 ms
76,668 KB
testcase_60 AC 396 ms
76,732 KB
testcase_61 AC 317 ms
76,328 KB
testcase_62 AC 448 ms
77,008 KB
testcase_63 AC 517 ms
76,856 KB
testcase_64 AC 555 ms
76,444 KB
testcase_65 AC 318 ms
76,788 KB
testcase_66 AC 469 ms
76,948 KB
testcase_67 AC 258 ms
76,688 KB
testcase_68 AC 292 ms
76,308 KB
testcase_69 AC 239 ms
76,312 KB
testcase_70 AC 280 ms
76,516 KB
testcase_71 AC 307 ms
76,464 KB
testcase_72 AC 423 ms
76,388 KB
testcase_73 AC 369 ms
76,360 KB
testcase_74 AC 474 ms
76,472 KB
testcase_75 AC 474 ms
76,748 KB
testcase_76 AC 424 ms
76,684 KB
testcase_77 AC 433 ms
76,320 KB
testcase_78 AC 513 ms
76,728 KB
testcase_79 AC 455 ms
76,200 KB
testcase_80 AC 523 ms
76,744 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

# https://tjkendev.github.io/procon-library/python/geometry/circles_associated_with_triangle.html
def circumcircle(P1, P2, P3):
    x1, y1 = P1
    x2, y2 = P2
    x3, y3 = P3
    a = 2 * (x1 - x2)
    b = 2 * (y1 - y2)
    p = x1**2 - x2**2 + y1**2 - y2**2
    c = 2 * (x1 - x3)
    d = 2 * (y1 - y3)
    q = x1**2 - x3**2 + y1**2 - y3**2
    det = a * d - b * c
    if det == 0:
        return None
    x = d * p - b * q
    y = a * q - c * p
    if det < 0:
        x = -x
        y = -y
        det = -det
    r2 = (x - x1 * det) ** 2 + (y - y1 * det) ** 2
    return x, y, r2, det


Q = int(input())
xa, ya, xb, yb, xc, yc = map(int, input().split())
res = circumcircle((xa, ya), (xb, yb), (xc, yc))
if res is None:
    Ps = [(xa, ya), (xb, yb), (xc, yc)]
    Ps.sort()
    x1, y1 = Ps[0]
    x2, y2 = Ps[2]
    cx = x1 + x2
    cy = y1 + y2
    det = 2
    r2 = (cx - 2 * x1) ** 2 + (cy - 2 * y1) ** 2
else:
    cx, cy, r2, det = res

for i in range(3):
    xa, xb, xc = xc, xa, xb
    ya, yb, yc = yc, ya, yb
    xx = xa + xb
    yy = ya + yb
    dd = 2
    rr = (xx - 2 * xa) ** 2 + (yy - 2 * ya) ** 2
    rr2 = (xx - 2 * xc) ** 2 + (yy - 2 * yc) ** 2
    if rr2 > rr:
        continue

    if rr * det * det < r2 * dd * dd:
        cx = xx
        cy = yy
        r2 = rr
        det = dd


for _ in range(Q):
    x, y = map(int, input().split())
    x *= det
    y *= det
    d2 = (x - cx) ** 2 + (y - cy) ** 2
    if d2 <= r2:
        print("Yes")
    else:
        print("No")
0