d,l,r,k=map(int,input().split())

TWO=[2**i for i in range(23)]

mod=10**9+7

FACT=[1]
for i in range(1,(1<<d)+1):
    FACT.append(FACT[-1]*i%mod)

for i in range(23):
    if l<TWO[i]:
        lh=i
        break

for i in range(23):
    if r<TWO[i]:
        rh=i
        break

if l==1:
    if rh-lh==k:
        ANS=1
        for i in range(d):
            ANS=ANS*FACT[1<<i]%mod

        print(ANS)

    else:
        print(0)
    
    
elif (lh+rh)%2==k%2 and abs(lh-rh)<=k<=lh+rh:
    ANS=1
    for i in range(d):
        ANS=ANS*FACT[1<<i]%mod

    #print(ANS)

    kei=(k-(rh-lh))//2

    #print(kei)

    if lh==rh:
        ANS=ANS*pow((1<<(lh-1))-1,mod-2,mod)*pow(2,kei-1,mod)%mod

        print(ANS)

    else:
        ANS=ANS*pow((1<<(lh-1)),mod-2,mod)*pow(2,max(kei,1)-1,mod)%mod

        print(ANS)

else:
    print(0)