#include using namespace std; using ll=long long; using ull=unsigned long long; using P=pair; templateusing minque=priority_queue,greater>; templatebool chmax(T &a,const T &b){return (abool chmin(T &a,const T &b){return (a>b?(a=b,true):false);} templateistream &operator>>(istream &is,pair&p){is>>p.first>>p.second;return is;} templateistream &operator>>(istream &is,tuple&a){is>>std::get<0>(a)>>std::get<1>(a)>>std::get<2>(a);return is;} templateistream &operator>>(istream &is,array&a){for(auto&i:a)is>>i;return is;} templateistream &operator>>(istream &is,vector &a){for(auto &i:a)is>>i;return is;} templatevoid operator++(pair&a,int n){a.first++,a.second++;} templatevoid operator--(pair&a,int n){a.first--,a.second--;} templatevoid operator++(vector&a,int n){for(auto &i:a)i++;} templatevoid operator--(vector&a,int n){for(auto &i:a)i--;} #define overload3(_1,_2,_3,name,...) name #define rep1(i,n) for(int i=0;i<(int)(n);i++) #define rep2(i,l,r) for(int i=(int)(l);i<(int)(r);i++) #define rep(...) overload3(__VA_ARGS__,rep2,rep1)(__VA_ARGS__) #define reps(i,l,r) rep2(i,l,r) #define all(x) x.begin(),x.end() #define pcnt(x) __builtin_popcountll(x) #define fin(x) return cout<<(x)<<'\n',static_cast(0) #define yn(x) cout<<((x)?"Yes\n":"No\n") #define uniq(x) sort(all(x)),x.erase(unique(all(x)),x.end()) template inline int fkey(vector&z,T key){return lower_bound(z.begin(),z.end(),key)-z.begin();} ll myceil(ll a,ll b){return (a+b-1)/b;} template auto vec(const int (&d)[n],const T &init=T()){ if constexpr (id(d,init)); else return init; } #ifdef LOCAL #include #define SWITCH(a,b) (a) #else #define debug(...) static_cast(0) #define debugg(...) static_cast(0) #define SWITCH(a,b) (b) templateostream &operator<<(ostream &os,const pair&p){os<>testcase; for(int i=0;i template struct ConstantMatrix{ private: std::array,N>dat; public: ConstantMatrix(){ dat.fill([](){std::arrayres;res.fill(T());return res;}()); } static ConstantMatrix E(){ static_assert(N==M); ConstantMatrix res; for(int i=0;ifriend ConstantMatrix operator*(const ConstantMatrix&lhs,const ConstantMatrix&rhs){ ConstantMatrixres; for(int i=0;i &operator[](int i){return dat[i];} const std::array &operator[](int i)const{return dat[i];} friend bool operator==(const ConstantMatrix&lhs,const ConstantMatrix&rhs){ for(int i=0;i ConstantMatrix pow(U k)const{ static_assert(N==M); ConstantMatrix res=ConstantMatrix::E(),a(*this); while(k){ if(k&1)res*=a; a*=a; k>>=1; } return res; } std::vector>to_vector()const{ std::vector>res(N,std::vector(M)); for(int i=0;i>(std::istream&is,ConstantMatrix&rhs){ for(int i=0;i>rhs[i][j]; return is; } friend std::ostream &operator<<(std::ostream &os,const ConstantMatrix&rhs){ for(int i=0;iM,M2; cin>>M>>M2; auto ans=M*M2*M*M2; rep(i,2)rep(j,2)cout<