#include<bits/stdc++.h>
using namespace std;
const int N=5005,mod=998244353;
int n,c[N][N],cnt[2][N],pre[2][N],suf[2][N],ans;
char a[N][N];
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)cin>>a[i][j];
	}
	c[0][0]=1;
	for(int i=1;i<=n;i++)
	{
		c[i][0]=1;
		for(int j=1;j<=i;j++)c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
	}
	for(int i=1;i<=n;i++)
	{
		int s=0;
		for(int j=1;j<=n;j++)s+=a[i][j]=='1';
		cnt[0][i]=s;
	}
	for(int j=1;j<=n;j++)
	{
		int s=0;
		for(int i=1;i<=n;i++)s+=a[i][j]=='1';
		cnt[1][j]=s;
	}
	for(int i=0;i<2;i++)sort(cnt[i]+1,cnt[i]+n+1),reverse(cnt[i]+1,cnt[i]+n+1);
	for(int i=0;i<2;i++)
	{
		for(int j=1;j<=n;j++)pre[i][j]=pre[i][j-1]+cnt[i][j];
		for(int j=n;j;j--)suf[i][j]=suf[i][j+1]+n-cnt[i][j];
	}
	for(int i=0;i<=n;i++)
	{
		int l1=i;
		while(l1<n&&cnt[0][l1]==cnt[0][l1+1])l1++;
		for(int j=0;j<=n;j++)
		{
			int l2=j;
			while(l2<n&&cnt[1][l2]==cnt[1][l2+1])l2++;
			for(int k=i;k<=l1;k++)
			{
				for(int k1=j;k1<=l2;k1++)if(pre[0][k]+pre[1][k1]+suf[0][k+1]+suf[1][k1+1]-k*k1-(n-k)*(n-k1)==n*n)
				{
					ans=(ans+1ll*c[l1-i+1][k-i+1]*c[l2-j+1][k1-j+1])%mod; 
				}
			}
			j=l2;
		}
		i=l1;
	}
	cout<<ans;
	return 0;
}