#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 biggest_stadium(int N, vector < vi > F) {
	L(i, 0, N - 1) {
		int l = 0;
		L(j, 0, N) {
			if(j == N || F[i][j] == 1) {
				if(l < j) {
					f[i][i][make_pair(l, j - 1)] = j - l;
				}
				l = j + 1;
			}
			if(j < N) lst[i][j] = l;
		}
	}
	
	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;
						chk(f[nl - 1][nr][make_pair(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;
						chk(f[nl][nr + 1][make_pair(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}});
//}