#include <bits/stdc++.h>
using namespace std;
int a[510][510];
const long long mod = 1e9 + 7;
long long f[510][510][2];
long long dp(int l, int r, int type) {
  long long &ans = f[l][r][type];
  if (~ans) return ans;
  ans = 0;
  int i;
  if (!type) {
    if (r - l < 2) return ans = r - l ^ 1;
    for (i = l + 1; i < r; i++) {
      ans += dp(l, i, 1) * (dp(i, r, 0) + dp(i, r, 1)) % mod;
    }
    return ans %= mod;
  } else {
    if (!a[l][r]) return ans = 0;
    if (r - l < 2) return ans = r - l;
    for (i = l; i < r; i++) {
      ans += (dp(l, i, 0) + dp(l, i, 1)) * (dp(i + 1, r, 0) + dp(i + 1, r, 1)) %
             mod;
    }
    return ans %= mod;
  }
}
int main() {
  memset(f, 255, sizeof(f));
  ios::sync_with_stdio(0);
  cin.tie(0), cout.tie(0);
  int n;
  cin >> n;
  int i, j;
  for (i = 1; i <= n; i++) {
    for (j = 1; j <= n; j++) {
      cin >> a[i][j];
    }
  }
  cout << (dp(1, n, 0) + dp(1, n, 1)) % mod << endl;
  return 0;
}
