#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#include<queue>
#include<vector>
#include <bitset>
#include <cmath>
#include <limits>
#include <iostream>
using namespace std;
#define INF 100000
#define MAX 100000
#define MOD 1000000007
typedef long long ll;
typedef pair<int,int> P;
typedef pair< pair<int,int>, int> p; 
#define bit(n,k) ((n>>k)&1) /*nのk bit目*/

int myplus[30];

int main(){
    ll fact[2000000];
    fact[0]=1;
    for(ll i=1;i<2000000;i++){
        fact[i]=i*fact[i-1];
        fact[i]%=MOD;
    }

    int d,ld,rd;
    ll l,r,k,s[30],t[30];
    cin>>d>>l>>r>>k;

    s[1]=1LL; t[1]=1LL;
    for(int i=2;i<=d;i++){
        s[i]=s[i-1]*2;
        t[i]=t[i-1]*2+1;
    }
    for(int i=1;i<=d;i++){
        if(s[i]<=l && l<=t[i]) ld=i;
        if(s[i]<=r && r<=t[i]) rd=i;
    }
    bool flag=false;
    int temp=rd-ld;
    ll ans=1;
    int unko=1;
    for(int i=0;i<ld;i++){
        if(i>1) unko*=2;
        if(temp+i*2==k){
            flag=true;
            if(rd==ld){ //正しく機能
                if(i==1){
                    ans*=2*(t[rd]-s[rd]+1)/2;
                }else{
                    ans*=(t[rd]-s[rd]+1)*unko;
                }
            }else{
                //unko/=2;
                if(unko<1) unko=1;
                ans*=unko*(t[rd]-s[rd]+1);
            }
        } 
    }
    if(!flag){
        cout<<0<<endl;
        return 0;
    }
    myplus[ld]--; myplus[rd]--;
    for(int i=1;i<=d;i++){
        ans*=fact[t[i]-s[i]+1+myplus[i]];
        ans%=MOD;
    }
    cout<<ans<<endl;
}