#include using namespace std; const int mod=998244353; class mint{ private: int x; long long ksm(long long w,int y) { if(!y)return 1; if(y%2==0)return ksm(w*w%mod,y>>1); return w*ksm(w*w%mod,y>>1)%mod; } public: void read() { int fl=0;x=0; char ch=getchar(); while(ch<'0'||ch>'9') { if(ch=='-')fl=1; ch=getchar(); } while('0'<=ch&&ch<='9') { x=(10ll*x+ch-'0')%mod; ch=getchar(); } if(fl&&x)ch=mod-x; } void print(){cout<=mod)a.x-=mod;return a;} mint operator~(void){mint a;a.x=x;if(a.x)a.x=mod-a.x;return a;} mint operator-(mint y)const{mint a;a=x;return a+(~y);} mint operator*(mint y)const{mint a;a.x=1ll*x*y.x%mod;return a;} mint inv(){mint a;a.x=ksm(x,mod-2);return a;} mint operator/(mint y)const{mint a;a=x;return a*(y.inv());} int val(){return x;} mint &operator++(void){x++;if(x==mod)x=0;return *this;} mint &operator--(void){x--;if(x==-1)x=mod-1;return *this;} mint &operator+=(mint y){x+=y.x;if(x>=mod)x-=mod;return *this;} mint &operator-=(mint y){x-=y.x;if(x<0)x+=mod;return *this;} mint &operator*=(mint y){x=1ll*x*y.x%mod;return *this;} mint &operator/=(mint y){x=1ll*x*(y.inv()).x%mod;return *this;} friend istream &operator >>(istream&,mint&); friend ostream &operator <<(ostream&,mint); }; istream &operator >>(istream& input,mint& x){ x.read(); return input; } ostream &operator <<(ostream& output,mint x){ x.print(); return output; } mint m0[505][505],m1[505][505]; mint a[505][505]; mint pl[505][505]; mint ans[505]; //求a的特征多项式 void char_poly(int n) { //hessenberg for(int i=1;i<=n;i++) { int p=-1; for(int j=i+1;j<=n;j++) { if(a[j][i].val()) { p=j; break; } } if(p==-1)continue; swap(a[i+1],a[p]); for(int j=1;j<=n;j++)swap(a[j][i+1],a[j][p]); for(int j=i+2;j<=n;j++) { mint v=a[j][i]/a[i+1][i]; for(int k=1;k<=n;k++)a[j][k]-=a[i+1][k]*v; for(int k=1;k<=n;k++)a[k][i+1]+=a[k][j]*v; } } pl[0][0]=1; for(int i=1;i<=n;i++) { for(int j=1;j=1;i--) { for(int j=i-1;j>=1;j--) { mint v=g[j][i]; for(int k=1;k<=n;k++) { g[j][k]-=g[i][k]*v; h[j][k]-=h[i][k]*v; } } } } mint r[505][505]; //r=g*h void product(int n) { for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { r[i][j]=0; } } for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { for(int k=1;k<=n;k++) { r[i][j]+=g[i][k]*h[k][j]; } } } } //det of h mint det(int n) { mint ret; ret=1; for(int i=1;i<=n;i++) { int p=-1; for(int j=1;j<=n;j++) { if(h[i][j].val()) { p=j; break; } } if(p==-1)continue; if(p!=i) { swap(h[i],h[p]); ret=~ret; } for(int j=i+1;j<=n;j++) { mint v=h[j][i]/h[i][i]; for(int k=1;k<=n;k++) { h[j][k]-=h[i][k]*v; } } } for(int i=1;i<=n;i++)ret*=h[i][i]; return ret; } mt19937 rnd(time(0)); mint sum[505]; mint cho[505][505]; mint mi[505]; int main() { // freopen("ex_5.in","r",stdin); int n; std::ios::sync_with_stdio(false); std::cin.tie(nullptr); cin>>n; int x; for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { cin>>x; m0[i][j]=x; } } for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { cin>>x; m1[i][j]=x; } } mint w; w=rnd()%(mod-1)+1; // w=1; for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { g[i][j]=m0[i][j]-w*m1[i][j]; } } matriv(n); for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { g[i][j]=m1[i][j]; } } product(n); for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { a[i][j]=~r[i][j]; } } char_poly(n); mint _v=det(n).inv(); for(int i=0;i<=n;i++)sum[i]=ans[n-i]*_v; // for(int i=0;i<=n;i++)cout<