from collections import *
import heapq
import bisect

INF = float("inf")
MOD = 998244353
mod = 998244353


def fur(n, r):
    p, q = 1, 1
    for i in range(r):
        p = p * (n - i) % mod
        q = q * (i + 1) % mod
    return p * pow(q, mod - 2, mod) % mod


X, Y, Z, W = map(int, input().split())
if W == 0:
    ans = 1
    for i in range(1, X - Z + Y - W):
        ans *= i
        ans %= MOD
    ans *= Y - W
    ans %= MOD
    ans *= fur(Y, Y - W)
    ans %= MOD
    ans *= fur(X, X - Z)
    ans %= MOD
else:
    ans = 1
    for i in range(1, X - Z + Y - W):
        ans *= i
        ans %= MOD
    ans *= X - Z
    ans %= MOD
    ans *= fur(Y, Y - W)
    ans %= MOD
    ans *= fur(X, X - Z)
    ans %= MOD

print(ans)