#include <cstdio>
#include <algorithm>
using namespace std;
int A, B, C, D;
long long dptable[3001][3001];
const long long mod = 998244353;
long long dp(int i, int j) {
	if (i == A && j == B) return 1;
	if (dptable[i][j] != 0) return dptable[i][j];
	if (i == A) return dptable[i][j] = i * dp(i, j - 1) % mod;
	if (j == B) return dptable[i][j] = j * dp(i - 1, j) % mod;
	return dptable[i][j] = (i * dp(i, j - 1) + j * dp(i - 1, j) - (i - 1) * (j - 1) * dp(i - 1, j - 1)) % mod;
}
int main() {
	scanf("%d %d %d %d", &A, &B, &C, &D);
	printf("%lld", (dp(C, D) + mod) % mod);
	return 0;
}