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

#define rep(i, N) for (int i = 0; i < (N); i++)
#define all(a) (a).begin(), (a).end()
#define pb push_back

using ll = long long;

const int MOD = 998244353;

short rmin[3010][3010];
short rmax[3010][3010];
short cmin[3010][3010];
short cmax[3010][3010];

int row(int i, int j) {
    return max(rmax[i][j] - rmin[i][j] + 1, 0);
}

int col(int i, int j) {
    return max(cmax[i][j] - cmin[i][j] + 1, 0);
}

int solve(vector<string> a) {
    int N = a.size();

    for (int i = 1; i <= N; i++) {
        for (int j = 1; j <= N; j++) {
            rmin[i][j] = rmin[i][j - 1];
            cmin[i][j] = cmin[i - 1][j];
            if (a[i - 1][j - 1] == '1') {
                rmin[i][j] = j;
                cmin[i][j] = i;
            }
        }
    }

    for (int i = 1; i <= N; i++) {
        int z = INT_MAX;
        for (int j = 1; j <= N; j++) {
            if (a[i - 1][j - 1] == '0') {
                z = min(z, j - 1);
            }
            rmax[i][j] = min(j, z);
        }
    }

    for (int j = 1; j <= N; j++) {
        int z = INT_MAX;
        for (int i = 1; i <= N; i++) {
            if (a[i - 1][j - 1] == '0') {
                z = min(z, i - 1);
            }
            cmax[i][j] = min(i, z);
        }
    }

    vector<int> dp(N + 1), ep(N + 1);
    dp[0] = ep[0] = 1;

    for (int i = 1; i <= N; i++) {
        vector<int> _dp(N + 1), _ep(N + 1);
        _dp[0] = _ep[0] = 1;
        for (int j = 1; j <= N; j++) {
            _ep[j] = ((ll)_ep[j - 1] * col(i - 1, j) + dp[j]) % MOD;
            _dp[j] = (ll)dp[j] * row(i, j) % MOD;
            if (j <= rmax[i][j] && (i >= 2 && a[i - 2][j - 1] != '1')) {
                _dp[j] = (_dp[j] + (ll)_ep[j - 1] * col(i - 2, j)) % MOD;
            }
        }
        dp = _dp;
        ep = _ep;
    }

    return dp[N];
}

int main() {
    int N;
    cin >> N;
    vector<string> a(N);
    rep(i, N) {
        cin >> a[i];
    }
    rep(i, N) {
        a[i].pb('0');
    }
    a.pb(string(N + 1, '1'));
    N++;

    cout << solve(a) << endl;
}