結果

問題 No.2968 Final MIGISITA Strike
ユーザー ねしんねしん
提出日時 2024-11-03 16:13:28
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,580 ms / 2,525 ms
コード長 3,994 bytes
コンパイル時間 291 ms
コンパイル使用メモリ 82,552 KB
実行使用メモリ 114,276 KB
最終ジャッジ日時 2024-11-06 22:00:26
合計ジャッジ時間 42,512 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
53,248 KB
testcase_01 AC 41 ms
53,248 KB
testcase_02 AC 41 ms
52,864 KB
testcase_03 AC 39 ms
53,120 KB
testcase_04 AC 40 ms
53,376 KB
testcase_05 AC 40 ms
52,864 KB
testcase_06 AC 40 ms
52,864 KB
testcase_07 AC 40 ms
52,864 KB
testcase_08 AC 41 ms
52,736 KB
testcase_09 AC 40 ms
52,992 KB
testcase_10 AC 706 ms
94,632 KB
testcase_11 AC 768 ms
95,876 KB
testcase_12 AC 347 ms
82,916 KB
testcase_13 AC 506 ms
87,456 KB
testcase_14 AC 314 ms
82,632 KB
testcase_15 AC 1,276 ms
104,824 KB
testcase_16 AC 1,255 ms
104,664 KB
testcase_17 AC 1,109 ms
100,016 KB
testcase_18 AC 744 ms
90,364 KB
testcase_19 AC 988 ms
96,096 KB
testcase_20 AC 1,045 ms
101,820 KB
testcase_21 AC 482 ms
83,500 KB
testcase_22 AC 708 ms
87,412 KB
testcase_23 AC 1,343 ms
107,840 KB
testcase_24 AC 1,393 ms
109,700 KB
testcase_25 AC 1,580 ms
114,276 KB
testcase_26 AC 675 ms
87,956 KB
testcase_27 AC 935 ms
98,724 KB
testcase_28 AC 714 ms
86,812 KB
testcase_29 AC 857 ms
92,368 KB
testcase_30 AC 178 ms
78,208 KB
testcase_31 AC 177 ms
77,824 KB
testcase_32 AC 193 ms
78,848 KB
testcase_33 AC 175 ms
78,336 KB
testcase_34 AC 704 ms
88,432 KB
testcase_35 AC 619 ms
87,172 KB
testcase_36 AC 246 ms
78,936 KB
testcase_37 AC 582 ms
88,576 KB
testcase_38 AC 550 ms
87,780 KB
testcase_39 AC 551 ms
86,136 KB
testcase_40 AC 766 ms
90,140 KB
testcase_41 AC 614 ms
81,584 KB
testcase_42 AC 432 ms
83,660 KB
testcase_43 AC 706 ms
86,812 KB
testcase_44 AC 963 ms
93,108 KB
testcase_45 AC 409 ms
80,392 KB
testcase_46 AC 1,165 ms
97,784 KB
testcase_47 AC 816 ms
90,120 KB
testcase_48 AC 763 ms
91,972 KB
testcase_49 AC 78 ms
75,648 KB
testcase_50 AC 464 ms
84,416 KB
testcase_51 AC 677 ms
88,612 KB
testcase_52 AC 618 ms
88,440 KB
testcase_53 AC 945 ms
87,544 KB
testcase_54 AC 371 ms
82,388 KB
testcase_55 AC 1,066 ms
97,392 KB
testcase_56 AC 568 ms
87,480 KB
testcase_57 AC 986 ms
95,180 KB
testcase_58 AC 785 ms
92,828 KB
testcase_59 AC 573 ms
86,392 KB
testcase_60 AC 463 ms
84,888 KB
testcase_61 AC 410 ms
83,028 KB
testcase_62 AC 915 ms
92,448 KB
testcase_63 AC 766 ms
92,752 KB
testcase_64 AC 42 ms
52,864 KB
testcase_65 AC 220 ms
83,720 KB
testcase_66 AC 237 ms
85,760 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#遭難者さんのコード、入力方式だけ変更
import sys

input = sys.stdin.readline
# import pypyjit
# pypyjit.set_param('max_unroll_recursion=-1')


