#include #define REP(i,n) for(int i=0,i##_len=(n);i=0) _num%=_mod; else _num+=(1-(_num+1)/_mod)*_mod; return *this; } ll _mpow(ll x, ll n){ //x^n(mod) ←普通にpow(x,n)では溢れてしまうため,随時mod計算 2分累乗法だから早い ll ans = 1; while(n != 0){ if(n&1) ans = ans*x % _mod; x = x*x % _mod; n = n >> 1; } return ans; } ll imod(ll n){return _mpow(n , _mod-2);} public: mint(){ _num = 0;_mod=mod; } mint(ll num){ _mod = mod; _num = (num+(1LL<<25)*mod) % mod; } mint(ll num,ll M){ _mod=M;_num=(num+(1LL<<25)*mod)%_mod; } mint(const mint &cp){_num=cp._num;_mod=cp._mod;} mint operator= (const ll x){ return set(x); } mint operator+ (const ll x){ return mint(_num + (x % _mod) , _mod); } mint operator- (const ll x){ return mint(_num - (x % _mod), _mod); } mint operator* (const ll x){ return mint(_num * (x % _mod) , _mod); } mint operator/ (ll x){ return mint(_num * imod(x) , _mod);} mint operator+=(const ll x){ return set(_num + (x % _mod)); } mint operator-=(const ll x){ return set(_num - (x % _mod)); } mint operator*=(const ll x){ return set(_num * (x % _mod)); } mint operator/=(ll x){ return set(_num * imod(x));} mint operator+ (const mint &x){ return mint(_num + x._num , _mod); } mint operator- (const mint &x){ return mint(_num - x._num , _mod);} mint operator* (const mint &x){ return mint(_num * x._num , _mod); } mint operator/ (mint x){ return mint(_num * imod(x._num) , _mod);} mint operator+=(const mint &x){ return set(_num + x._num); } mint operator-=(const mint &x){ return set(_num - x._num); } mint operator*=(const mint &x){ return set(_num * x._num); } mint operator/=(mint x){ return set(_num * imod(x._num));} bool operator<(const mint &x)const{return _num(const mint &x)const{return _num>x._num;} friend mint operator+(ll x,const mint &m){return mint(m._num + (x % m._mod) , m._mod);} friend mint operator-(ll x,const mint &m){return mint( (x % m._mod) - m._num , m._mod);} friend mint operator*(ll x,const mint &m){return mint(m._num * (x % m._mod) , m._mod);} friend mint operator/(ll x,mint m){return mint(m.imod(m._num) * x , m._mod);} explicit operator ll() { return _num; } explicit operator int() { return (int)_num; } friend ostream& operator<<(ostream &os, const mint &x){ os << x._num; return os; } friend istream& operator>>(istream &is, mint &x){ll val; is>>val; x.set(val); return is;} }; template class MAT{ private: int row,col; vector> _A; MAT set(vector> A){_A = A ; return *this;} public: MAT(){ } MAT(int n,int m){ if(n<1 || m<0){cout << "err Matrix::Matrix" < a(col,0.0); _A.push_back(a); } // 値の初期化 if(m==0) REP(i,n) _A[i][i]=1.0; } MAT(const MAT &cp){_A=cp._A;row=cp.row;col=cp.col;} T* operator[] (int i){return _A[i].data();} MAT operator= (vector> x) {return set(x);} MAT operator+ (MAT x) { if(row!=x.row || col!=x.col){ cout << "err Matrix::operator+" <>=1; } return r; } MAT operator+= (MAT x) { if(row!=x.row || col!=x.col){ cout << "err Matrix::operator+=" <>(){return _A;} friend ostream &operator<<(ostream &os,const MAT &x){ REP(i,x.row) REP(j,x.col) os<>(istream &is,MAT &x){REP(i,x.row) REP(j,x.col) is>>x._A[i][j];return is;} int size_row(){return row;} int size_col(){return col;} MAT transpose(){ MAT r(col,row); REP(i,col) REP(j,row) r[i][j]=_A[j][i]; return r; } MAT inverse(){ T buf; MAT inv_a(row,0); vector> a=_A; //掃き出し法 REP(i,row){ buf=1/a[i][i]; REP(j,row){ a[i][j]*=buf; inv_a[i][j]*=buf; } REP(j,row){ if(i!=j){ buf=a[j][i]; REP(k,row){ a[j][k]-=a[i][k]*buf; inv_a[j][k]-=inv_a[i][k]*buf; } } } } return inv_a; } // O( n^3 ). int rank() { vector> A=_A; const int n = row, m = col; int r = 0; for(int i = 0; r < n && i < m; ++i) { int pivot = r; for(int j = r+1; j < n; ++j) if(fabs(A[j][i]) > fabs(A[pivot][i])) pivot = j; swap(A[pivot], A[r]); if(fabs(A[r][i]) < eps) continue; for (int k = m-1; k >= i; --k) A[r][k] /= A[r][i]; rep(j,r+1,n) rep(k,i,m) A[j][k] -= A[r][k] * A[j][i]; ++r; } return r; } }; int main(){ MAT x(2,1),A(2,2); ll a,b,n;cin>>a>>b>>n; A[0][0]=a;A[0][1]=b;A[1][0]=1; x[0][0]=1;x[1][0]=0; x=(A^n)*x; cout<