/*

                                _/                                      _/
       _/_/        _/_/_/      _/_/_/        _/_/_/      _/  _/_/               _/_/_/
    _/    _/    _/            _/    _/    _/    _/      _/_/          _/       _/    _/
   _/    _/    _/            _/    _/    _/    _/      _/            _/       _/    _/
    _/_/        _/_/_/      _/    _/      _/_/_/      _/            _/       _/    _/

*/
#include<bits/stdc++.h>
#define int long long

using namespace std;
const int mod=998244353;

void solve(){
    int n;
    cin>>n;
    vector<string>s(n);
    for(auto &i:s) cin>>i,i=i+' ';
    vector f(n+1,vector(n+1,0ll));
    vector mi0(n,n),ma1(n,-1ll);
    f[0][0]=1;
    for(int i=0;i<n;++i){
        int max1=-1,min0=n;
        for(int j=0;j<=n;++j){
            int cl=max1+1,cr=min(j,min0);
            if(cl<=cr) (f[i+1][j]+=f[i][j]*(cr-cl+1)%mod)%=mod;
            if(s[i][j]=='1') max1=max(max1,j);
            if(s[i][j]=='0') min0=min(min0,j);
        }
        for(int j=0;j<n;++j){
            int cl=ma1[j]+1,cr=min(i,mi0[j]);
            if(cl<=cr) (f[i][j+1]+=f[i][j]*(cr-cl+1)%mod)%=mod;
            cr=min(cr,i-1);
            if(cl<=cr && min0>j) (f[i+1][j+1]+=f[i][j]*(cr-cl+1)%mod)%=mod;
            if(s[i][j]=='1') ma1[j]=max(ma1[j],i);
            if(s[i][j]=='0') mi0[j]=min(mi0[j],i);
        }
    }
    for(int i=0;i<n;++i){
        int cl=ma1[i]+1,cr=min(n,mi0[i]);
        if(cl<=cr) (f[n][i+1]+=f[n][i]*(cr-cl+1)%mod)%=mod;
    }
    cout<<f[n][n]<<"\n";
}

signed main(){
    ios::sync_with_stdio(false);
    cin.tie(0);solve();
}