#include <bits/stdc++.h>
using namespace std;
long long mod = 1000000007;
long long gcdExtended(long long a, long long b, long long *x, long long *y);
long long modInverse(long long a, long long m) {
  long long x, y;
  long long g = gcdExtended(a, m, &x, &y);
  long long res = (x % m + m) % m;
  return res;
}
long long gcdExtended(long long a, long long b, long long *x, long long *y) {
  if (a == 0) {
    *x = 0, *y = 1;
    return b;
  }
  long long x1, y1;
  long long gcd = gcdExtended(b % a, a, &x1, &y1);
  *x = y1 - (b / a) * x1;
  *y = x1;
  return gcd;
}
long long power(long long x, long long y) {
  long long temp;
  if (y == 0) return 1;
  temp = (power(x, y / 2)) % mod;
  if (y % 2 == 0) {
    long long p = ((temp % mod) * (temp % mod)) % mod;
    return p;
  } else {
    long long p = ((temp % mod) * (temp % mod)) % mod;
    long long c = x % mod;
    long long s = (x * p) % mod;
    return s;
  }
}
int main() {
  long long a, b, n, x;
  cin >> a >> b >> n >> x;
  long long q = power(a, n) % mod;
  long long p = (q * x) % mod;
  if (a > 1)
    q = (((((q - 1) % mod) * (modInverse(a - 1, mod)) % mod) % mod) *
         (b % mod)) %
        mod;
  else
    q = ((n % mod) * (b % mod)) % mod;
  cout << (p + q) % mod << "\n";
}
