#include<bits/stdc++.h>
#define int long long
using namespace std;
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }
template<class T>bool chmin(T &a, const T &b) { if (a>b) { a=b; return 1; } return 0; }

const int mod=1e9+7;

signed main()
{
    cin.tie(0);
    ios::sync_with_stdio(false);

    int d,l,r,k; cin>>d>>l>>r>>k;

    if(l>r) swap(l,r);
    l=32-__builtin_clz(l)-1;
    r=32-__builtin_clz(r)-1;
    int lca=-1;
    if(l+r-k>=0 and (l+r-k)%2==0) lca=(l+r-k)/2;
    if(lca==-1 or lca>l or lca>r){
        cout<<0<<endl;
        return 0;
    }
    // cout<<l<<' '<<r<<' '<<lca<<endl;

    vector<int> fact((1<<(d-1))+1);
    fact[0]=1;
    for(int i=1;i<=1<<(d-1);i++) fact[i]=fact[i-1]*i%mod;

    int ans=1;
    if(lca==l){
        ans=ans*(1<<lca)%mod;
        ans=ans*(1<<(r-lca))%mod;
    }else{
        ans=ans*(1<<lca)%mod;
        ans=ans*(1<<(l-lca))%mod;
        ans=ans*(1<<(r-lca-1))%mod;
    }
    // cout<<ans<<endl;
    for(int i=0;i<d;i++){
        int cnt=1<<i;
        if(i==l) cnt--;
        if(i==r) cnt--;
        ans=ans*fact[cnt]%mod;
    }
    cout<<ans<<endl;
    return 0;
}