結果
問題 | 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 itertoolsimport mathdef inv(a, mod):b, x, y=mod, 1, 0while b:t = a // ba -= t * bx -= t * ya, b = b, ax, y = y, xreturn x % modQ=int(input())for _ in range(Q):w, h, d, mx, my, hx, hy, vx, vy=map(int, input().split())if vx < 0:vx = -vxhx, mx = w - hx, w - mxif vy < 0:vy = -vyhy, my = h - hy, h - mya = 2 * h * vxb = 2 * w * vyg = math.gcd(a, b)a //= gb //= gfor cx, cy in itertools.product((mx, 2 * w - mx), (my, 2 * h - my)):c = (cx - hx) * vy - (cy - hy) * vxif c % g != 0:continuec //= gif b:r = inv(a, b) * c % bx1 = -2 * h * r - cy + hyif (x1 + 2 * h * b - 1) // (2 * h * b) <= (x1 + d * vy) // (2 * h * b):print("Hit")breakelse:s = inv(b, a) * c % ax1 = -2 * w * s - cx + hxif (x1 + 2 * w * a - 1) // (2 * w * a) <= (x1 + d * vx)//(2 * w * a):print("Hit")breakelse:print("Miss")