#include<bits/stdc++.h>
#define L(i, j, k) for(int i = (j); i <= (k); i++)
#define R(i, j, k) for(int i = (j); i >= (k); i--)
#define ll long long
#define ull unsigned long long 
#define sz(a) ((int) a.size())
#define vi vector<int>
#define me(a, x) memset(a, x, sizeof(a))
using namespace std;
const int N = 2007;

int n;
int lst[N][N];
map < pair < int, int >, int > f[N][N]; 

void chk(int &x, int y) {
	if(y > x) x = y;
}

int up[N][N], dw[N][N], tms;
void pb(int nl, int nr, int tl, int tr, int ans) {
	++tms;
	if(tms > 6e7) {
		assert(false);
	}
	int clen = nr - nl; 
	nl = *max_element(up[nl] + tl, up[nl] + tr + 1);
	nr = *min_element(dw[nr] + tl, dw[nr] + tr + 1);
	while(nl > 0 && lst[nl - 1][tr] <= tl) --nl;
	while(nr < n - 1 && lst[nr + 1][tr] <= tl) ++nr;
	chk(f[nl][nr][make_pair(tl, tr)], ans + (nr - nl - clen) * (tr - tl + 1));
}

int biggest_stadium(int N, vector < vi > F) {
	n = N;
	L(i, 0, N - 1) {
		int l = 0;
		L(j, 0, N) {
			if(j == N || F[i][j] == 1) l = j + 1;
			if(j < N) lst[i][j] = l;
		}
	}
	L(i, 0, N - 1) {
		L(j, 0, N - 1) {
			if(F[i][j] == 1) up[i][j] = i + 1;
			else up[i][j] = i == 0 ? i : up[i - 1][j];
		}
	}
	R(i, N - 1, 0) {
		L(j, 0, N - 1) {
			if(F[i][j] == 1) dw[i][j] = i - 1;
			else dw[i][j] = i == N - 1 ? i : dw[i + 1][j];
		}
	}
	
	L(i, 0, N - 1) {
		int l = 0;
		L(j, 0, N) {
			if(j == N || F[i][j] == 1) {
				if(l < j) pb(i, i, l, j - 1, j - l);
				l = j + 1;
			}
		}
	}
	
	int ans = 0;
	L(len, 1, N) {
		L(l, 0, N - len) {
			int r = l + len - 1;
			for(auto &o : f[l][r]) {
				ans = max(ans, o.second);
				int tl = o.first.first;
				int tr = o.first.second;
				
				int nl = l;
				int nr = r;
				while(nl > 0 && lst[nl - 1][tr] <= tl) --nl;
				while(nr < N - 1 && lst[nr + 1][tr] <= tl) ++nr;
				
				if(nl != l || nr != r) {
					chk(f[nl][nr][make_pair(tl, tr)], 
						o.second + ((nr - nl) - (r - l)) * (tr - tl + 1));
					continue;
				}
				
				if(nl > 0) {
					for(int r = tr, l; r >= tl; r = l - 2) {
						l = max(tl, lst[nl - 1][r]);
						if(l > r) continue;
						pb(nl - 1, nr, l, r, o.second + r - l + 1);
					}
				}
				
				if(nr < N - 1) {
					for(int r = tr, l; r >= tl; r = l - 2) {
						l = max(tl, lst[nr + 1][r]);
						if(l > r) continue;
						pb(nl, nr + 1, l, r, o.second + r - l + 1);
					}
				}
			}
		}
	}
	return ans;
}

//int main() {
//	cout << biggest_stadium(5, {vi{0, 0, 0, 0, 0}, 
// vi{1, 0, 0, 0, 0}, 
// vi{0, 0, 0, 0, 0}, 
// vi{0, 0, 0, 0, 0}, 
// vi{0, 0, 1, 0, 0}});
//}