MOD = 10**9+7
N = 2**21
fac,p2 = [1]*N,[1]*N
for i in range(1,N):
    fac[i] = fac[i-1]*i%MOD
    p2[i] = p2[i-1]*2%MOD
d,l,r,k = map(int,input().split())
a,b = l.bit_length(), r.bit_length()
c = b-a
ans = 1-(k-c)%2
for i in range(1,d+1):
    if b != i != a: ans *= fac[1<<i-1]
k = (k-c)//2
ans *= p2[c] if k==0 else p2[c+k-1] if 0 <= k <= a-1 else 0
ans *= [fac[(1<<a-1)]*fac[(1<<b-1)-1],p2[a-1]*fac[(1<<a-1)-2]][a==b] 
print(ans%MOD)