bzoj1057 [ZJOI2007]棋盘制作

Description

---

　　国际象棋是世界上最古老的博弈游戏之一，和中国的围棋、象棋以及日本的将棋同享盛名。据说国际象棋起源
于易经的思想，棋盘是一个8*8大小的黑白相间的方阵，对应八八六十四卦，黑白对应阴阳。而我们的主人公小Q，
正是国际象棋的狂热爱好者。作为一个顶尖高手，他已不满足于普通的棋盘与规则，于是他跟他的好朋友小W决定
将棋盘扩大以适应他们的新规则。小Q找到了一张由N*M个正方形的格子组成的矩形纸片，每个格子被涂有黑白两种
颜色之一。小Q想在这种纸中裁减一部分作为新棋盘，当然，他希望这个棋盘尽可能的大。不过小Q还没有决定是找
一个正方形的棋盘还是一个矩形的棋盘（当然，不管哪种，棋盘必须都黑白相间，即相邻的格子不同色），所以他
希望可以找到最大的正方形棋盘面积和最大的矩形棋盘面积，从而决定哪个更好一些。于是小Q找到了即将参加全
国信息学竞赛的你，你能帮助他么？

Input

---

　　第一行包含两个整数N和M，分别表示矩形纸片的长和宽。接下来的N行包含一个N * M的01矩阵，表示这张矩形
纸片的颜色（0表示白色，1表示黑色）。

Output

---

　　包含两行，每行包含一个整数。第一行为可以找到的最大正方形棋盘的面积，第二行为可以找到的最大矩形棋
盘的面积（注意正方形和矩形是可以相交或者包含的）。

HINT

---

N, M ≤ 2000

Solution

---

显然的可以想到隔一个异或1，这样就变成了求最大的全相同的子矩阵。这里通过记录一个点向上、向左、向右分别能扩展最多多少个，递推维护一下枚举答案即可

第一问的正方形一定在最大矩阵里面

Code

---

#include <stdio.h>
#include <string.h>
#include <math.h>
#define rep(i, st, ed) for (int i = st; i <= ed; i += 1) 
#define drp(i, st, ed) for (int i = st; i >= ed; i -= 1)
#define fill(x, t) memset(x, t, sizeof(x))
#define max(x, y) ((x)>(y)?(x):(y))
#define min(x, y) ((x)<(y)?(x):(y))
#define N 2005
int a[N][N], u[N][N], l[N][N], r[N][N];
int n, m, ans1, ans2;
int read() {
    int x = 0, v = 1; char ch = getchar();
    for (; ch < '0' || ch > '9'; v = ((ch == '-')?(-1):(v)), ch = getchar());
    for (; ch <= '9' && ch >= '0'; x = x * 10 + ch - '0', ch = getchar());
    return x * v;
}
int sqr(int x) {
    return x * x;
}
void solve() {
    rep(i, 1, n) {
        rep(j, 1, m) {
            l[i][j] = l[i][j - 1] * (a[i][j - 1] == a[i][j]) + 1;
        }
        drp(j, m, 1) {
            r[i][j] = r[i][j + 1] * (a[i][j + 1] == a[i][j]) + 1;
        }
    }
    rep(i, 1, m) {
        u[1][i] = 1;
    }
    rep(i, 1, n - 1) {
        rep(j, 1, m) {
            u[i + 1][j] = u[i][j] * (a[i + 1][j] == a[i][j]) + 1;
            if (a[i + 1][j] == a[i][j]) {
                l[i + 1][j] = min(l[i][j], l[i + 1][j]);
                r[i + 1][j] = min(r[i][j], r[i + 1][j]);
            }
        }
    }
    rep(i, 1, n) {
        rep(j, 1, m) {
            ans1 = max(ans1, sqr(min(u[i][j], (l[i][j] + r[i][j] - 1))));
            ans2 = max(ans2, u[i][j] * (l[i][j] + r[i][j] - 1));
        }
    }
}
int main(void) {
    n = read();
    m = read();
    fill(a, 31);
    rep(i, 1, n) {
        rep(j, 1, m) {
            a[i][j] = read();
        }
    }
    rep(i, 1, n) {
        for (int j = 1 + (i & 1); j <= m; j += 2) {
            a[i][j] ^= 1;
        }
    }
    ans1 = 0;
    ans2 = 0;
    solve();
    printf("%d\n%d\n", ans1, ans2);
    return 0;
}