結果
| 問題 | No.62 リベリオン(Extra) | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2024-05-05 08:48:23 | 
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 321 ms / 5,000 ms | 
| コード長 | 1,168 bytes | 
| コンパイル時間 | 105 ms | 
| コンパイル使用メモリ | 12,672 KB | 
| 実行使用メモリ | 11,008 KB | 
| 最終ジャッジ日時 | 2024-11-27 04:06:11 | 
| 合計ジャッジ時間 | 1,364 ms | 
| ジャッジサーバーID (参考情報) | judge2 / judge4 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 2 | 
| other | AC * 3 | 
ソースコード
import itertools
import math
def inv(a, mod):
    b, x, y=mod, 1, 0
    while b:
        t = a // b
        a -= t * b
        x -= t * y
        a, b = b, a
        x, y = y, x
    return x % mod
Q=int(input())
for _ in range(Q):
    w, h, d, mx, my, hx, hy, vx, vy=map(int, input().split())
    if vx < 0:
        vx = -vx
        hx, mx = w - hx, w - mx
    if vy < 0:
        vy = -vy
        hy, my = h - hy, h - my
    a = 2 * h * vx
    b = 2 * w * vy
    g = math.gcd(a, b)
    a //= g
    b //= g
    for cx, cy in itertools.product((mx, 2 * w - mx), (my, 2 * h - my)):
        c = (cx - hx) * vy - (cy - hy) * vx
        if c % g != 0:
            continue
        c //= g
        if b:
            r = inv(a, b) * c % b
            x1 = -2 * h * r - cy + hy
            if (x1 + 2 * h * b - 1) // (2 * h * b) <= (x1 + d * vy) // (2 * h * b):
                print("Hit")
                break
        else:
            s = inv(b, a) * c % a
            x1 = -2 * w * s - cx + hx
            if (x1 + 2 * w * a - 1) // (2 * w * a) <= (x1 + d * vx)//(2 * w * a):
                print("Hit")
                break
    else:
        print("Miss")
        
            
            
            
        