#include <bits/stdc++.h>
using namespace std;
#define ALL(x) x.begin(),x.end()
#define SIZE(x) (int)x.size()
#define forsn(i,s,n) for(int i=int(s);i<int(n);i++)
#define forn(i,n) for(int i=0;i<int(n);i++)
#define dforsn(i,s,n) for(int i=int(n)-1;i>=int(s);i--)
#define dforn(i,n) for(int i=int(n)-1;i>=0;i--)
#define fst first
#define snd second
typedef long long ll;
typedef pair<int,int> ii;
const int MAXN = 500;
int f[MAXN][MAXN];
pair<ll,ii> dp[MAXN][MAXN][MAXN];
int L[MAXN][MAXN], R[MAXN][MAXN];
int n;
pair<ll,ii> solve(int l, int r, int p) {
if(l>r) return {0LL,{-1,n}};
if(dp[l][r][p].fst!=-1) return dp[l][r][p];
auto calculate = [&](pair<int,ii> prev, int i) {
auto [val,range] = prev;
int mini = max(range.fst, L[i][p]), maxi = min(range.snd, R[i][p]);
return (pair<ll,ii>){val+max(maxi-mini-1,0),{mini,maxi}};
};
return dp[l][r][p] = max(calculate(solve(l+1,r,p),l),
calculate(solve(l,r-1,p),r));
}
int biggest_stadium(int N, vector<vector<int>> F) {
n=N;
forn(i,n) forn(j,n) {
f[i][j]=F[i][j];
}
forn(i,n) forn(j,n) forn(k,n) {
dp[i][j][k].fst = -1;
}
forn(i,n) {
int k=-1;
forn(j,n) {
if(f[i][j]==1) k=j;
L[i][j]=k;
}
}
forn(i,n) {
int k=n;
dforn(j,n) {
if(f[i][j]==1) k=j;
R[i][j]=k;
}
}
ll ans = 0LL;
forn(i,n) {
ll curr = solve(0,n-1,i).fst;
ans=max(ans,curr);
}
return int(ans);
}
