#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const int mod=998244353;

int a,b,c,d;
ll dp[3005][3005];

int main(){

    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);

    cin>>a>>b>>c>>d;
    dp[a][b]=1;
    for(int i=a+1;i<=c;i++)dp[i][b]=dp[i-1][b]*b%mod;
    for(int i=b+1;i<=d;i++)dp[a][i]=dp[a][i-1]*a%mod;
    for(int i=a+1;i<=c;i++){
        for(int j=b+1;j<=d;j++){
            dp[i][j]=((dp[i-1][j]*j+dp[i][j-1]*i-dp[i-1][j-1]*(i-1)*(j-1))%mod+mod)%mod;
        }
    }
    cout<<dp[c][d]<<endl;

    return 0;
}