結果

問題 No.62 リベリオン(Extra)
ユーザー mkawa2mkawa2
提出日時 2020-01-09 18:20:30
言語 Python3
(3.11.6 + numpy 1.26.0 + scipy 1.11.3)
結果
TLE  
実行時間 -
コード長 2,341 bytes
コンパイル時間 226 ms
コンパイル使用メモリ 11,000 KB
実行使用メモリ 8,100 KB
最終ジャッジ日時 2023-08-15 01:22:14
合計ジャッジ時間 8,150 ms
ジャッジサーバーID
(参考情報)
judge14 / judge11
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 18 ms
8,036 KB
testcase_01 AC 17 ms
8,100 KB
testcase_02 AC 180 ms
8,060 KB
testcase_03 AC 184 ms
8,064 KB
testcase_04 TLE -
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys

sys.setrecursionlimit(10 ** 6)
int1 = lambda x: int(x) - 1
p2D = lambda x: print(*x, sep="\n")
def II(): return int(sys.stdin.readline())
def MI(): return map(int, sys.stdin.readline().split())
def LI(): return list(map(int, sys.stdin.readline().split()))
def LLI(rows_number): return [LI() for _ in range(rows_number)]

def gcd(a, b):
    a, b = abs(a), abs(b)
    while b: a, b = b, a % b
    return a

def main():
    q = II()
    for _ in range(q):
        w, h, d, gx, gy, sx, sy, vx, vy = MI()
        # vx,vyが負の時、対称移動で正負を逆転し、平行移動で第1象限に戻す
        if vx < 0: gx, sx, vx = w - gx, w - sx, - vx
        if vy < 0: gy, sy, vy = h - gy, h - sy, - vy
        # vxかvyが0のときは、別処理
        if vx == 0:
            print("Hit") if sx == gx and (sy <= gy <= sy + d * vy or sy <= 2 * h - gy <= sy + d * vy) else print("Miss")
            continue
        if vy == 0:
            print("Hit") if sy == gy and (sx <= gx <= sx + d * vx or sx <= 2 * w - gx <= sx + d * vx) else print("Miss")
            continue
        # vx,vyが互いに素でないとき、gcdで割って、時間をgcd倍する
        g = gcd(vx, vy)
        vx, vy, d = vx // g, vy // g, d * g
        # 弾丸を反射させるのではなく、辺を軸とした線対称な長方形が無限に並んでいると考える
        # つまりターゲットも対称な位置に無限にいる
        # x,yを独立に考える
        # 弾丸が通る格子点のx座標はsx+vx*i (0<=i<=d)で表される
        # ターゲットのx座標は2*w*j-gx , 2*w*j+gxで表される
        # Hitするのはsx+vx*i=2*w*j+-gxのとき
        # jは何でもいいので、2*wでmodをとるとsx+vx*i≡+-gx (mod 2*w)
        # 移項してvx*i≡+-gx-sx (mod 2*w)
        # yでも同様にして、どちらも満たすような最小のiがd以下だったらHit
        xr1, xr2 = (gx - sx) % (2 * w), (-gx - sx) % (2 * w)
        yr1, yr2 = (gy - sy) % (2 * h), (-gy - sy) % (2 * h)
        for i in range(1, min(d, 2 * h * w // gcd(h, w)) + 1):
            xl, yl = (vx * i) % (2 * w), (vy * i) % (2 * h)
            if (xl - xr1) * (xl - xr2) == 0 and (yl - yr1) * (yl - yr2) == 0:
                print("Hit")
                break
        else:
            print("Miss")

main()
0