#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
constexpr int N=3e3+7,mod=998244353;

int n;
char s[N][N];
int f[N][N][2],h1[N][N],h2[N][N];
int w1[N][N],w2[N][N];

int main(){
	cin>>n;
	for(int i=1;i<=n;i++)scanf("%s",s[i]+1);
	for(int i=1;i<=n;i++){
		h1[i][0]=w1[i][0]=1;
		bool flg=true;
		for(int j=1;j<=n;j++){
			if(s[i][j]=='0')flg=false;
			if(s[i][j]!='1')h1[i][j]=h1[i][j-1];
			if(s[i][j]!='0'&&flg)h1[i][j]++,w1[i][j]=1;
		}
	}
	for(int i=1;i<=n;i++){
		h2[0][i]=w2[0][i]=1;
		bool flg=true;
		for(int j=1;j<=n;j++){
			if(s[j][i]=='0')flg=false;
			if(s[j][i]!='1')h2[j][i]=h2[j-1][i];
			if(s[j][i]!='0'&&flg)h2[j][i]++,w2[j][i]=1;
		}
	}
	f[n][n][0]=1;
	for(int i=n;i>=0;i--)
		for(int j=n;j>=0;j--){
			int v;
			if(f[i][j][0]){
				v=f[i][j][0];
				if(j)f[i][j-1][1]=(f[i][j-1][1]+(ll)v*h2[i][j])%mod;
				if(i)f[i-1][j][0]=(f[i-1][j][0]+(ll)v*h1[i][j])%mod;
				if(i&&j)f[i-1][j-1][0]=(f[i-1][j-1][0]+(ll)v*w2[i][j]*(h1[i][j-1]-w1[i][j-1]))%mod;
			}
			if(f[i][j][1]){
				v=f[i][j][1];
				if(j)f[i][j-1][1]=(f[i][j-1][1]+(ll)v*h2[i][j])%mod;
				if(i&&j)f[i-1][j-1][0]=(f[i-1][j-1][0]+(ll)v*w2[i][j]*(h1[i][j-1]-w1[i][j-1]))%mod;
			}
		}
	printf("%d\n",(f[0][0][0]+f[0][0][1])%mod);
	return 0;
}