#define PROBLEM "https://judge.yosupo.jp/problem/system_of_linear_equations"
//#include"../../template/template.hpp"//#include"../../modint/montgomery-modint.hpp"
#include"../../modulo/gauss-elimination-fast.hpp"usingmint=LazyMontgomeryModInt<998244353>;usingnamespaceNyaan;voidNyaan::solve(){ini(N,M);V<V<mint>>A(N,V<mint>(M));in(A);V<mint>B(N);in(B);autov=LinearEquation<mint>(A,B);out(sz(v)-1);each(x,v)out(x);}
#line 1 "verify/verify-yosupo-math/yosupo-linear-equation.test.cpp"
#define PROBLEM "https://judge.yosupo.jp/problem/system_of_linear_equations"
//#line 2 "template/template.hpp"
usingnamespacestd;// intrinstic#include<immintrin.h>#include<algorithm>
#include<array>
#include<bitset>
#include<cassert>
#include<cctype>
#include<cfenv>
#include<cfloat>
#include<chrono>
#include<cinttypes>
#include<climits>
#include<cmath>
#include<complex>
#include<cstdarg>
#include<cstddef>
#include<cstdint>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<deque>
#include<fstream>
#include<functional>
#include<initializer_list>
#include<iomanip>
#include<ios>
#include<iostream>
#include<istream>
#include<iterator>
#include<limits>
#include<list>
#include<map>
#include<memory>
#include<new>
#include<numeric>
#include<ostream>
#include<queue>
#include<random>
#include<set>
#include<sstream>
#include<stack>
#include<streambuf>
#include<string>
#include<tuple>
#include<type_traits>
#include<typeinfo>
#include<unordered_map>
#include<unordered_set>
#include<utility>
#include<vector>// utility#line 1 "template/util.hpp"
namespaceNyaan{usingll=longlong;usingi64=longlong;usingu64=unsignedlonglong;usingi128=__int128_t;usingu128=__uint128_t;template<typenameT>usingV=vector<T>;template<typenameT>usingVV=vector<vector<T>>;usingvi=vector<int>;usingvl=vector<longlong>;usingvd=V<double>;usingvs=V<string>;usingvvi=vector<vector<int>>;usingvvl=vector<vector<longlong>>;template<typenameT>usingminpq=priority_queue<T,vector<T>,greater<T>>;template<typenameT,typenameU>structP:pair<T,U>{template<typename...Args>P(Args...args):pair<T,U>(args...){}usingpair<T,U>::first;usingpair<T,U>::second;P&operator+=(constP&r){first+=r.first;second+=r.second;return*this;}P&operator-=(constP&r){first-=r.first;second-=r.second;return*this;}P&operator*=(constP&r){first*=r.first;second*=r.second;return*this;}template<typenameS>P&operator*=(constS&r){first*=r,second*=r;return*this;}Poperator+(constP&r)const{returnP(*this)+=r;}Poperator-(constP&r)const{returnP(*this)-=r;}Poperator*(constP&r)const{returnP(*this)*=r;}template<typenameS>Poperator*(constS&r)const{returnP(*this)*=r;}Poperator-()const{returnP{-first,-second};}};usingpl=P<ll,ll>;usingpi=P<int,int>;usingvp=V<pl>;constexprintinf=1001001001;constexprlonglonginfLL=4004004004004004004LL;template<typenameT>intsz(constT&t){returnt.size();}template<typenameT,typenameU>inlineboolamin(T&x,Uy){return(y<x)?(x=y,true):false;}template<typenameT,typenameU>inlineboolamax(T&x,Uy){return(x<y)?(x=y,true):false;}template<typenameT>inlineTMax(constvector<T>&v){return*max_element(begin(v),end(v));}template<typenameT>inlineTMin(constvector<T>&v){return*min_element(begin(v),end(v));}template<typenameT>inlinelonglongSum(constvector<T>&v){returnaccumulate(begin(v),end(v),0LL);}template<typenameT>intlb(constvector<T>&v,constT&a){returnlower_bound(begin(v),end(v),a)-begin(v);}template<typenameT>intub(constvector<T>&v,constT&a){returnupper_bound(begin(v),end(v),a)-begin(v);}constexprlonglongTEN(intn){longlongret=1,x=10;for(;n;x*=x,n>>=1)ret*=(n&1?x:1);returnret;}template<typenameT,typenameU>pair<T,U>mkp(constT&t,constU&u){returnmake_pair(t,u);}template<typenameT>vector<T>mkrui(constvector<T>&v,boolrev=false){vector<T>ret(v.size()+1);if(rev){for(inti=int(v.size())-1;i>=0;i--)ret[i]=v[i]+ret[i+1];}else{for(inti=0;i<int(v.size());i++)ret[i+1]=ret[i]+v[i];}returnret;};template<typenameT>vector<T>mkuni(constvector<T>&v){vector<T>ret(v);sort(ret.begin(),ret.end());ret.erase(unique(ret.begin(),ret.end()),ret.end());returnret;}template<typenameF>vector<int>mkord(intN,Ff){vector<int>ord(N);iota(begin(ord),end(ord),0);sort(begin(ord),end(ord),f);returnord;}template<typenameT>vector<int>mkinv(vector<T>&v){intmax_val=*max_element(begin(v),end(v));vector<int>inv(max_val+1,-1);for(inti=0;i<(int)v.size();i++)inv[v[i]]=i;returninv;}vector<int>mkiota(intn){vector<int>ret(n);iota(begin(ret),end(ret),0);returnret;}template<typenameT>Tmkrev(constT&v){Tw{v};reverse(begin(w),end(w));returnw;}template<typenameT>boolnxp(T&v){returnnext_permutation(begin(v),end(v));}// 返り値の型は入力の T に依存// i 要素目 : [0, a[i])template<typenameT>vector<vector<T>>product(constvector<T>&a){vector<vector<T>>ret;vector<T>v;autodfs=[&](autorc,inti)->void{if(i==(int)a.size()){ret.push_back(v);return;}for(intj=0;j<a[i];j++)v.push_back(j),rc(rc,i+1),v.pop_back();};dfs(dfs,0);returnret;}// F : function(void(T&)), mod を取る操作// T : 整数型のときはオーバーフローに注意するtemplate<typenameT>TPower(Ta,longlongn,constT&I,constfunction<void(T&)>&f){Tres=I;for(;n;f(a=a*a),n>>=1){if(n&1)f(res=res*a);}returnres;}// T : 整数型のときはオーバーフローに注意するtemplate<typenameT>TPower(Ta,longlongn,constT&I=T{1}){returnPower(a,n,I,function<void(T&)>{[](T&)->void{}});}template<typenameT>TRev(constT&v){Tres=v;reverse(begin(res),end(res));returnres;}template<typenameT>vector<T>Transpose(constvector<T>&v){usingU=typenameT::value_type;if(v.empty())return{};intH=v.size(),W=v[0].size();vectorres(W,T(H,U{}));for(inti=0;i<H;i++){for(intj=0;j<W;j++){res[j][i]=v[i][j];}}returnres;}template<typenameT>vector<T>Rotate(constvector<T>&v,intclockwise=true){usingU=typenameT::value_type;intH=v.size(),W=v[0].size();vectorres(W,T(H,U{}));for(inti=0;i<H;i++){for(intj=0;j<W;j++){if(clockwise){res[W-1-j][i]=v[i][j];}else{res[j][H-1-i]=v[i][j];}}}returnres;}}// namespace Nyaan#line 58 "template/template.hpp"
// bit operation#line 1 "template/bitop.hpp"
namespaceNyaan{__attribute__((target("popcnt")))inlineintpopcnt(constu64&a){return__builtin_popcountll(a);}inlineintlsb(constu64&a){returna?__builtin_ctzll(a):64;}inlineintctz(constu64&a){returna?__builtin_ctzll(a):64;}inlineintmsb(constu64&a){returna?63-__builtin_clzll(a):-1;}template<typenameT>inlineintgbit(constT&a,inti){return(a>>i)&1;}template<typenameT>inlinevoidsbit(T&a,inti,boolb){if(gbit(a,i)!=b)a^=T(1)<<i;}constexprlonglongPW(intn){return1LL<<n;}constexprlonglongMSK(intn){return(1LL<<n)-1;}}// namespace Nyaan#line 61 "template/template.hpp"
// inout#line 1 "template/inout.hpp"
namespaceNyaan{template<typenameT,typenameU>ostream&operator<<(ostream&os,constpair<T,U>&p){os<<p.first<<" "<<p.second;returnos;}template<typenameT,typenameU>istream&operator>>(istream&is,pair<T,U>&p){is>>p.first>>p.second;returnis;}template<typenameT>ostream&operator<<(ostream&os,constvector<T>&v){ints=(int)v.size();for(inti=0;i<s;i++)os<<(i?" ":"")<<v[i];returnos;}template<typenameT>istream&operator>>(istream&is,vector<T>&v){for(auto&x:v)is>>x;returnis;}istream&operator>>(istream&is,__int128_t&x){stringS;is>>S;x=0;intflag=0;for(auto&c:S){if(c=='-'){flag=true;continue;}x*=10;x+=c-'0';}if(flag)x=-x;returnis;}istream&operator>>(istream&is,__uint128_t&x){stringS;is>>S;x=0;for(auto&c:S){x*=10;x+=c-'0';}returnis;}ostream&operator<<(ostream&os,__int128_tx){if(x==0)returnos<<0;if(x<0)os<<'-',x=-x;stringS;while(x)S.push_back('0'+x%10),x/=10;reverse(begin(S),end(S));returnos<<S;}ostream&operator<<(ostream&os,__uint128_tx){if(x==0)returnos<<0;stringS;while(x)S.push_back('0'+x%10),x/=10;reverse(begin(S),end(S));returnos<<S;}voidin(){}template<typenameT,class...U>voidin(T&t,U&...u){cin>>t;in(u...);}voidout(){cout<<"\n";}template<typenameT,class...U,charsep=' '>voidout(constT&t,constU&...u){cout<<t;if(sizeof...(u))cout<<sep;out(u...);}structIoSetupNya{IoSetupNya(){cin.tie(nullptr);ios::sync_with_stdio(false);cout<<fixed<<setprecision(15);cerr<<fixed<<setprecision(7);}}iosetupnya;}// namespace Nyaan#line 64 "template/template.hpp"
// debug#line 1 "template/debug.hpp"
namespaceDebugImpl{template<typenameU,typename=void>structis_specialize:false_type{};template<typenameU>structis_specialize<U,typenameconditional<false,typenameU::iterator,void>::type>:true_type{};template<typenameU>structis_specialize<U,typenameconditional<false,decltype(U::first),void>::type>:true_type{};template<typenameU>structis_specialize<U,enable_if_t<is_integral<U>::value,void>>:true_type{};voiddump(constchar&t){cerr<<t;}voiddump(conststring&t){cerr<<t;}voiddump(constbool&t){cerr<<(t?"true":"false");}voiddump(__int128_tt){if(t==0)cerr<<0;if(t<0)cerr<<'-',t=-t;stringS;while(t)S.push_back('0'+t%10),t/=10;reverse(begin(S),end(S));cerr<<S;}voiddump(__uint128_tt){if(t==0)cerr<<0;stringS;while(t)S.push_back('0'+t%10),t/=10;reverse(begin(S),end(S));cerr<<S;}template<typenameU,enable_if_t<!is_specialize<U>::value,nullptr_t>=nullptr>voiddump(constU&t){cerr<<t;}template<typenameT>voiddump(constT&t,enable_if_t<is_integral<T>::value>*=nullptr){stringres;if(t==Nyaan::inf)res="inf";ifconstexpr(is_signed<T>::value){if(t==-Nyaan::inf)res="-inf";}ifconstexpr(sizeof(T)==8){if(t==Nyaan::infLL)res="inf";ifconstexpr(is_signed<T>::value){if(t==-Nyaan::infLL)res="-inf";}}if(res.empty())res=to_string(t);cerr<<res;}template<typenameT,typenameU>voiddump(constpair<T,U>&);template<typenameT>voiddump(constpair<T*,int>&);template<typenameT>voiddump(constT&t,enable_if_t<!is_void<typenameT::iterator>::value>*=nullptr){cerr<<"[ ";for(autoit=t.begin();it!=t.end();){dump(*it);cerr<<(++it==t.end()?"":", ");}cerr<<" ]";}template<typenameT,typenameU>voiddump(constpair<T,U>&t){cerr<<"( ";dump(t.first);cerr<<", ";dump(t.second);cerr<<" )";}template<typenameT>voiddump(constpair<T*,int>&t){cerr<<"[ ";for(inti=0;i<t.second;i++){dump(t.first[i]);cerr<<(i==t.second-1?"":", ");}cerr<<" ]";}voidtrace(){cerr<<endl;}template<typenameHead,typename...Tail>voidtrace(Head&&head,Tail&&...tail){cerr<<" ";dump(head);if(sizeof...(tail)!=0)cerr<<",";trace(std::forward<Tail>(tail)...);}}// namespace DebugImpl#ifdef NyaanDebug
#define trc(...) \
do { \
cerr << "## " << #__VA_ARGS__ << " = "; \
DebugImpl::trace(__VA_ARGS__); \
} while (0)
#else
#define trc(...) (void(0))
#endif
#ifdef NyaanLocal
#define trc2(...) \
do { \
cerr << "## " << #__VA_ARGS__ << " = "; \
DebugImpl::trace(__VA_ARGS__); \
} while (0)
#else
#define trc2(...) (void(0))
#endif
#line 67 "template/template.hpp"
// macro#line 1 "template/macro.hpp"
#define each(x, v) for (auto&& x : v)
#define each2(x, y, v) for (auto&& [x, y] : v)
#define all(v) (v).begin(), (v).end()
#define rep(i, N) for (long long i = 0; i < (long long)(N); i++)
#define repr(i, N) for (long long i = (long long)(N)-1; i >= 0; i--)
#define rep1(i, N) for (long long i = 1; i <= (long long)(N); i++)
#define repr1(i, N) for (long long i = (N); (long long)(i) > 0; i--)
#define reg(i, a, b) for (long long i = (a); i < (b); i++)
#define regr(i, a, b) for (long long i = (b)-1; i >= (a); i--)
#define fi first
#define se second
#define ini(...) \
int __VA_ARGS__; \
in(__VA_ARGS__)
#define inl(...) \
long long __VA_ARGS__; \
in(__VA_ARGS__)
#define ins(...) \
string __VA_ARGS__; \
in(__VA_ARGS__)
#define in2(s, t) \
for (int i = 0; i < (int)s.size(); i++) { \
in(s[i], t[i]); \
}
#define in3(s, t, u) \
for (int i = 0; i < (int)s.size(); i++) { \
in(s[i], t[i], u[i]); \
}
#define in4(s, t, u, v) \
for (int i = 0; i < (int)s.size(); i++) { \
in(s[i], t[i], u[i], v[i]); \
}
#define die(...) \
do { \
Nyaan::out(__VA_ARGS__); \
return; \
} while (0)
#line 70 "template/template.hpp"
namespaceNyaan{voidsolve();}intmain(){Nyaan::solve();}#line 4 "verify/verify-yosupo-math/yosupo-linear-equation.test.cpp"
//#line 2 "modint/montgomery-modint.hpp"
template<uint32_tmod>structLazyMontgomeryModInt{usingmint=LazyMontgomeryModInt;usingi32=int32_t;usingu32=uint32_t;usingu64=uint64_t;staticconstexpru32get_r(){u32ret=mod;for(i32i=0;i<4;++i)ret*=2-mod*ret;returnret;}staticconstexpru32r=get_r();staticconstexpru32n2=-u64(mod)%mod;static_assert(mod<(1<<30),"invalid, mod >= 2 ^ 30");static_assert((mod&1)==1,"invalid, mod % 2 == 0");static_assert(r*mod==1,"this code has bugs.");u32a;constexprLazyMontgomeryModInt():a(0){}constexprLazyMontgomeryModInt(constint64_t&b):a(reduce(u64(b%mod+mod)*n2)){};staticconstexpru32reduce(constu64&b){return(b+u64(u32(b)*u32(-r))*mod)>>32;}constexprmint&operator+=(constmint&b){if(i32(a+=b.a-2*mod)<0)a+=2*mod;return*this;}constexprmint&operator-=(constmint&b){if(i32(a-=b.a)<0)a+=2*mod;return*this;}constexprmint&operator*=(constmint&b){a=reduce(u64(a)*b.a);return*this;}constexprmint&operator/=(constmint&b){*this*=b.inverse();return*this;}constexprmintoperator+(constmint&b)const{returnmint(*this)+=b;}constexprmintoperator-(constmint&b)const{returnmint(*this)-=b;}constexprmintoperator*(constmint&b)const{returnmint(*this)*=b;}constexprmintoperator/(constmint&b)const{returnmint(*this)/=b;}constexprbooloperator==(constmint&b)const{return(a>=mod?a-mod:a)==(b.a>=mod?b.a-mod:b.a);}constexprbooloperator!=(constmint&b)const{return(a>=mod?a-mod:a)!=(b.a>=mod?b.a-mod:b.a);}constexprmintoperator-()const{returnmint()-mint(*this);}constexprmintoperator+()const{returnmint(*this);}constexprmintpow(u64n)const{mintret(1),mul(*this);while(n>0){if(n&1)ret*=mul;mul*=mul;n>>=1;}returnret;}constexprmintinverse()const{intx=get(),y=mod,u=1,v=0,t=0,tmp=0;while(y>0){t=x/y;x-=t*y,u-=t*v;tmp=x,x=y,y=tmp;tmp=u,u=v,v=tmp;}returnmint{u};}friendostream&operator<<(ostream&os,constmint&b){returnos<<b.get();}friendistream&operator>>(istream&is,mint&b){int64_tt;is>>t;b=LazyMontgomeryModInt<mod>(t);return(is);}constexpru32get()const{u32ret=reduce(a);returnret>=mod?ret-mod:ret;}staticconstexpru32get_mod(){returnmod;}};#line 2 "modulo/gauss-elimination-fast.hpp"
#line 2 "modint/simd-montgomery.hpp"
#line 4 "modint/simd-montgomery.hpp"
__attribute__((target("sse4.2")))inline__m128imy128_mullo_epu32(const__m128i&a,const__m128i&b){return_mm_mullo_epi32(a,b);}__attribute__((target("sse4.2")))inline__m128imy128_mulhi_epu32(const__m128i&a,const__m128i&b){__m128ia13=_mm_shuffle_epi32(a,0xF5);__m128ib13=_mm_shuffle_epi32(b,0xF5);__m128iprod02=_mm_mul_epu32(a,b);__m128iprod13=_mm_mul_epu32(a13,b13);__m128iprod=_mm_unpackhi_epi64(_mm_unpacklo_epi32(prod02,prod13),_mm_unpackhi_epi32(prod02,prod13));returnprod;}__attribute__((target("sse4.2")))inline__m128imontgomery_mul_128(const__m128i&a,const__m128i&b,const__m128i&r,const__m128i&m1){return_mm_sub_epi32(_mm_add_epi32(my128_mulhi_epu32(a,b),m1),my128_mulhi_epu32(my128_mullo_epu32(my128_mullo_epu32(a,b),r),m1));}__attribute__((target("sse4.2")))inline__m128imontgomery_add_128(const__m128i&a,const__m128i&b,const__m128i&m2,const__m128i&m0){__m128iret=_mm_sub_epi32(_mm_add_epi32(a,b),m2);return_mm_add_epi32(_mm_and_si128(_mm_cmpgt_epi32(m0,ret),m2),ret);}__attribute__((target("sse4.2")))inline__m128imontgomery_sub_128(const__m128i&a,const__m128i&b,const__m128i&m2,const__m128i&m0){__m128iret=_mm_sub_epi32(a,b);return_mm_add_epi32(_mm_and_si128(_mm_cmpgt_epi32(m0,ret),m2),ret);}__attribute__((target("avx2")))inline__m256imy256_mullo_epu32(const__m256i&a,const__m256i&b){return_mm256_mullo_epi32(a,b);}__attribute__((target("avx2")))inline__m256imy256_mulhi_epu32(const__m256i&a,const__m256i&b){__m256ia13=_mm256_shuffle_epi32(a,0xF5);__m256ib13=_mm256_shuffle_epi32(b,0xF5);__m256iprod02=_mm256_mul_epu32(a,b);__m256iprod13=_mm256_mul_epu32(a13,b13);__m256iprod=_mm256_unpackhi_epi64(_mm256_unpacklo_epi32(prod02,prod13),_mm256_unpackhi_epi32(prod02,prod13));returnprod;}__attribute__((target("avx2")))inline__m256imontgomery_mul_256(const__m256i&a,const__m256i&b,const__m256i&r,const__m256i&m1){return_mm256_sub_epi32(_mm256_add_epi32(my256_mulhi_epu32(a,b),m1),my256_mulhi_epu32(my256_mullo_epu32(my256_mullo_epu32(a,b),r),m1));}__attribute__((target("avx2")))inline__m256imontgomery_add_256(const__m256i&a,const__m256i&b,const__m256i&m2,const__m256i&m0){__m256iret=_mm256_sub_epi32(_mm256_add_epi32(a,b),m2);return_mm256_add_epi32(_mm256_and_si256(_mm256_cmpgt_epi32(m0,ret),m2),ret);}__attribute__((target("avx2")))inline__m256imontgomery_sub_256(const__m256i&a,const__m256i&b,const__m256i&m2,const__m256i&m0){__m256iret=_mm256_sub_epi32(a,b);return_mm256_add_epi32(_mm256_and_si256(_mm256_cmpgt_epi32(m0,ret),m2),ret);}#line 4 "modulo/gauss-elimination-fast.hpp"
namespaceGauss{constexprintMAT_SIZE=4096;uint32_ta_buf_[MAT_SIZE][MAT_SIZE]__attribute__((aligned(64)));// return value: (rank, (-1) ^ (number of swap time))template<typenamemint>__attribute__((target("avx2")))pair<int,mint>GaussianElimination(constvector<vector<mint>>&m,intLinearEquation=false){mint(&a)[MAT_SIZE][MAT_SIZE]=*reinterpret_cast<mint(*)[MAT_SIZE][MAT_SIZE]>(a_buf_);intH=m.size(),W=m[0].size(),rank=0;mintdet=1;for(inti=0;i<H;i++)for(intj=0;j<W;j++)a[i][j].a=m[i][j].a;__m256ir=_mm256_set1_epi32(mint::r);__m256im0=_mm256_set1_epi32(0);__m256im1=_mm256_set1_epi32(mint::get_mod());__m256im2=_mm256_set1_epi32(mint::get_mod()<<1);for(intj=0;j<(LinearEquation?(W-1):W);j++){// find basisif(rank==H)break;intidx=-1;for(inti=rank;i<H;i++){if(a[i][j].get()!=0){idx=i;break;}}if(idx==-1){det=0;continue;}// swapif(rank!=idx){det=-det;for(intl=j;l<W;l++)swap(a[rank][l],a[idx][l]);}det*=a[rank][j];// normalizeif(LinearEquation){if(a[rank][j].get()!=1){mintcoeff=a[rank][j].inverse();__m256iCOEFF=_mm256_set1_epi32(coeff.a);for(inti=j/8*8;i<W;i+=8){__m256iR=_mm256_load_si256((__m256i*)(a[rank]+i));__m256iRmulC=montgomery_mul_256(R,COEFF,r,m1);_mm256_store_si256((__m256i*)(a[rank]+i),RmulC);}}}// eliminationfor(intk=(LinearEquation?0:rank+1);k<H;k++){if(k==rank)continue;if(a[k][j].get()!=0){mintcoeff=a[k][j]/a[rank][j];__m256iCOEFF=_mm256_set1_epi32(coeff.a);for(inti=j/8*8;i<W;i+=8){__m256iR=_mm256_load_si256((__m256i*)(a[rank]+i));__m256iK=_mm256_load_si256((__m256i*)(a[k]+i));__m256iRmulC=montgomery_mul_256(R,COEFF,r,m1);__m256iKmnsR=montgomery_sub_256(K,RmulC,m2,m0);_mm256_store_si256((__m256i*)(a[k]+i),KmnsR);}}}rank++;}return{rank,det};}// calculate determinanttemplate<typenamemint>mintdeterminant(constvector<vector<mint>>&mat){returnGaussianElimination(mat).second;}// return V<V<mint>>// 0 column ... one of solutions// 1 ~ (W - rank) column ... bases// if not exist, return empty vectortemplate<typenamemint>vector<vector<mint>>LinearEquation(vector<vector<mint>>A,vector<mint>B){intH=A.size(),W=A[0].size();for(inti=0;i<H;i++)A[i].push_back(B[i]);autop=GaussianElimination(A,true);mint(&a)[MAT_SIZE][MAT_SIZE]=*reinterpret_cast<mint(*)[MAT_SIZE][MAT_SIZE]>(a_buf_);intrank=p.first;// check if solutions existfor(inti=rank;i<H;++i)if(a[i][W]!=0)returnvector<vector<mint>>{};vector<vector<mint>>res(1,vector<mint>(W));vector<int>pivot(W,-1);for(inti=0,j=0;i<rank;++i){while(a[i][j]==0)++j;res[0][j]=a[i][W],pivot[j]=i;}for(intj=0;j<W;++j){if(pivot[j]==-1){vector<mint>x(W);x[j]=1;for(intk=0;k<j;++k)if(pivot[k]!=-1)x[k]=-a[pivot[k]][j];res.push_back(x);}}returnres;}}// namespace GaussusingGauss::determinant;usingGauss::LinearEquation;#line 7 "verify/verify-yosupo-math/yosupo-linear-equation.test.cpp"
usingmint=LazyMontgomeryModInt<998244353>;usingnamespaceNyaan;voidNyaan::solve(){ini(N,M);V<V<mint>>A(N,V<mint>(M));in(A);V<mint>B(N);in(B);autov=LinearEquation<mint>(A,B);out(sz(v)-1);each(x,v)out(x);}