根据奇偶性，把x，y坐标加起来是奇数的位置翻转一下，就变成求最大01子矩阵了。
然后预处理出每个位置往上最长的01路径，每行用一个单调栈解决。
单调栈从栈底到栈顶维护每一个矩阵的高度和宽度，按高度递增维护。
当前高度如果小于栈顶的高度，那就将栈顶的矩阵更新答案，并且出栈，但是它的宽度是可以为当前高度的矩阵做贡献的，就是当前要加入栈的矩阵的宽度得加上这个要出栈的矩阵。
最后一列做完之后再清一下栈就行了。
#include <bits/stdc++.h> #define pii pair<int, int> #define fi first #define se second const int N = 2222; int a[N][N], dp[N][N], n, m, ans1, ans2; std::pii st[N]; void solve() { for (int i = 1; i <= n; i++) for (int j = 1; j <= m; j++) dp[i][j] = (a[i][j] ? dp[i - 1][j] + 1 : 0); for (int i = 1; i <= n; i++) { int top = 0; for (int j = 1; j <= m + 1; j++) { int w = 0; while (top && st[top].fi > dp[i][j]) { int h = st[top].fi; w += st[top].se; top--; ans1 = std::max(ans1, std::min(w, h) * std::min(w, h)); ans2 = std::max(ans2, w * h); } st[++top] = std::pii(dp[i][j], w + 1); } } } int main() { scanf("%d%d", &n, &m); for (int i = 1; i <= n; i++) for (int j = 1; j <= m; j++) { scanf("%d", a[i] + j); if ((i + j) % 2) a[i][j] ^= 1; } solve(); for (int i = 1; i <= n; i++) for (int j = 1; j <= m; j++) a[i][j] ^= 1; solve(); printf("%d\n%d\n", ans1, ans2); return 0; }