#define PROBLEM "https://yukicoder.me/problems/1973" #include using namespace std; #define call_from_test #ifndef call_from_test #include using namespace std; #endif //BEGIN CUT HERE struct FastIO{ FastIO(){ cin.tie(0); ios::sync_with_stdio(0); } }fastio_beet; //END CUT HERE #ifndef call_from_test signed main(){ return 0; } #endif #ifndef call_from_test #include using namespace std; #endif //BEGIN CUT HERE template struct Mint{ static constexpr T mod = MOD; T v; Mint():v(0){} Mint(signed v):v(v){} Mint(long long t){v=t%MOD;if(v<0) v+=MOD;} Mint pow(long long k){ Mint res(1),tmp(v); while(k){ if(k&1) res*=tmp; tmp*=tmp; k>>=1; } return res; } static Mint add_identity(){return Mint(0);} static Mint mul_identity(){return Mint(1);} Mint inv(){return pow(MOD-2);} Mint& operator+=(Mint a){v+=a.v;if(v>=MOD)v-=MOD;return *this;} Mint& operator-=(Mint a){v+=MOD-a.v;if(v>=MOD)v-=MOD;return *this;} Mint& operator*=(Mint a){v=1LL*v*a.v%MOD;return *this;} Mint& operator/=(Mint a){return (*this)*=a.inv();} Mint operator+(Mint a) const{return Mint(v)+=a;} Mint operator-(Mint a) const{return Mint(v)-=a;} Mint operator*(Mint a) const{return Mint(v)*=a;} Mint operator/(Mint a) const{return Mint(v)/=a;} Mint operator-() const{return v?Mint(MOD-v):Mint(v);} bool operator==(const Mint a)const{return v==a.v;} bool operator!=(const Mint a)const{return v!=a.v;} bool operator <(const Mint a)const{return v constexpr T Mint::mod; template ostream& operator<<(ostream &os,Mint m){os<>h>>w>>k; using M = Mint; M ans{0}; for(int d=1;d using namespace std; #endif //BEGIN CUT HERE template struct Kitamasa{ using VR = vector; VR cs; vector rs; int m; Kitamasa(const VR &C) :cs(C),rs(1),m(C.size()){ rs[0].assign(2*m+1,R::add_identity()); rs[0][1]=R::mul_identity(); } VR merge(const VR &xs,const VR &ys){ VR zs(2*m+1,R::add_identity()); for(int i=1;i<=m;i++) for(int j=1;j<=m;j++) zs[i+j]=zs[i+j]+(xs[i]*ys[j]); for(int i=2*m;i>m;zs[i--]=R::add_identity()) for(int j=1;j<=m;j++) zs[i-j]=zs[i-j]+(cs[m-j]*zs[i]); return zs; } R calc(const VR &A,long long n){ assert((int)A.size()==m); VR res(rs[0]); for(int i=0;n;i++,n>>=1){ if(i>=(int)rs.size()) rs.emplace_back(merge(rs[i-1],rs[i-1])); if(~n&1) continue; res=merge(res,rs[i]); } R ans=R::add_identity(); for(int i=1;i<=m;i++) ans=ans+(res[i]*A[i-1]); return ans; } }; //END CUT HERE #ifndef call_from_test #define call_from_test #include "../mod/mint.cpp" #undef call_from_test //INSERT ABOVE HERE signed TDPC_F(){ using ll = long long; int n; ll k; scanf("%d %lld",&n,&k); using M = Mint; vector as(n,1),cs(n,1); Kitamasa fib(cs); printf("%d\n",fib.calc(as,--k).v); return 0; } /* verified on 2019/12/10 https://atcoder.jp/contests/tdpc/tasks/tdpc_fibonacci */ signed ABC009_D(){ int k,m; scanf("%d %d",&k,&m); struct M{ uint32_t v; M(){*this=add_identity();} M(uint32_t v):v(v){} M operator+(const M &a)const{return M(v^a.v);} M operator*(const M &a)const{return M(v&a.v);} static M add_identity(){return M(0);} static M mul_identity(){return M(~0u);} }; vector as(k),cs(k); for(int i=0;i v(cs); printf("%u\n",v.calc(as,--m).v); return 0; } /* verified on 2019/12/10 https://atcoder.jp/contests/abc009/tasks/abc009_4 */ signed main(){ //TDPC_F(); //ABC009_D(); return 0; } #endif #undef call_from_test signed main(){ using M = Mint; const int n = 4; vector as(n,0),cs(n,1); as[n-1]=1; Kitamasa fib(cs); int q; cin>>q; for(int i=0;i>p; cout<