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

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef long long lint;

const int MAXN = 3000;

int N;
char A[MAXN+10][MAXN+10];

const int MOD = 998244353;
struct mint
{
    int x;
    mint() : x(0) {}
    mint(int _x) : x(_x) {}
    mint operator + (int ot) const { return x+ot<MOD ? x+ot : x+ot-MOD; }
    mint operator - (int ot) const { return x<ot ? x+MOD-ot : x-ot; }
    mint operator * (int ot) const { return 1ll*x*ot%MOD; }
    mint operator += (int ot) { return *this = *this + ot; }
    mint operator -= (int ot) { return *this = *this - ot; }
    mint operator *= (int ot) { return *this = *this * ot; }
    operator int() const { return x; }
};

mint dp[MAXN+10][2];
mint dp3[MAXN+10][MAXN+10], dp4[MAXN+10][MAXN+10]; 
mint dp5[MAXN+10][MAXN+10], dp6[MAXN+10][MAXN+10]; 
mint dp1[MAXN+10][MAXN+10], dp2[MAXN+10][MAXN+10];

int main()
{
    ios_base::sync_with_stdio(false); cin.tie(NULL);

    cin >> N;
    for(int i=1; i<=N; i++)
    {
        string S;
        cin >> S;
        for(int j=0; j<N; j++) A[i][j+1]=S[j];
    }

    for(int i=1; i<=N; i++)
    {
        dp[0][0]=1; dp[0][1]=0;
        dp3[i][0]=1;
        for(int j=1; j<=N; j++)
        {
            dp[j][0]=0; dp[j][1]=0;
            if(A[i][j]!='0') dp[j][0]+=dp[j-1][0];
            if(A[i][j]!='1') dp[j][1]+=dp[j-1][0]+dp[j-1][1];
            dp3[i][j]=dp[j][0]+dp[j][1];
            dp5[i][j]=dp[j][1];
        }
    }

    for(int i=1; i<=N; i++)
    {
        dp[0][0]=1; dp[0][1]=0;
        dp4[0][i]=1;
        for(int j=1; j<=N; j++)
        {
            dp[j][0]=0; dp[j][1]=0;
            if(A[j][i]!='0') dp[j][0]+=dp[j-1][0];
            if(A[j][i]!='1') dp[j][1]+=dp[j-1][0]+dp[j-1][1];
            dp4[j][i]=dp[j][0]+dp[j][1];
        }
    }

    for(int i=1; i<=N; i++)
    {
        dp6[0][i]=1;
        for(int j=1; j<=N; j++)
        {
            if(A[j][i]=='0') dp6[j][i]=0;
            else dp6[j][i]=dp6[j-1][i];
        }
    }

    for(int i=0; i<=N; i++) dp1[0][i]=dp1[i][0]=1;
    for(int i=1; i<=N; i++)
    {
        for(int j=1; j<=N; j++)
        {
            dp2[i][j]=dp2[i][j-1]*dp4[i][j]+dp1[i-1][j-1]*dp5[i][j-1]*dp6[i][j];
            dp1[i][j]=dp1[i-1][j]*dp3[i][j]+dp2[i][j];

            // cout << dp2[i][j] << " ";
        }
        // cout << "\n";
    }
    cout << dp1[N][N] << "\n";
}