#define _USE_MATH_DEFINES

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#include <map>
#include <set>
#include <unordered_map> 
#include <unordered_set>
#include <functional>

using namespace std;

typedef pair<long long int, long long int> P;
typedef tuple<int, int, int> T;

long long int INF = 1e18;
long long int MOD = 1e9 + 7;

long long int rev_mod(long long int num){
	long long int pow_num = MOD - 2;
	long long int mod_num = MOD;
	long long int pow_cur = num;
	long long int ret = 1;
	while(pow_num > 0){
		if(pow_num % 2 == 1){
			ret *= pow_cur;
			ret %= mod_num;
		}
		pow_cur *= pow_cur;
		pow_cur %= mod_num;
		pow_num /= 2;
	}
	return ret;
}

int main(){
	long long int d, l, r, k;
	cin >> d >> l >> r >> k;
	long long int n1 = 0, n2 = 0;
	while(l != 0){
		l /= 2;
		n1 += 1;
	}
	while(r != 0){
		r /= 2;
		n2 += 1;
	}
	n1 -= 1;
	n2 -= 1;
	if(n2 - n1 > k || n1 + n2 < k || (n1 + n2) % 2 != k % 2){
		cout << 0 << endl;
		return 0;
	}
	long long int ans = 1;
	for(int i = 0; i < d; i++){
		for(int j = 1; j <= (1 << i); j++){
			ans *= j;
			ans %= MOD;
		}
	}
	int p = (n1 + n2 - k) / 2;
	if(n1 != n2){
		ans *= rev_mod(1 << p);
		ans %= MOD;
	}else{
		ans *= (1 << (n1 - p - 1));
		ans %= MOD;
		ans *= rev_mod((1 << n1) - 1);
		ans %= MOD;
	}
	cout << ans << endl;
	return 0;
}