a, b, c = map(int, input().split())
k = int(input())
mod = 10 ** 9 + 7

x = pow(2, k, mod - 1)

ans = 1
for i in [a, b, c]:
    ans *= pow(i, x, mod)
    ans %= mod
print(ans)