#include<map>
#include<set>
#include<cmath>
#include<ctime>
#include<queue>
#include<stack>
#include<cstdio>
#include<vector>
#include<string>
#include<bitset>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXN=5e3+10;
const int N=5e3;
const int INF=0x3f3f3f3f;
const long long LINF=0x3f3f3f3f3f3f3f3f;
const int mod=998244353;
int n,ans=0;
int cx[MAXN],cy[MAXN];
int prex[MAXN],sufx[MAXN];
int prey[MAXN],sufy[MAXN];
int cntx[MAXN],cnty[MAXN];
int firx[MAXN],firy[MAXN];
int comb[MAXN][MAXN];
signed main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            int x;
            scanf("%1d",&x);
            cx[i]+=x;
            cy[j]+=x;
        }
    }
    sort(cx+1,cx+1+n,greater<int>());
    sort(cy+1,cy+1+n,greater<int>());
    for(int i=1;i<=n;i++)
    {
        prex[i]=prex[i-1]+cx[i];
        prey[i]=prey[i-1]+cy[i];
    }
    for(int i=n;i;i--)
    {
        sufx[i]=sufx[i+1]+(n-cx[i]);
        sufy[i]=sufy[i+1]+(n-cy[i]);
        cntx[cx[i]]++;
        cnty[cy[i]]++;
        firx[cx[i]]=i;
        firy[cy[i]]=i;
    }
    comb[0][0]=1;
    for(int i=1;i<=n;i++)
    {
        comb[i][0]=1;
        for(int j=1;j<=i;j++)
        {
            comb[i][j]=comb[i-1][j-1]+comb[i-1][j];
            (comb[i][j]>=mod)&&(comb[i][j]-=mod);
        }
    }
    for(int i=0;i<=n;i++)
    {
        for(int j=0;j<=n;j++)
        {
            #define fx() (i?comb[cntx[cx[i]]][i-firx[cx[i]]+1]:1)
            #define fy() (j?comb[cnty[cy[j]]][j-firy[cy[j]]+1]:1)
            if((prex[i]+sufx[i+1])+(prey[j]+sufy[j+1])-i*j-(n-i)*(n-j)==n*n){
                ans=(ans+1ll*fx()*fy())%mod;
            }
        }
    }
    printf("%d\n",ans);
    return 0;
}