#include #include using namespace std; using ll = long long; constexpr ll MOD = 1000000007; using Mat = vector>; Mat mul(const Mat& m1, const Mat& m2) { Mat ans(2, vector(2, 0LL)); for (int i = 0; i < 2; i++) { for (int j = 0; j < 2; j++) { for (int k = 0; k < 2; k++) { (ans[i][j] += m1[i][k] * m2[k][j] % MOD) %= MOD; } } } return ans; } Mat unit() { Mat ans(2, vector(2, 0)); ans[0][0] = 1; ans[1][1] = 1; return ans; } Mat power(const Mat& m, const ll n) { if (n == 0) { return unit(); } if (n % 2 == 1) { return mul(power(m, n - 1), m); } else { const Mat pp = power(m, n / 2); return mul(pp, pp); } } int main() { ll B, C, D; cin >> B >> C >> D; B %= MOD, C %= MOD; Mat mat(2, vector(2, 0)); mat[0][0] = C; mat[0][1] = B * C % MOD; mat[1][1] = 1; const Mat ans = power(mat, D); cout << ans[0][1] << endl; return 0; }