#include<bits/stdc++.h>
#define ll long long
using namespace std;

const ll mod=998244353;

int n,a[5005],b[5005],d[5005],e[5005],f,g,h;
ll fac[5005],el[5005],base[5005],ans;
char c[5005][5005];
vector<int> adj[5005];

inline ll calc(int x,int y)
{
    if(y<0||y>x)
        return 0;
    return fac[x]*el[y]%mod*el[x-y]%mod;
}

int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0),cout.tie(0);

    cin>>n;

    base[0]=fac[0]=el[0]=el[1]=1;
    for(int i=2;i<=n;i++)
        el[i]=(mod-mod/i)*el[mod%i]%mod;
    for(int i=1;i<=n;i++)
        base[i]=(base[i-1]<<1)%mod,fac[i]=fac[i-1]*i%mod,el[i]=el[i-1]*el[i]%mod;

    for(int i=1;i<=n;i++)
    {
        cin>>(c[i]+1);

        for(int j=1;j<=n;j++)
        {
            if(c[i][j]=='1')
                b[i]++;
            else
                e[j]++;
        }

        adj[b[i]].push_back(i);
    }

    for(int i=0;i<=n;i++)
    {
        for(int j:adj[i])
        {
            for(int k=1;k<=n;k++)
            {
                if(c[j][k]=='1')
                    d[k]++;
            }
        }

        f=n,g=i,h=1;

        for(int j=1;j<=n;j++)
        {
            if(d[j]&&e[j])
                h=0;

            if(d[j])
                f--,g--;
            if(e[j])
                f--;
        }

        ans=(ans+h*base[adj[i].size()]*calc(f,g))%mod;

        for(int j:adj[i])
        {
            for(int k=1;k<=n;k++)
            {
                if(c[j][k]=='0')
                    e[k]--;
            }
        }
    }

    cout<<ans<<'\n';

    return 0;
}