def egcd(a, b):
    if a == 0:
        return b, 0, 1
    else:
        g, y, x = egcd(b % a, a)
        return g, x - (b // a) * y, y


def crt(bList, mList):
    r, m = 0, 1
    for i in range(len(bList)):
        d, x, y = egcd(m, mList[i])
        if (bList[i] - r) % d != 0:
            return [0, -1]
        tmp = (bList[i] - r) // d * x % (mList[i] // d)
        r += m * tmp
        m *= mList[i] // d
    return (r, m)


MOD = 998244353
H, W, N, M, S, C, X, Y = map(int, input().split())
Y = W - Y
H2, W2 = 2 * H, 2 * W
A = [0]*N
B = [0]*M
x, y = [0] * N, [0] * N
for i in range(N):
    x[i], y[i], A[i] = map(int, input().split())
    y[i] = W - y[i]
a, b = [0] * M, [0] * M
for i in range(M):
    a[i], b[i], B[i] = map(int, input().split())
    b[i] = W - b[i]
_, syu = crt([0, 0], [H2, W2])


def xy(t):
    xx, yy = (X + t) % H2, (Y + t) % W2
    return min(xx, H2 - xx), min(yy, W2 - yy)


L = []
for i in range(N):
    res, rres = crt([(x[i] - X) % H2, (y[i] - Y) % W2], [H2, W2])
    if rres != -1:
        L.append((res, A[i]))
    res, rres = crt([(-x[i] - X) % H2, (y[i] - Y) % W2], [H2, W2])
    if rres != -1:
        L.append((res, A[i]))
    res, rres = crt([(x[i] - X) % H2, (-y[i] - Y) % W2], [H2, W2])
    if rres != -1:
        L.append((res, A[i]))
    res, rres = crt([(-x[i] - X) % H2, (-y[i] - Y) % W2], [H2, W2])
    if rres != -1:
        L.append((res, A[i]))
for i in range(M):
    res, rres = crt([(a[i] - X) % H2, (b[i] - Y) % W2], [H2, W2])
    if rres != -1:
        L.append((res, -B[i]))
    res, rres = crt([(-a[i] - X) % H2, (b[i] - Y) % W2], [H2, W2])
    if rres != -1:
        L.append((res, -B[i]))
    res, rres = crt([(a[i] - X) % H2, (-b[i] - Y) % W2], [H2, W2])
    if rres != -1:
        L.append((res, -B[i]))
    res, rres = crt([(-a[i] - X) % H2, (-b[i] - Y) % W2], [H2, W2])
    if rres != -1:
        L.append((res, -B[i]))
L.append((0, 0))
L.append((syu, 0))
L.sort()
n = len(L)
now_hp = 0
min_hp = 0
ss, tt = 0, 0
for i in range(1, n):
    real_hp = now_hp + L[i][1] - C * ((L[i][0] + X) // H + (L[i][0] + Y) // W)
    if S + real_hp - L[i][1] <= 0 or S + real_hp <= 0:
        if S + real_hp - L[i][1] > 0:
            x, y = xy(L[i][0])
            assert L[i][1] < 0
            tt += 1
            print(x, W - y, ss % MOD, tt % MOD)
            exit()
        assert i != 0
        ok, ng = L[i - 1][0], L[i][0]
        while ng - ok != 1:
            vs = (ok + ng) >> 1
            real_hp2 = now_hp - C * ((vs + X) // H + (vs + Y) // W)
            if S + real_hp2 <= 0:
                ng = vs
            else:
                ok = vs
        x, y = xy(ng)
        print(x, W - y, ss % MOD, tt % MOD)
        exit()
    now_hp += L[i][1]
    min_hp = min(min_hp, real_hp)
    if L[i][1] > 0:
        ss += 1
    if L[i][1] < 0:
        tt += 1
now_hp -= C * ((syu + X) // H + (syu + Y) // W)
if now_hp >= 0:
    assert S + min_hp > 0
    print(-1)
    exit()
k = (S + min_hp - 1) // -now_hp + 1
assert k > 0
S += k * now_hp
now_hp = 0
ss, tt = ss * k % MOD, tt * k % MOD
for i in range(1, n):
    real_hp = now_hp + L[i][1] - C * ((L[i][0] + X) // H + (L[i][0] + Y) // W)
    if S + real_hp - L[i][1] <= 0 or S + real_hp <= 0:
        if S + real_hp - L[i][1] > 0:
            x, y = xy(L[i][0])
            assert L[i][1] < 0
            tt += 1
            print(x, W - y, ss % MOD, tt % MOD)
            exit()
        ok, ng = L[i - 1][0], L[i][0]
        while ng - ok != 1:
            vs = (ok + ng) >> 1
            real_hp2 = now_hp - C * ((vs + X) // H + (vs + Y) // W)
            if S + real_hp2 <= 0:
                ng = vs
            else:
                ok = vs
        x, y = xy(ng)
        print(x, W - y, ss % MOD, tt % MOD)
        exit()
    now_hp += L[i][1]
    if L[i][1] > 0:
        ss += 1
    if L[i][1] < 0:
        tt += 1
exit(1)
0