#pragma GCC optimize("O3") #include using namespace std; using ll=long long; //typedef pair P; template using V=vector; #define fi first #define se second #define all(v) (v).begin(),(v).end() const ll inf=(1e18); //const ll mod=998244353; const ll mod=1000000007; ll gcd(ll a,ll b) {return b ? gcd(b,a%b):a;} ll lcm(ll c,ll d){return c/gcd(c,d)*d;} struct __INIT{__INIT(){cin.tie(0);ios::sync_with_stdio(false);cout< bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (a>b) { a=b; return 1; } return 0; } struct mint{ using ull=unsigned long long int; ull v; mint(ll vv=0){s(vv%mod+mod);} mint& s(ull vv){ v=vv>=1ll; } return res; } mint inv()const{return pow(mod-2);} //拡張ユークリッドの互除法 /* mint inv()const{ int x,y; int g=extgcd(v,mod,x,y); assert(g==1); if(x<0)x+=mod; return mint(x); }*/ friend ostream& operator<<(ostream&os,const mint&val){ return os<(const mint&val)const{return v>val.v;} }; template struct Matrix{ vector> A; Matrix() {} Matrix(size_t n,size_t m):A(n,vector(m,0)) {} Matrix(size_t n):A(n,vector(n,0)) {} size_t height() const{return (A.size());} size_t width() const{return (A[0].size());} inline const vector &operator[](int k)const{return (A.at(k));} inline vector &operator[](int k){return (A.at(k));} //単位行列 static Matrix I(size_t n){ Matrix mat(n); for(int i=0;i> C(n,vector(m,0)); for(int i=0;i0){ 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 &mat){ size_t n=mat.height(),m=mat.width(); for(int i=0;i A(2),x(2,1); ll a,b,n; cin>>a>>b>>n; A.A[0][0]=mint(a); A.A[0][1]=mint(b); A.A[1][0]=mint(1); x.A[0][0]=mint(1); if(n==0){ cout<<0<<"\n"; return 0; } if(n==1){ cout<<1<<"\n"; return 0; } A^=(n-1); A*=x; cout<