結果

問題 No.61 リベリオン
ユーザー mkawa2mkawa2
提出日時 2020-01-09 18:17:30
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 203 ms / 5,000 ms
コード長 2,341 bytes
コンパイル時間 149 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 10,752 KB
最終ジャッジ日時 2024-10-13 16:16:57
合計ジャッジ時間 1,141 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 29 ms
10,624 KB
testcase_01 AC 28 ms
10,752 KB
testcase_02 AC 27 ms
10,752 KB
testcase_03 AC 198 ms
10,752 KB
testcase_04 AC 203 ms
10,752 KB
testcase_05 AC 28 ms
10,624 KB
権限があれば一括ダウンロードができます

ソースコード

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