P1169 [ZJOI2007]棋盘制作 - 单调栈

这题真是难想。。。按奇偶性分类其实是套路。。。
把横纵坐标和为奇数的^=1一下，然后满足要求的棋盘要么全1要么全0，然后求解答案，把矩阵全^=1一遍再求一次答案就好了，两次求解总有一次棋盘是全0
对每个点预处理以下他往右最多有几个0，枚举一下矩形左边界所在的列，然后枚举矩形下边界的行，用预处理的东西往上跑一下取个可行的，但是这样有点慢，可以用单调栈维护一下每个点的预处理值。。。因为当栈里进了个小于栈顶的值，那么再往下的满足要求的矩阵肯定没有栈顶那么“长/宽”

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
#define debug(x) cerr << #x << "=" << x << endl;
const int MAXN = 2000 + 10;
int n,w,a[MAXN][MAXN],m,s[MAXN][MAXN],ans1,ans2,top,qx[MAXN],q[MAXN];

int calc(int a) {
    return a*a;
}

void update(int s, int line) {
    for(int i=line; top && qx[top] > s; line=q[top--]) {
        ans1 = max(ans1, calc(min(i-q[top], qx[top])));
        ans2 = max(ans2, (i-q[top]) * qx[top]);
    }
    q[++top] = line, qx[top] = s;
}

void solve() {
    for(int i=1; i<=n; i++) {
        for(int j=m; j; j--) {
            s[i][j] = !a[i][j] ? s[i][j+1] + 1 : 0;
        }
    }
    for(int j=1; j<=m; j++) {
        for(int i=1; i<=n; i++) {
            update(s[i][j], i);
        }
        update(0, n+1);
        top = 0;
    }
}

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", ans1, ans2);
    return 0;
}