#include #include #include #define llint long long #define mod 1000000007 using namespace std; struct Matrix{ int h, w; vector mat; Matrix(){ h = w = 0; } Matrix(int h, int w){ this->h = h, this->w = w; mat.resize(h*w); } llint& at(int i, int j){ return mat[w*(i-1)+(j-1)]; } static Matrix ident(int size){ Matrix ret(size, size); for(int i = 1; i <= size; i++) ret.at(i, i) = 1; return ret; } Matrix operator*(Matrix& ope){ Matrix ret(h, ope.w); for(int i = 1; i <= h; i++){ for(int j = 1; j <= ope.w; j++){ for(int k = 1; k <= w; k++){ ret.at(i, j) += at(i, k) * ope.at(k, j) % mod; ret.at(i, j) %= mod; } } } return ret; } Matrix pow(llint n){ if(n == 0) return ident(h); if(n % 2){ return pow(n-1) * (*this); } else{ Matrix tmp = pow(n/2); return tmp * tmp; } } void print(){ for(int i = 1; i <= h; i++){ for(int j = 1; j <= w; j++){ cout << at(i, j) << " "; } cout << endl; } } }; llint a, b, n; int main(void) { ios::sync_with_stdio(0); cin.tie(0); cin >> a >> b >> n; Matrix mat(3, 3); for(int i = 1; i <= 3; i++){ for(int j = 1; j <= 3; j++){ mat.at(i, j) = 0; } } mat.at(1, 1) = a, mat.at(1, 2) = b, mat.at(2, 1) = mat.at(3, 2) = 1; Matrix pmat = mat.pow(n); llint ans = pmat.at(3, 1)*a%mod + pmat.at(3, 2); ans %= mod; cout << ans << endl; return 0; }