from sympy import Symbol,solve
[a,b,c,d,n] = list(map(int,input().split()))
MOD = 10**9+7
x = Symbol('x')
[x1,x2,x3,x4] = solve(x**4-x**3-x-1,x)
q = Symbol('q')
w = Symbol('w')
e = Symbol('e')
r = Symbol('r')
d = solve([-a+q*x1+w*x2+e*x3+r*x4,-b+q*(x1**2)+w*(x2**2)+e*(x3**2)+r*(x4**2),
-c+q*(x1**3)+w*(x2**3)+e*(x3**3)+r*(x4**3),-d+q*(x1**4)+w*(x2**4)+e*(x3**4)+r*(x4**4)],[q,w,e,r])
print(int(d[q]*(x1**n)+d[w]*(x2**n)+d[e]*(x3**n)+d[r]*(x4**n))%MOD)