#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 3000 + 10, mod = 998244353;
ll f[N][N][2];
void add(ll &x,ll y){
    (x += y) %= mod;
}
int main(){
    ios::sync_with_stdio(0);
    int a,b,c,d; cin>>a>>b>>c>>d;
    f[a][b][0] = 1;
    for(int i=a;i<=c;i++){
        for(int j=b;j<=d;j++){
            add(f[i+1][j][0],f[i][j][0]*j);
            add(f[i][j+1][1],f[i][j][0]*i);
            add(f[i+1][j][0],f[i][j][1]);
            add(f[i][j+1][1],f[i][j][1]*i);
        }
    }
    ll ans = 0;
    add(ans,f[c][d][0]);
    add(ans,f[c][d][1]);
    cout<<ans<<endl;
}
