// #pragma GCC target("avx2")  // CPU 処理並列化
// #pragma GCC optimize("O3")  // CPU 処理並列化
// #pragma GCC optimize("unroll-loops")  // 条件処理の呼び出しを減らす
#include<stdio.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<deque>
#include<stack>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<bitset>
#include<stdlib.h>
#include<cassert>
#include<time.h>
#include<bitset>
#include<numeric>
using namespace std;
const long long mod=1000000007;
const long long inf=mod*mod;
const long long d2=(mod+1)/2;
const double EPS=1e-10;
const double INF=1e+10;
const double PI=acos(-1.0);
const int C_SIZE = 3121000;
long long fact[C_SIZE];
long long finv[C_SIZE];
long long inv[C_SIZE];
long long Comb(int a,int b){
 	if(a<b||b<0)return 0;
 	return fact[a]*finv[b]%mod*finv[a-b]%mod;
}
void init_C(int n){
	fact[0]=finv[0]=inv[1]=1;
	for(int i=2;i<n;i++){
		inv[i]=(mod-(mod/i)*inv[mod%i]%mod)%mod;
	}
	for(int i=1;i<n;i++){
		fact[i]=fact[i-1]*i%mod;
		finv[i]=finv[i-1]*inv[i]%mod;
	}
}
long long pw(long long a,long long b){
	if(a<0LL)return 0;
	if(b<0LL)return 0;
	long long ret=1;
	while(b){
		if(b%2)ret=ret*a%mod;
		a=a*a%mod;
		b/=2;
	}
	return ret;
}
int ABS(int a){return max(a,-a);}
long long ABS(long long a){return max(a,-a);}
double ABS(double a){return max(a,-a);}
int sig(double r) { return (r < -EPS) ? -1 : (r > +EPS) ? +1 : 0; }
// ここから編集しろ

int main(){
	init_C(1100000);
	int a,b,c,d;
	scanf("%d%d%d%d",&a,&b,&c,&d);
	long long ret=0;
	for(int i=0;i<=d;i++){
		for(int x=0;x<2;x++){
			for(int y=0;y<=a-b;y++){
				if(d-(long long)b*i-x*b-(long long)y*i<0)break;
				long long ks=Comb(a,b);
				ks=ks*Comb(a-b,y)%mod;
				long long L=max(0LL,c-(long long)b*i-x*b-(long long)y*i);
				long long R=d-(long long)b*i-x*b-(long long)y*i;
				long long tmp=Comb(R+a,a);
				if(L)tmp=(tmp+mod-Comb(L-1+a,a))%mod;
				ks=ks*tmp%mod;
				if((x+y)%2==0)ret=(ret+ks)%mod;
				else ret=(ret+mod-ks)%mod;
	//			printf("%d %d %d: %lld\n",i,x,y,ks);
			}
		}
	//	printf("%d: %lld\n",i,ret);
	}
	long long sum=0;
	for(int i=c;i<=d;i++){
		sum=(sum+Comb(i+a-1,a-1))%mod;
	}
	ret=(sum+mod-ret)%mod;
	printf("%lld\n",ret);
}