#include using namespace std; template class Matrix { protected: vector< vector > _mat; int _col, _row; const bool usingModulo = true; const Numeric modulo = 1000000007; // using modulo void resize(int col, int row) { _col = col; _row = row; _mat.assign(col, vector(row, 0)); } public: Matrix(int col, int row) : _col(col), _row(row) { _mat.assign(col, vector(row, 0)); } Matrix(int col, int row, int num) : _col(col), _row(row) { _mat.assign(col, vector(row, num)); } int col() const { return _col; } int row() const { return _row; } vector> mat() const {return _mat; } vector &operator[](int i) { return _mat[i]; } Numeric &operator() (int i, int j) { return _mat[i][j]; } Matrix &operator=(Matrix x) { this->_mat = x.mat(); return *this; } Matrix operator+(Matrix &x) { if (this->col() != x.col() || this->row() != x.row()) { cout << "matrices are different sizes .. "; this->size(false); cout << "+"; x.size(true); exit(0); } Matrix ret(this->col(), this->row()); for (int i = 0; i < this->col(); ++i) { for (int j = 0; j < this->row(); ++j) { ret[i][j] = _mat[i][j] + x[i][j]; if (usingModulo) ret[i][j] %= modulo; } } return ret; } Matrix &operator+=(Matrix &x) { if (this->col() != x.col() || this->row() != x.row()) { cout << "matrices are different sizes .. "; this->size(false); cout << "+"; x.size(true); exit(0); } for (int i = 0; i < this->col(); ++i) { for (int j = 0; j < this->row(); ++j) { _mat[i][j] += x[i][j]; if (usingModulo) _mat[i][j] %= modulo; } } return *this; } Matrix operator-(Matrix &x) { if (this->col() != x.col() || this->row() != x.row()) { cout << "matrices are different sizes .. "; this->size(false); cout << "-"; x.size(true); exit(0); } Matrix ret(this->col(), this->row()); for (int i = 0; i < this->col(); ++i) { for (int j = 0; j < this->row(); ++j) { ret[i][j] = _mat[i][j] - x[i][j]; if (usingModulo) ret[i][j] %= modulo; } } return ret; } Matrix &operator-=(Matrix &x) { if (this->col() != x.col() || this->row() != x.row()) { cout << "matrices are different sizes .. "; this->size(false); cout << "-"; x.size(true); exit(0); } for (int i = 0; i < this->col(); ++i) { for (int j = 0; j < this->row(); ++j) { _mat[i][j] -= x[i][j]; if (usingModulo) _mat[i][j] %= modulo; } } return *this; } Matrix operator*(Matrix &x) { if (this->row() != x.col()) { cout << "cannot multiply matrices .. "; this->size(false); cout << "*"; x.size(true); exit(0); } Matrix ret(this->col(), x.row()); for (int i = 0; i < this->col(); ++i) { for (int j = 0; j < x.row(); ++j) { for (int k = 0; k < this->row(); ++k) { ret[i][j] += _mat[i][k] * x[k][j]; if (usingModulo) ret[i][j] %= modulo; } } } return ret; } Matrix &operator*(Numeric k) { Matrix ret(this->col(), this->row()); for (int i = 0; i < this->col(); ++i) { for (int j = 0; j < this->row(); ++j) { ret(i, j) = _mat[i][j] * k; if (usingModulo) ret(i, j) %= modulo; } } return ret; } Matrix &operator*=(Matrix &x) { if (this->row() != x.col()) { cout << "cannot multiply matrices .. "; this->size(false); cout << "*"; x.size(true); exit(0); } Matrix ret(this->col(), x.row()); for (int i = 0; i < this->col(); ++i) { for (int j = 0; j < x.row(); ++j) { for (int k = 0; k < this->row(); ++k) { ret[i][j] += _mat[i][k] * x[k][j]; if (usingModulo) ret[i][j] %= modulo; } } } this->resize(this->col(), x.row()); for (int i = 0; i < this->col(); ++i) { for (int j = 0; j < x.row(); ++j) { _mat[i][j] = ret[i][j]; } } return *this; } Matrix &operator*=(Numeric k) { Matrix ret(this->col(), this->row()); for (int i = 0; i < this->col(); ++i) { for (int j = 0; j < this->row(); ++j) { _mat[i][j] *= k; if (usingModulo) _mat %= modulo; } } return *this; } Numeric sum() { Numeric ret = 0; for (int i = 0; i < this->col(); ++i) { for (int j = 0; j < this->row(); ++j) { ret += _mat[i][j]; if (usingModulo) ret %= modulo; } } return ret; } void print() { cout << "-- print "; this->size(false); cout << " --" << endl; for (int i = 0; i < _col; ++i) { for (int j = 0; j < _row; ++j) { if (j != 0) cout << "\t"; cout << _mat[i][j]; } cout << endl; } cout << "------------------" << endl; } void print(string str) { cout << "-- print \"" << str << "\" "; this->size(false); cout << " --" << endl; for (int i = 0; i < _col; ++i) { for (int j = 0; j < _row; ++j) { if (j != 0) cout << "\t"; cout << _mat[i][j]; } cout << endl; } cout << "------------------" << endl; } void size (bool endline = true) const { cout << "(" << this->col() << ", " << this->row() << ")"; if (endline) cout << endl; } }; int main() { long long a, b, n; cin >> a >> b >> n; Matrix ans(2, 1), x(2, 2); ans[0][0] = 1; ans[1][0] = 0; x[0][0] = a; x[0][1] = b; x[1][0] = 1; x[1][1] = 0; while (n) { if (n & 1) ans = x * ans; x *= x; n /= 2; } cout << ans[1][0] << endl; return 0; }