#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;
struct rmq {
	int n;
	int bcnt;
	vi lg;
	unsigned int sb[N];
	vi a, suf, pre;
	vector < vi > s;
	rmq(vi p = vi{}) {
		n = sz(p); 
		if(!n) return ;
		a = p;
		suf.resize(n + 1);
		pre.resize(n + 1);
		s.resize(1);
		for(int l = 0, r; l < n; l = r + 1) {
			r = min(l + 31, n - 1);
			L(i, l, r) pre[i] = suf[i] = a[i];
			L(i, l + 1, r) pre[i] = max(pre[i], pre[i - 1]);
			R(i, r - 1, l) suf[i] = max(suf[i], suf[i + 1]);
			s[0].emplace_back(pre[r]);
			unsigned int cur = 0;
			vi stk;
			L(i, 0, r - l) {
				while(sz(stk) && a[stk.back() + l] < a[i + l]) 
					cur -= 1u << stk.back(), stk.pop_back();
				stk.emplace_back(i), cur += 1u << i, sb[l + i] = cur;
			}
		}
		bcnt = sz(s[0]);
		lg.resize(bcnt + 1);
		L(i, 2, bcnt) lg[i] = lg[i >> 1] + 1; 
		s.resize(lg[bcnt] + 1);
		L(i, 1, lg[bcnt]) s[i].resize(bcnt);
		L(i, 1, lg[bcnt]) 
			L(j, 0, bcnt - (1 << i)) 
				s[i][j] = max(s[i - 1][j], s[i - 1][j + (1 << (i - 1))]);
	}
	
	inline int qry(int l, int r) {
		int bl = l / 32;
		int br = r / 32;
		if(bl == br) {
			return a[l + __builtin_ctz(sb[r] >> (l - bl * 32))];
		}
		int ans = max(suf[l], pre[r]);
		++bl, --br;
		if(bl <= br) {
			int p = lg[br - bl + 1];
			ans = max(ans, max(s[p][bl], s[p][br - (1 << p) + 1]));
		}
		return ans;
	}
} U[N], D[N];

void pb(int nl, int nr, int tl, int tr, int ans) {
	++tms;
	if(tms > 2e7) {
		assert(false);
	}
	int clen = nr - nl; 
	nl = U[nl].qry(tl, tr); 
	nr = -D[nr].qry(tl, tr);
	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];
		}
		vi a(N);
		L(j, 0, n - 1) a[j] = up[i][j];
		U[i] = rmq(a);
	}
	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];
		}
		vi a(N);
		L(j, 0, n - 1) a[j] = -dw[i][j];
		D[i] = rmq(a);
	}
	
	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;
}