#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=5010,mod=998244353;
int n,c[N][N],a[N],b[N];
int ed1[N],ed2[N];
int jc[N],invjc[N];
int ksm(int x,int y){int ans=1;while(y){if(y&1)ans=ans*x%mod;x=x*x%mod;y>>=1;}return ans;}
int inv(int x){return ksm(x,mod-2);}
void init(){
	jc[0]=1;for(int i=1;i<=n;i++)jc[i]=jc[i-1]*i%mod;
	invjc[n]=inv(jc[n]);
	for(int i=n-1;i>=0;i--)invjc[i]=invjc[i+1]*(i+1)%mod;
}
int C(int n,int m){
	if(n<m)return 0;
	return jc[n]*invjc[m]%mod*invjc[n-m]%mod;
}
char ch[N];
int sum1[N],sum2[N];
int ans=0;
signed main(){
	scanf("%lld",&n);init();
	ans=1;
	for(int i=1;i<=n;i++){
		scanf("%s",ch+1);
		for(int j=1;j<=n;j++){
			c[i][j]=ch[j]-'0';
			if(c[i][j])ans=0;
			a[i]+=c[i][j];
			b[j]+=c[i][j];
		}
	}
	sort(a+1,a+n+1);reverse(a+1,a+n+1);
	sort(b+1,b+n+1);reverse(b+1,b+n+1);
	for(int i=n;i>=1;i--)sum1[i]=sum1[i+1]+n-a[i];
	for(int i=n;i>=1;i--)sum2[i]=sum2[i+1]+n-b[i];
	int s1=0,s2=0,l1=1,l2=1;
	for(int i=1;i<=n;i++){
		if(a[i]!=a[l1]){
			ed1[l1]=i-1-l1+1;l1=i;
		}
	}ed1[l1]=n-l1+1;
	for(int i=1;i<=n;i++){
		if(b[i]!=b[l2]){
			ed2[l2]=i-1-l2+1;l2=i;
		}
	}ed2[l2]=n-l2+1;
	// for(int i=1;i<=n;i++)cout<<ed1[i]<<" ";cout<<endl;
	// for(int i=1;i<=n;i++)cout<<ed2[i]<<" ";cout<<endl;cout<<endl;
	l1=l2=1;
	for(int i=1;i<=n;i++){
		s2=0;s1+=a[i];if(a[i]!=a[l1])l1=i;
		l2=1;
		for(int j=1;j<=n;j++){
			s2+=b[j];if(b[j]!=b[l2])l2=j;
			int s=s1+s2+sum1[i+1]+sum2[j+1]-(i*j+(n-i)*(n-j));
			// cout<<i<<" "<<j<<"  "<<s1<<" "<<s2<<" "<<sum1[i+1]<<" "<<sum2[j+1]<<" "<<s<<endl;
			if(s==n*n){
				// cout<<i<<" "<<j<<" "<<C(ed1[l1],i-l1+1)*C(ed2[l2],j-l2+1)%mod<<" "<<l1<<" "<<l2<<endl;
				ans+=C(ed1[l1],i-l1+1)*C(ed2[l2],j-l2+1)%mod;ans%=mod;
			}
		}
	}
	s2=0;l2=1;
	for(int j=1;j<=n;j++){
		s2+=b[j];if(b[j]!=b[l2])l2=j;
		int s=s2+sum1[1]+sum2[j+1]-n*(n-j);
		if(s==n*n){
			ans+=C(ed2[l2],j-l2+1)%mod;ans%=mod;
		}
	}
	s1=0,l1=1;
	for(int i=1;i<=n;i++){
		s1+=a[i];if(a[i]!=a[l1])l1=i;
		int s=s1+sum1[i+1]+sum2[1]-(n-i)*n;
		if(s==n*n){
			ans+=C(ed1[l1],i-l1+1)%mod;ans%=mod;
		}
	}
	cout<<ans;
	return 0;
}