#include <algorithm>
#include <bitset>
#include <cstdint>
#include <cstring>
#include <iostream>
#include <limits.h>
#include <math.h>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <vector>
//#define int long long
#define loop(X, N) for(int X = 0; X < (N); X++)
#define all(V) V.begin(), V.end()
#define rall(V) V.rbegin(), V.rend()
using namespace std;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef pair<int, int> ii;
typedef vector<ii> vii;
typedef vector<vector<ii>> vvii;
typedef vector<bool> vb;
typedef vector<vector<bool>> vvb;
int biggest_stadium(int n, vvi f) {
vvi up(n, vi(n)), down(n, vi(n)), left(n, vi(n)), right(n, vi(n));
loop(x, n) {
up[0][x] = 1 - f[0][x];
for (int y = 1; y < n; y++) {
if (f[y][x])
up[y][x] = 0;
else
up[y][x] = up[y - 1][x] + 1;
}
down[n - 1][x] = 1 - f[n - 1][x];
for (int y = n - 2; y >= 0; y--) {
if (f[y][x])
down[y][x] = 0;
else
down[y][x] = down[y + 1][x] + 1;
}
}
loop(y, n) {
left[y][0] = 1 - f[y][0];
for (int x = 1; x < n; x++) {
if (f[y][x])
left[y][x] = 0;
else
left[y][x] = left[y][x - 1] + 1;
}
right[y][n - 1] = 1 - f[y][n - 1];
for (int x = n - 2; x >= 0; x--) {
if (f[y][x])
right[y][x] = 0;
else
right[y][x] = right[y][x + 1] + 1;
}
//loop(y, n)
// cout << high[x][y] << low[x][y] << ' ';
//cout << endl;
}
//cout << endl;
int result = 0;
loop(y, n) {
loop(x, n) {
if (f[y][x]) {
//cout << "00 ";
continue;
}
int count1 = 0;
{
int l = 0, r = INT_MAX;
vii intervals;
for (int y2 = y; y2 < y + down[y][x]; y2++) {
l = max(l, x - left[y2][x] + 1);
r = min(r, x + right[y2][x] - 1);
intervals.push_back({ l, r });
count1 += r - l + 1;
}
l = x - left[y][x] + 1;
r = x + right[y][x] - 1;
for (int y2 = y - 1; y2 > y - up[y][x]; y2--) {
l = max(l, x - left[y2][x] + 1);
r = min(r, x + right[y2][x] - 1);
for (auto [li, ri] : intervals) {
if (r < ri)
l = max(l, li);
if (l > li)
r = min(r, ri);
}
if (r < l)
break;
count1 += r - l + 1;
}
}
int count2 = 0;
{
int l = 0, r = INT_MAX;
vii intervals;
for (int y2 = y; y2 > y - up[y][x]; y2--) {
l = max(l, x - left[y2][x] + 1);
r = min(r, x + right[y2][x] - 1);
intervals.push_back({ l, r });
count2 += r - l + 1;
}
l = x - left[y][x] + 1;
r = x + right[y][x] - 1;
for (int y2 = y + 1; y2 < y + down[y][x]; y2++) {
l = max(l, x - left[y2][x] + 1);
r = min(r, x + right[y2][x] - 1);
for (auto [li, ri] : intervals) {
if (r < ri)
l = max(l, li);
if (l > li)
r = min(r, ri);
}
if (r < l)
break;
count2 += r - l + 1;
}
}
int count = max(count1, count2);
result = max(result, count);
//cout << count / 10 << count % 10 << ' ';
}
//cout << endl;
}
return result;
}
