#include<bits/stdc++.h>
using namespace std;
#define ll long long 
const int mod=998244353;
const int maxn=5005;
char s[maxn][maxn];
int n;
int cnta[maxn],cntb[maxn];
int c[maxn][maxn];
inline ll qpow(ll a,ll b){
	ll ret=1;
	for(;b;b>>=1,a=a*a%mod)
		if(b&1) ret=ret*a%mod;
	return ret;
}
ll F[maxn],S[maxn];
int main(){
	scanf("%d",&n);
	c[0][0]=c[1][0]=c[1][1]=1;
	for(int i=2;i<=n;i++){
		c[i][0]=1;
		for(int j=1;j<=n;j++) c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
	}
	for(int i=1;i<=n;i++) scanf("%s",s[i]+1);
	for(int i=1;i<=n;i++){
		int c1=0,c2=0;
		for(int j=1;j<=n;j++){
			if(s[i][j]=='1') c1++;
			if(s[j][i]=='1') c2++;
		}
		cnta[c1]++,cntb[c2]++;
	}
	ll sum=0,sum2=0,f=1,ans=0;
	for(int i=0;i<=n;i++) sum+=(ll)(n-i)*cnta[i];
	ll f2=1;
	for(int j=0;j<=n;j++) sum2+=(ll)(n-j)*cntb[j];
	for(int y=0,j=n,cb=0;y<=n;y++){
		F[y]=f2,S[y]=sum2;
		if(y==n) continue;
		if(cntb[j] && cntb[j]>cb){
			f2=f2*qpow(c[cntb[j]][cb],mod-2)%mod;
			sum2-=(n-j);
			sum2+=j;
			cb++;
			f2=f2*(ll)c[cntb[j]][cb]%mod;
		}
		else{
			j--,cb=1;
			while(cntb[j]==0) j--;
			sum2-=(n-j);
			sum2+=j;
			f2=f2*(ll)c[cntb[j]][cb]%mod;
		}
	}
	for(int x=0,i=n,ca=0;x<=n;x++){
		for(int y=0;y<=n;y++){
			if(sum+S[y]==(ll)n*n+(ll)x*y+(ll)(n-x)*(n-y)){
				ans=(ans+f*F[y]%mod)%mod;
			}
		}
		if(x==n) continue;
		if(cnta[i] && cnta[i]>ca){
			f=f*qpow(c[cnta[i]][ca],mod-2)%mod;
			sum-=(n-i);
			sum+=i;
			ca++;
			f=f*(ll)c[cnta[i]][ca]%mod;
		}
		else{
			i--,ca=1;
			while(cnta[i]==0) i--;
			sum-=(n-i);
			sum+=i;
			f=f*(ll)c[cnta[i]][ca]%mod;
		}
	}
	printf("%lld",ans);
	return 0;
}