#include<iostream> #include<vector> #include<algorithm> using namespace std; using ll=long long; #define rep(i,n) for(int i=0;i<n;i++) #define rrep(i,n) for(int i=(n)-1;i>=0;i--) #define all(v) v.begin(),v.end() #define rall(v) v.rbegin(),v.rend() template<class T> bool chmax(T &a, T b){if (a < b){a = b;return true;} else return false;} template<class T> bool chmin(T &a, T b){if (a > b){a = b;return true;} else return false;} const int MAX=1e6,MOD=998244353; vector<ll>fac(MAX),finv(MAX),inv(MAX); void set_fac(){ fac[0]=fac[1]=1; finv[0]=finv[1]=1; inv[1]=1; for (int i=2;i<MAX;i++){ fac[i]=fac[i-1]*i%MOD; inv[i]=MOD-inv[MOD%i]*(MOD/i)%MOD; finv[i]=finv[i-1]*inv[i]%MOD; } } ll cmb(int n,int r){ if (r<0||n<r)return 0; return fac[n]*(finv[r]*finv[n-r]%MOD)%MOD; } int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); set_fac(); ll n,m,a,b; cin>>n>>m>>a>>b; if(a*(n-1)>b){ cout<<0<<"\n"; return 0; } m-=a*(n-1); b-=a*(n-1); ll ans=0; for(int i=0;i<=b;i++){ ll cnt=cmb(i+n-2,i)*(m-i)%MOD; ans+=cnt; ans%=MOD; } rep(i,n)ans=ans*(i+1)%MOD; cout<<ans<<"\n"; }