import sys

sys.setrecursionlimit(10 ** 6)
input = sys.stdin.readline

def main():
    md = 10 ** 9 + 7
    d, l, r, k = map(int, input().split())
    mx = pow(2, 20) + 5
    fac = [1] * mx
    for i in range(1, mx):
        fac[i] = (fac[i - 1] * i) % md

    d0 = d1 = 0
    while l:
        l //= 2
        d0 += 1
    while r:
        r //= 2
        d1 += 1
    c = k - (d1 - d0)
    if c < 0 or c % 2:
        print(0)
        exit()
    h = k - c // 2
    if c > 0: h -= 1
    ans = 1
    for i in range(d):
        if i + 1 == d1:
            w = pow(2, h, md)
            if d0 == d1:
                pw = pow(2, i, md)
                ans = (ans * pw * w * fac[pw - 2]) % md
            else:
                ans = (ans * fac[pow(2, i, md) - w] * w) % md
        else:
            ans = (ans * fac[pow(2, i, md)]) % md
    print(ans)

main()