#include<bits/stdc++.h>
using namespace std;

constexpr long long mod = 1e9+7;
long long fac[1 << 20];

void fill_fac(){
    fac[0] = 1;
    for (int i = 1; i < (1 << 20); ++i){
        fac[i] = fac[i - 1] * i % mod;
    }
}

int find_depth(int x){
    int d = 1;
    while (x > (1 << (d - 1))){
        x -= 1 << (d - 1);
        ++d;
    }
    return d;
}

int solve_core(int d, int dr, int dl, int k){
    if (dl < dr) swap(dl, dr);
    int s = (k - (dl - dr)) / 2;
    long long ans = 0;
    if (s == 0){
        if (dl == dr) return 0;
        ans = (1 << (dl - dr)) % mod;
        ans = ans * fac[(1 << (dl - 1)) - 1] % mod;
        for (int i = 1; i <= d; ++i){
            if (i == dl) continue;
            ans = ans * fac[1 << (i - 1)] % mod;
        }
        return ans;
    }else if (dl == dr){
        ans = 1 << (dl - 1);
        ans = ans * (1 << (s - 1)) % mod;
        ans = ans * fac[(1 << (dl - 1)) - 2] % mod;
        for (int i = 1; i <= d; ++i){
            if (i == dl) continue;
            ans = ans * fac[1 << (i - 1)] % mod;
        }
        return ans;
    }else{
        ans = 1 << (s - 1);
        ans = ans * fac[(1 << (dr - 1)) - 1] % mod;
        for (int i = 1; i <= d; ++i){
            if (i == dr) continue;
            ans = ans * fac[1 << (i - 1)] % mod;
        }
        return ans;
    }
}

int solve(int d, int l, int r, int k){
    int dl = find_depth(l);
    int dr = find_depth(r);
    if (abs(dl - dr) > k){
        return 0;
    }else if (dl + dr < k){
        return 0;
    }else if ((k - abs(dl - dr)) % 2){
        return 0;
    }else if ((k - abs(dl - dr)) / 2 > min(dl, dr) - 1){
        return 0;
    }
    return solve_core(d, dl, dr, k);
}

int main(){
    fill_fac();
    int d, l, r, k;
    cin >> d >> l >> r >> k;
    cout << solve(d, l, r, k) << endl;
    return 0;
}