#include<bits/stdc++.h>
using lint=long long;
using real=long double;
int main(){
    std::cin.tie(nullptr);std::ios_base::sync_with_stdio(false);
    std::cout.setf(std::ios_base::fixed);std::cout.precision(15);
    lint d,l,r,k;std::cin>>d>>l>>r>>k;
    l=std::numeric_limits<lint>::digits-__builtin_clzll(l);
    r=std::numeric_limits<lint>::digits-__builtin_clzll(r);
    bool kyoen_ng=l==r;
    k-=std::max(l,r)-std::min(l,r);
    l=r=std::min(l,r);
    if(k<0||k==0&&kyoen_ng||k%2||2*l<k){
        std::cout<<0<<'\n';
        return 0;
    }
    lint mod=1'000'000'007;
    lint D=1ll<<d;
    std::vector<lint>fact(D+1,1);
    for(lint i=1;i<=D;i++)fact.at(i)=fact.at(i-1)*i%mod;
    lint ans=1;
    for(lint i=0;i<d-1;i++)ans=ans*fact.at(1ll<<i)%mod;
    for(lint i=1;i<=D/2;i++)if(i!=(1ll<<l)-kyoen_ng)ans=ans*i%mod;
    ans=ans*(k==0?1ll:1ll<<k/2-1)%mod;
    std::cout<<ans<<'\n';
}