#include <bits/stdc++.h>
using namespace std;
long int modpow(long int a, long int n, long int p) {
  long int b=n,disit=0;
  while (b>0) {
    b/=2;
    disit++;
  }
  long int c=1;
  vector<int>N(disit),M(disit);
  for (int i=0; i<disit; i++) {
    N[i]=n/c%2;
    c*=2;
  }
  long int d=a%p;
  for (int i=0; i<disit; i++) {
    M[i]=d;
    d=d*d;
    d%=p;
  }
  long int ans=1;
  for (int i=0; i<disit; i++) {
    if (N[i]==1) {
      ans*=N[i]*M[i];
      ans%=p;
    }
  }
  return ans;
}

int main() {
    long int A,B,C,K;
    cin >> A >> B >> C >> K;
    long int p=1000000007;
    long int r=modpow(2,K,p-1);
    long int a=modpow(A,r,p);
    long int b=modpow(B,r,p);
    long int c=modpow(C,r,p);
    cout << a*b%p*c%p << endl;
}