d,L,R,k = map(int,input().split())

MOD = 10**9+7
N = 2**d

fact = [1]*(N+1)
rfact = [1]*(N+1)
r = 1
for i in range(1, N+1):
    fact[i] = r = r * i % MOD
rfact[N] = r = pow(fact[N], MOD-2, MOD)
for i in range(N, 0, -1):
    rfact[i-1] = r = r * i % MOD

# nPk (mod MOD) を求める
def perm(n, k):
    return fact[n] * rfact[n-k] % MOD

rui = []
val = 1
for i in range(d):
    rui.append(val)
    if L >= val:
        lidx = i
    if R >= val:
        ridx = i
    val *= 2

ans = 1

f = False
if lidx == ridx:
    f = True
    if k % 2 != 0 or k > 2*lidx:
        print(0)
        exit()
else:
    if k % 2 != (ridx-lidx) % 2 or k > lidx+ridx:
        print(0)
        exit()
print(lidx,ridx)
if f == True:
    for i in range(d):
        if i == lidx:
            tmp = pow(2,i,MOD)
            if tmp == 2:
                ans *= tmp*1*perm(tmp-2,tmp-2)
            else:
                ans *= tmp*2*perm(tmp-2,tmp-2)
            ans %= MOD
            continue
        tmp = pow(2,i,MOD)
        ans *= perm(tmp,tmp)
        ans %= MOD
    print(ans%MOD)
else:
    for i in range(d):
        if i == lidx:
            tmp = pow(2,i,MOD)
            aa = k // 2
            ans *= (2**(aa-1))*perm(tmp-1,tmp-1)
            ans %= MOD
            continue
        tmp = pow(2,i,MOD)
        ans *= perm(tmp,tmp)
        ans %= MOD
    print(ans%MOD)