#include<bits/stdc++.h>
#define int long long
using namespace std;
constexpr int mod=998244353;
int dp[3000][3000];
signed main(){
    int a,b,c,d;
    cin>>a>>b>>c>>d;
    for(int i=0;i<3000;i++){
        for(int j=0;j<3000;j++){
            if(i==0){
                if(j==0)
                    dp[i][j]=1;
                else
                    dp[i][j]=dp[i][j-1]*a%mod;
            }else{
                if(j==0)
                    dp[i][j]=dp[i-1][j]*b%mod;
                else
                    dp[i][j]=(dp[i][j-1]*(a+i)+dp[i-1][j]*(b+j)-dp[i-1][j-1]*(a+i-1)*(b+j-1)%mod+mod)%mod;
            }
        }
    }
    cout<<dp[c-a][d-b]<<endl;
}