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

int main(){
  long long a,b,c,d,i,j;
  cin >> a >> b >> c >> d;
  vector<vector<long long>> dp(c+1,vector<long long>(d+1));
  long long p=998244353;
  
  dp.at(a).at(b)=1;
  for(i=a+1;i<=c;i++){
    dp.at(i).at(b)=(dp.at(i-1).at(b)*b)%p;
  }
  for(j=b+1;j<=d;j++){
    dp.at(a).at(j)=(dp.at(a).at(j-1)*a)%p;
  }
  
  for(i=a+1;i<=c;i++){
    for(j=b+1;j<=d;j++){
      dp.at(i).at(j)=(dp.at(i-1).at(j)*j+dp.at(i).at(j-1)*i-dp.at(i-1).at(j-1)*(i-1)*(j-1))%p;
    }
  }
  
  if(dp.at(c).at(d)<0){
    dp.at(c).at(d) += p;
  }
  
  cout << dp.at(c).at(d) << endl;
}