x, y, z, w = map(int, input().split())
mod = 998244353

dif = x + y - z - w - 1
fact = [1] * (dif + 1)
inv = [1] * (dif + 1)
finv = [1] * (dif + 1)
for i in range(2, dif + 1):
    fact[i] = fact[i - 1] * i % mod
    inv[i] = mod - inv[mod % i] * (mod // i) % mod
    finv[i] = finv[i - 1] * inv[i] % mod
def comb(x, y):
    return fact[x] * finv[y] % mod * finv[x - y] % mod

ans = 1
for i in range(x, z, -1):
    ans *= i
    ans %= mod
for i in range(y, w, -1):
    ans *= i
    ans %= mod
if z == 0:
    ans *= comb(dif, y - w)
else:
    ans *= comb(dif, x - z)
ans %= mod
print(ans)