BZOJ 1057 棋盘制作 极大子矩阵（悬线法）

Description

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

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

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

3 3

1 0 1

0 1 0

1 0 0
Sample Output

4

6
HINT

N, M ≤ 2000

分析：这道题就是求极大子矩阵的裸题，首先求极大子矩阵有2种方案，1是枚举障碍点来不断扩展求出最大子矩阵，复杂度与障碍点的个数有关。适用于障碍点稀疏的图。2是悬线法，枚举每条悬线，因为悬线包含的矩阵集合一定包含极大子矩阵，于是可以枚举悬线，因为悬线只有（n-1）*m条所以复杂度较为稳定，n*m，适用于障碍点稠密的图。
悬线法大致是定义三个数组：h[i][j],l[i][j],r[i][j]，h表示能扩展的最大高度，l表示能向左扩展的最远点，r表示能向右扩展的最远点。逐步向下地推，若下面的格子于当前格子颜色相同，则h加1，并且l=max（li-1），r=min（ri-1），因为你要满组上面所有的矩阵都能与下面格子组成新的矩阵，所以l取max，r取min，画个图就好了，最后每次更新长度，得到答案。

# include <cstdio>
# include <iostream>
# include <cstring>
# include <algorithm>
# include <cmath>
# include <vector>
using namespace std;
int read()
{
    register int f=1,i=0;char c=getchar();
    while(c<'0'||c>'9') {if(c=='-'){f=-1;}c=getchar();}
    while(c>='0'&&c<='9') {i=(i<<3)+(i<<1)+c-'0';c=getchar();}
    return i*f;
}
const int N=2005;
int n,m,f[N][N],l[N][N],r[N][N],h[N][N],a[N][N];
inline void workzheng()
{
    int Ans=0;
    for(int i=1;i<=n;++i) 
        for(int j=1;j<=m;++j) f[i][j]=1;
    for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j)
        {
            if(a[i][j]==a[i-1][j]&&a[i-1][j]==a[i-1][j-1]&&a[i][j]==a[i][j-1]) 
                f[i][j]=min(f[i-1][j],min(f[i][j-1],f[i-1][j-1]))+1;
            Ans=max(f[i][j],Ans);
        }
    printf("%d\n",Ans*Ans);
}
inline void workchang()
{
    for(int i=1;i<=n;++i)
        for(int j=2;j<=m;++j)
            if(a[i][j-1]==a[i][j]) l[i][j]=l[i][j-1];
    for(int i=1;i<=n;++i)
        for(int j=m-1;j;--j)
            if(a[i][j+1]==a[i][j]) r[i][j]=r[i][j+1];
    int Ans=0; 
    for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j)
        {
            if(i>1&&a[i][j]==a[i-1][j])
            {
                l[i][j]=max(l[i][j],l[i-1][j]);
                r[i][j]=min(r[i][j],r[i-1][j]);
                h[i][j]=h[i-1][j]+1;
            }
            int len=r[i][j]-l[i][j]+1;
            Ans=max(Ans,len*h[i][j]);
        }
    printf("%d\n",Ans);
}
int main()
{
    freopen("lx.in","r",stdin);
    n=read(),m=read();
    for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j) 
        {
            a[i][j]=read();
            if((j+i)&1) a[i][j]=!a[i][j];
            l[i][j]=j;r[i][j]=j;h[i][j]=1;
        }
    workzheng();workchang();
}

第一次写时第二个问有点懵，于是分开写的。第一个work就是求最大正方形矩阵，第二个work就是悬线法了。