#include<bits/stdc++.h>
using namespace std;
template <typename A, typename B>
string to_string(pair<A, B> p);

template <typename A, typename B, typename C>
string to_string(tuple<A, B, C> p);

template <typename A, typename B, typename C, typename D>
string to_string(tuple<A, B, C, D> p);

string to_string(const string &s)
{
    return '"' + s + '"';
}

string to_string(const char *s)
{
    return to_string((string)s);
}

string to_string(bool b)
{
    return (b ? "true" : "false");
}

string to_string(vector<bool> v)
{
    bool first = true;
    string res = "{";
    for (int i = 0; i < static_cast<int>(v.size()); i++)
    {
        if (!first)
        {
            res += ", ";
        }
        first = false;
        res += to_string(v[i]);
    }
    res += "}";
    return res;
}

template <size_t N>
string to_string(bitset<N> v)
{
    string res = "";
    for (size_t i = 0; i < N; i++)
    {
        res += static_cast<char>('0' + v[i]);
    }
    return res;
}

template <typename A>
string to_string(A v)
{
    bool first = true;
    string res = "{";
    for (const auto &x : v)
    {
        if (!first)
        {
            res += ", ";
        }
        first = false;
        res += to_string(x);
    }
    res += "}";
    return res;
}

template <typename A, typename B>
string to_string(pair<A, B> p)
{
    return "(" + to_string(p.first) + ", " + to_string(p.second) + ")";
}

template <typename A, typename B, typename C>
string to_string(tuple<A, B, C> p)
{
    return "(" + to_string(get<0>(p)) + ", " + to_string(get<1>(p)) + ", " + to_string(get<2>(p)) + ")";
}

template <typename A, typename B, typename C, typename D>
string to_string(tuple<A, B, C, D> p)
{
    return "(" + to_string(get<0>(p)) + ", " + to_string(get<1>(p)) + ", " + to_string(get<2>(p)) + ", " + to_string(get<3>(p)) + ")";
}

void debug_out() { cerr << endl; }

template <typename Head, typename... Tail>
void debug_out(Head H, Tail... T)
{
    cerr << " " << to_string(H);
    debug_out(T...);
}

typedef long long ll;
typedef pair<int, int> pii;
#define mp make_pair
#define fst first
#define snd second
#define fr(i, a, b) for(int i=a; i<b; i++)
template<class T = int> using V = vector<T>;
template<class T = int> using VV = V< V<T> >;
ll mod = 998244353;
ll dp[3005][3005];
int a, b, cc, d;

ll f(int r, int c) {
    ll &res = dp[r][c];
    if(res != -1) return res;
    if(r==a && c==b) return res = 1;
    res = 0;
    if(!r || !c) return res;
    res += f(r-1, c) * c %mod;
    res += f(r, c-1) * r %mod;
    res -= f(r-1,c-1) * (r-1) * (c-1) %mod;
    res += mod;
    res %= mod;
    // if(r==a) {
    //     for(int i=c; i>b; i--) if(r>1)res = res * (r-1) %mod;
    // } else if(c == b) {
    //     for(int i=r; i>a; i--) if(c>1) res = res *(c-1) % mod;
    // } else {
    //     res = res * (r+c-1) * (r+c-2) / 2 % mod * f(r-1, c-1) % mod;
    // }
    // for(int i=r; i>a; i--) {
    //     for(int j=c; j>b; j--) {
    //         res = res * 
    //     }
    // }
    return res;
}

int main() {
    cin.tie(nullptr); ios::sync_with_stdio(false);
    cin >> a >> b >> cc >> d;
    memset(dp, -1, sizeof dp);
    dp[a][b] = 1;
    cout << f(cc, d);
}


Measured memory_bytes: 74620