p,q,r=map(int,input().split())
k=int(input())
mod=10**9+7
t=pow(2,k,mod-1)
p=pow(p,t,mod)
q=pow(q,t,mod)
r=pow(r,t,mod)
print((p*q*r%mod))