H,W,N,P = map(int,input().split())
s = set()
for _ in range(N):
    x,y = map(int,input().split())
    s.add((x-1,y-1))
dp = [[0] * W for _ in range(H)]
memo = [[0] * W for _ in range(H)]
stack = [(0,0)]
while stack:
    x,y = stack.pop()
    for i,j in [(1,0),(0,1)]:
        if x + i < H and y + j < W:
            if memo[x + i][y + j] == 0:
                memo[x + i][y + j] = 1
                stack.append((x + i,y + j))
            if (x + i,y + j) in s:
                t = dp[x][y] + 1
            else:
                t = dp[x][y]
            if t > dp[x + i][y + j]:
                dp[x + i][y + j] = t
u = dp[-1][-1]
mod = 998244353
v = H + W - 3 - u
r = pow(P,mod-2,mod)
tmp = pow((1-r)%mod,v,mod)*pow((1-r * 2)%mod,u,mod)
ans = (1-tmp)%mod
print(ans)