#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1e5+3,MOD=1e9+7;
int n,k,l1,l2,ans1;
LL ans2,ans3,ans4;
LL p[N],q[N];
LL f[N][22],g[N][22];
void upd(LL &x,LL y){
	(x+=y)%=MOD;
}
LL qpow(LL x,int y=MOD-2){
	LL cnt=1;
	while(y){
		if(y&1) cnt=cnt*x%MOD;
		x=x*x%MOD;
		y>>=1;
	}
	return cnt;
}
LL C(int x,int y){
	if(x<0||y<0||x<y) return 0;
	return p[x]*q[y]%MOD*q[x-y]%MOD;
}
int Main()
{
    scanf("%d%d%d%d",&n,&k,&l1,&l2);
    p[0]=1;for(int i=1;i<=N-3;i++) p[i]=p[i-1]*i%MOD;
    q[N-3]=qpow(p[N-3]);for(int i=N-3;i>=1;i--) q[i-1]=q[i]*i%MOD; 
    for(int i=1;i<=k;i++) if(k%i==0) ans1++;
    f[1][1]=1,g[1][1]=1;
    for(int i=2;i<=n;i++){
    	for(int j=1;j*j<=i;j++)
    	    if(i%j==0){
    	    	for(int l=0;l<20;l++)
    	    	    upd(f[i][l+1],f[j][l]),upd(g[i][l+1],g[j][l]+f[j][l]*i);
    	    	if(j==1||j*j==i) continue;
    	    	for(int l=0;l<20;l++)
    	    	    upd(f[i][l+1],f[i/j][l]),upd(g[i][l+1],g[i/j][l]+f[i/j][l]*i);
			} 
	}
	for(int i=1;i<=n;i++)
	    for(int j=0;j<20;j++)
	        upd(f[0][j+1],f[i][j]),upd(g[0][j+1],g[i][j]);
	for(int i=1;i<20;i++){
		//cout<<i<<' '<<f[0][i]<<' '<<g[0][i]<<' '<<C(l1,i-1)<<endl;
		upd(ans2,f[0][i]*C(l1,i-1));
		if(i<=l2+1) upd(ans3,f[0][i]);
		upd(ans4,g[0][i]);
	}
	printf("%d %lld %lld %lld\n",ans1,ans2,ans3,ans4);
	return 0;
}
int main()
{
    int Case=1;
    while(Case--) Main();
	return 0;
}