#include typedef uint64_t u64; typedef int64_t i64; using namespace std; template void scan_(vector& a){ for(auto& v:a){ cin>>v; } } template void scan_(T& a){ cin>>a; } void scan(){} template void scan(T& n,Args&... args){ scan_(n),scan(args...); } template void println(T n){ cout< void println(T n,Args... args){ cout< void println(vector v){ for(size_t i=0;i struct modint{ u64 val; modint(i64 val_=0):val((val_%i64(mod)+i64(mod))%i64(mod)){} modint operator-(){ return (val==0)?0:mod-val; } modint operator+(modint rhs){ return modint(*this)+=rhs; } modint operator-(modint rhs){ return modint(*this)-=rhs; } modint operator*(modint rhs){ return modint(*this)*=rhs; } modint operator/(modint rhs){ return modint(*this)/=rhs; } modint operator^(i64 rhs){ return modint(*this)^=rhs; } modint &operator+=(modint rhs){ val+=rhs.val,val-=((val>=mod)?mod:0); return (*this); } modint &operator-=(modint rhs){ val+=((val>=1; } return (*this)=res; } bool operator==(modint rhs){ return val==rhs.val; } bool operator!=(modint rhs){ return val!=rhs.val; } friend std::ostream &operator<<(std::ostream& os,modint x){ return os<<(x.val); } friend std::istream &operator>>(std::istream& is,modint& x){ u64 t; is>>t,x=t; return is; } }; template struct matrix{ vector> val; int height,width; matrix(int H,int W,T init=0):val(H,vector (W,init)),height(H),width(W){} matrix(vector> val_):val(val_),height(val_.size()),width(val_[0].size()){} vector &operator[](int i){ return val[i]; } matrix operator+(matrix rhs){ return matrix(*this)+=rhs; } matrix operator*(matrix rhs){ return matrix(*this)*=rhs; } matrix operator^(i64 rhs){ return matrix(*this)^=rhs; } matrix operator+=(matrix rhs){ for(int i=0;i res(height,rhs.width,0); for(int i=0;iheight;++i){ for(int j=0;jwidth;++k){ res[i][j]+=(*this)[i][k]*rhs[k][j]; } } } return (*this)=res; } matrix operator^=(i64 k){ int n=(*this).height; matrix res=identity(n),now=(*this); while(k>0){ if(k%2==1){ res*=now; } now*=now; k/=2; } return (*this)=res; } matrix identity(int n){ matrix res(n,n); for(int i=0;i mint; mint A,B; scan(A,B); int N; scan(N); matrix C({{0,0,1},{1,0,0},{0,0,0}}),D(vector>{{1,0,0},{0,1,0},{A,B,1}}),E(vector>{{1},{1},{0}}); for(int i=0;i