//#define _AOJ_ /*vvv> zzzzzI .---. zzuzuI .vgggg&,. +++++= dAC:|I .WbbWo JMM9^```?TMB` ..&gNNg,. gggggggJ, qgggggggg] (&&&&&&&&[ c+OA&J, (&&&&&&+J, .cJeAA&-. (&&&&&&&&x .&AA&=-. +++++= dTqk|I Xpbpbp JM#` (M#^ ?MMp MM| +TMN. JMF ' |yk ` dVY 7Vk, Vy XV cVf ?Y! JM V$ ` +++++= dcf:|I Xppppp dMN .MM+ .MM MM| MM] JMMMMMM+ |@tqkoh) ,y0 (V$ yyyyyyyV7 VV JMWyZWr TWVVVVW&, ++++++ d7qk|0 Xppppp ^HMN, _.db WMm, .MMF MM| ..MM` JMF . |yk .WV&. .XW' yy 4yn. jyn +. JM #S `++++` ?ZZZX= ?WWWW= -THMMMMH9^ (TMMMMM9! MMMMMMM"" JMMMMMMMME |UU. ?TUUUUY= UU. (UU- ^7TUUUV7! JUUUUUUUU 7TUNKO*/ //Ricty Diminished #include "bits/stdc++.h" using namespace std;typedef long long lint;typedef vector liv; //#define rep(i,n) for(int i=0;i=0;--i) #define _rrep2(i,n) for(lint i=n-1;i>=0;--i) #define _rrep3(i,a,b) for(lint i=b-1;i>=a;--i) #define _rrep4(i,a,b,c) for(lint i=a+(b-a-1)/c*c;i>=a;i-=c) template istream& operator>>(istream& is,vector& vec); template istream& operator>>(istream& is,array& vec); template istream& operator>>(istream& is,pair& p); template ostream& operator<<(ostream& os,vector& vec); template ostream& operator<<(ostream& os,pair& p); template istream& operator>>(istream& is,vector& vec){ for(T& x: vec) is>>x;return is; } template istream& operator>>(istream& is,pair& p){ is>>p.first;is>>p.second;return is; } //template //ostream& operator<<(ostream& os,vector& vec){ os< //ostream& operator<<(ostream& os,deque& deq){ os< ostream& operator<<(ostream& os,pair& p){ os< inline void in(Head&& head,Tail&&... tail){ cin>>head;in(move(tail)...); } template inline bool out(T t){ cout< inline bool out(Head head,Tail... tail){ cout< name(size);in(name) #define fi e.first #define se e.second #define YES(c) cout<<((c)?"YES\n":"NO\n"),0 #define Yes(c) cout<<((c)?"Yes\n":"No\n"),0 #define o(p) cout< md_tmp class mint{ using u64=uint_fast64_t; public: u64 a; constexpr mint(const u64 x=0) noexcept: a(x%md){} constexpr u64 &value() noexcept{ return a; } constexpr const u64 &value() const noexcept{ return a; } constexpr mint operator+(const mint rhs) const noexcept{ return mint(*this)+=rhs; } constexpr mint operator-(const mint rhs) const noexcept{ return mint(*this)-=rhs; } constexpr mint operator*(const mint rhs) const noexcept{ return mint(*this)*=rhs; } constexpr mint operator^(const lint rhs) const noexcept{ return mint(*this)^=rhs; } constexpr mint operator/(const mint rhs) const noexcept{ return mint(*this)/=rhs; } constexpr mint &operator+=(const mint rhs) noexcept{ a+=rhs.a; if(a>=md)a-=md; return *this; } constexpr mint &operator-=(const mint rhs) noexcept{ if(aa; while(exp){ if(exp&1)*this*=base; base*=base; exp>>=1; } return *this; } constexpr mint &operator/=(const mint rhs) noexcept{ a=(*this*(rhs^(md-2))).a; return *this; } }; md_tmp istream& operator>>(istream& os,mint& m){ os>>m.a,m.a%=md; return os; } md_tmp ostream& operator<<(ostream& os,const mint& m){ return os< ncr(lint n,lint r){ if(n(0); mintncr_res=1,ncr_div=1; rep(r)ncr_res*=(n-i),ncr_div*=(r-i); return ncr_res/ncr_div; } #ifndef _AOJ_ mint<> operator""m(unsigned long long n){ return mint<>(n); } mint<998244353> operator""m9(unsigned long long n){ return mint<998244353>(n); } mint<1000003> operator""m3(unsigned long long n){ return mint<1000003>(n); } #endif using mi=mint<>; //P class P{ public:lint f,s; constexpr P(lint a,lint b):f(a),s(b){};constexpr P():f(0),s(0){}; }; istream& operator>>(istream& is,P& p){ is>>p.f>>p.s;return is; } ostream& operator<<(ostream& os,const P& p){ return os<(const P& l,const P& r){ return(l.f-r.f?l.f>r.f:l.s>r.s); } struct C{ lint f,s,t; }; bool operator<(const C& l,const C& r){ return l.t(const C& l,const C& r){ return l.t>r.t; } #define linf 1152921504606846976 #define inf linf//INT_MAX #define MAXN 200100 #define md_1e9_7 1000000007 #define md_998244353 998244353 //#define mod md_1e9_7 const int d4[5]={0,1,0,-1,0}; template< class T > struct Matrix{ vector< vector< T > > A; Matrix(){} Matrix(size_t n,size_t m): A(n,vector< T >(m,0)){} Matrix(size_t n): A(n,vector< T >(n,0)){}; size_t height() const{ return (A.size()); } size_t width() const{ return (A[0].size()); } inline const vector< T > &operator[](int k) const{ return (A.at(k)); } inline vector< T > &operator[](int k){ return (A.at(k)); } static Matrix I(size_t n){ Matrix mat(n); for(int i = 0; i < n; i++) mat[i][i] = 1; return (mat); } Matrix &operator+=(const Matrix &B){ size_t n = height(),m = width(); assert(n == B.height() && m == B.width()); for(int i = 0; i < n; i++) for(int j = 0; j < m; j++) (*this)[i][j] += B[i][j]; return (*this); } Matrix &operator-=(const Matrix &B){ size_t n = height(),m = width(); assert(n == B.height() && m == B.width()); for(int i = 0; i < n; i++) for(int j = 0; j < m; j++) (*this)[i][j] -= B[i][j]; return (*this); } Matrix &operator*=(const Matrix &B){ size_t n = height(),m = B.width(),p = width(); assert(p == B.height()); vector< vector< T > > C(n,vector< T >(m,0)); for(int i = 0; i < n; i++) for(int j = 0; j < m; j++) for(int k = 0; k < p; k++) C[i][j] = (C[i][j] + (*this)[i][k] * B[k][j]); A.swap(C); return (*this); } Matrix &operator^=(long long k){ Matrix B = Matrix::I(height()); while(k > 0){ if(k & 1) B *= *this; *this *= *this; k >>= 1LL; } A.swap(B.A); return (*this); } Matrix operator+(const Matrix &B) const{ return (Matrix(*this) += B); } Matrix operator-(const Matrix &B) const{ return (Matrix(*this) -= B); } Matrix operator*(const Matrix &B) const{ return (Matrix(*this) *= B); } Matrix operator^(const long long k) const{ return (Matrix(*this) ^= k); } friend ostream &operator<<(ostream &os,Matrix &p){ size_t n = p.height(),m = p.width(); for(int i = 0; i < n; i++){ os << "["; for(int j = 0; j < m; j++){ os << p[i][j] << (j + 1 == m ? "]\n" : ","); } } return (os); } T determinant(){ Matrix B(*this); assert(width() == height()); T ret = 1; for(int i = 0; i < width(); i++){ int idx = -1; for(int j = i; j < width(); j++){ if(B[j][i] != 0) idx = j; } if(idx == -1) return (0); if(i != idx){ ret *= -1; swap(B[i],B[idx]); } ret *= B[i][i]; T vv = B[i][i]; for(int j = 0; j < width(); j++){ B[i][j] /= vv; } for(int j = i + 1; j < width(); j++){ T a = B[j][i]; for(int k = 0; k < width(); k++){ B[j][k] -= B[i][k] * a; } } } return (ret); } }; int main(){ mint<> a,b; lint n; in(a,b,n); Matrix> A;A.A={{a,1m},{b,0m}}; if(n<=1)return o(n); o((A^(n-1))[0][0]); } //sub-EOF