#include <bits/stdc++.h>

long long mo=(long long)1e9+7;

long long fac[2000000];
long long ifac[2000000];

long long modpow(long long a,long long n){
  long long ret=1;
  for(;n>0;n>>=1,a=a*a%mo){
    if(n%2==1){
      ret=ret*a%mo;
    }
  }
  return ret;
}

long long inv(long a){
  return modpow(a,mo-2);
}

long long comb(int n,int k){
  return fac[n]*ifac[k]%mo*ifac[n-k]%mo;
}

void build(){
  fac[0]=1;
  for(int i=1;i<1000000;++i){
    fac[i]=fac[i-1]*i%mo;
  }
  for(int i=0;i<1000000;++i){
    ifac[i]=inv(fac[i]);
  }
}


int main(){
  build();
  long long N,M,D1,D2;
  std::cin>>N>>M>>D1>>D2;
  long long ans=0;
  for(int i=0;i<=N-1;++i){
    long long  res=M-1-(long long)D1*(N-1)-(long long)i*(D2-D1+1);
    if(res<0)continue;
    ans=(ans+modpow(mo-1,i)*comb(N-1,i)%mo*comb(N+res,res)%mo)%mo;
  }
  std::cout<<ans<<std::endl;
}