import sys input = sys.stdin.readline 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 = list(map(int, input().split())) B = list(map(int, input().split())) x, y = [0] * N, [0] * N for i in range(N): x[i], y[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] = 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(n): real_hp = now_hp + L[i][1] - C * ((L[i][0] + X) // H + (L[i][0] + Y) // W) if 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, tt) 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 + now_hp - 1) // -min_hp + 1 assert k > 0 S += k * min_hp now_hp = 0 ss, tt = ss * k % MOD, tt * k % MOD for i in range(n): real_hp = now_hp + L[i][1] - C * ((L[i][0] + X) // H + (L[i][0] + Y) // W) if 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) 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] if L[i][1] > 0: ss += 1 if L[i][1] < 0: tt += 1 exit(1)