MOD = 10**9+7
N = 2**20+2
fac = [1]*N
for i in range(2,N):
    fac[i] = fac[i-1]*i%MOD
p2 = [1]*N
for i in range(1,N):
    p2[i] = p2[i-1]*2%MOD

d,l,r,k = map(int,input().split())
a = l.bit_length()
b = r.bit_length()
ans = 1
for i in range(1,d+1):
    if i != a and i != b:
        ans *= fac[1<<i-1]
        ans %= MOD

if a==b:
    ans *= p2[a-1]
    ans %= MOD
    ans *= fac[(1<<a-1)-2]
    ans %= MOD
    if k%2: ans = 0
    k //= 2
    if 0 <= k <= a-1:
        ans *= p2[k-1]
        ans %= MOD
    else:
        ans = 0

else:
    ans *= fac[(1<<a-1)]
    ans %= MOD
    ans *= fac[(1<<b-1)-1]
    ans %= MOD
    c = b-a
    if (k-c)%2: ans = 0
    k = (k-c)//2
    if k==0:
        ans *= p2[c]
        ans %= MOD
    elif 0 <= k <= a-1:
        ans *= p2[c+k-1]
        ans %= MOD
    else:
        ans = 0

print(ans)