mod = 10 ** 9 + 7

def f(x):
    v = 1
    for i in range(1, 22):
        if 2 ** (i - 1) <= x <= 2 ** i - 1:
            v = i
    return v
            
def g(x):
    v = 1
    for i in range(1, x + 1):
        v *= i
        v %= mod
    return v

D, L, R, K = map(int, input().split())
ans = 1
for i in range(D):
    ans *= g(2 ** i)
    ans %= mod
    
dL, dR = f(L), f(R)
if dL >= dR:
    dL, dR = dR, dL
K -= dR - dL
# print("test", K, dR, dL, ans)
if K < 0 or K % 2:
    print(0)
    exit()
H = K // 2
if H >= dL:
    print(0)
    exit()
if dL == dR:
    ans *= pow(2 ** (dL - 1) - 1, mod - 2, mod)
    ans *= 2 ** max(0, (H - 1))
    ans %= mod
else:
    ans *= pow(2 ** (dL - 1), mod - 2, mod)
    ans *= 2 ** max(0, (H - 1))
    ans %= mod
print(ans)