#include<bits/stdc++.h>
using namespace std;
const int N=3e3+5;
const int p=998244353;
int n;
char grid[N][N];
int val1[N][N][2],val2[N][N][2];
int dp[N][N][3];
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>>grid[i][j];
        }
        grid[i][n+1]='1';
    }
    for(int i=1;i<=n;i++){
        val1[i][0][1]=1;
        for(int j=1;j<=n;j++){
            if(grid[i][j]=='0' || grid[i][j]=='?'){
                val1[i][j][0]=(val1[i][j-1][0]+val1[i][j-1][1])%p;
            }
            if(grid[i][j]=='1' || grid[i][j]=='?'){
                val1[i][j][1]=val1[i][j-1][1];
            }
        }
    }
    for(int j=1;j<=n+1;j++){
        val2[0][j][1]=1;
        for(int i=1;i<=n;i++){
            if(grid[i][j]=='0' || grid[i][j]=='?'){
                val2[i][j][0]=(val2[i-1][j][0]+val2[i-1][j][1])%p;
            }
            if(grid[i][j]=='1' || grid[i][j]=='?'){
                val2[i][j][1]=val2[i-1][j][1];
            }
        }
    }
    for(int i=0;i<=n;i++){
        dp[i][1][0]=val2[i][1][0];
        dp[i][1][1]=val2[i][1][1];
        for(int j=2;j<=n+1;j++){
            if(i>=1)dp[i][j][2]=(1ll*(val1[i][j-1][0]+val1[i][j-1][1])%p*dp[i-1][j][2]%p+1ll*val1[i][j-1][1]*dp[i-1][j-1][0]%p)%p;
            dp[i][j][0]=1ll*((dp[i][j-1][0]+dp[i][j-1][1])%p+dp[i][j][2])%p*val2[i][j][0]%p;
            dp[i][j][1]=1ll*((dp[i][j-1][0]+dp[i][j-1][1])%p+dp[i][j][2])%p*val2[i][j][1]%p;
        }
    }
    cout<<dp[n][n+1][1];
    return 0;
}