LGOJ1169-[ZJOI2007]棋盘制作（DP+悬线法）

写在前面

由于我太菜了，不知道悬线法，所以自己瞎搞一通，$WA$了，于是去抄学习了题解，所以才有了这篇博客。

悬线法

它是什么？

这是一种$DP$的思路，主要适用于在棋盘 $Or$ 矩阵上$DP$。

它能干什么？

解决给定矩阵中满足特定条件（具体看题目）的最大子矩阵

怎么用？

一般来说，用一条线左右移动$Or$上下移动直到不满足约束条件或者到达边界。

Description

传送门
国际象棋是世界上最古老的博弈游戏之一，和中国的围棋、象棋以及日本的将棋同享盛名。据说国际象棋起源于易经的思想，棋盘是一个$8\times 8$大小的黑白相间的方阵，对应八八六十四卦，黑白对应阴阳。

而我们的主人公小Q，正是国际象棋的狂热爱好者。作为一个顶尖高手，他已不满足于普通的棋盘与规则，于是他跟他的好朋友小W决定将棋盘扩大以适应他们的新规则。

小Q找到了一张由$N\times M$个正方形的格子组成的矩形纸片，每个格子被涂有黑白两种颜色之一。小Q想在这种纸中裁减一部分作为新棋盘，当然，他希望这个棋盘尽可能的大。

不过小Q还没有决定是找一个正方形的棋盘还是一个矩形的棋盘（当然，不管哪种，棋盘必须都黑白相间，即相邻的格子不同色），所以他希望可以找到最大的正方形棋盘面积和最大的矩形棋盘面积，从而决定哪个更好一些。

于是小Q找到了即将参加全国信息学竞赛的你，你能帮助他么？

输入输出格式

输入格式：
包含两个整数$N$和$M$，分别表示矩形纸片的长和宽。接下来的$N$行包含一个$N\times M$的0101矩阵，表示这张矩形纸片的颜色（$0$表示白色，$1$表示黑色）。

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

Solution

都说了是悬线法了，那么定义$DP$状态也就很容易了。
本题中，我们需要计算一个矩形的面积，所以说我们就直接考虑算长和宽，宽的计算就直接考虑从一个点$(i,j)$出发，向上能够到达的最远距离，设为$up[i][j]$。
长就可以用向右能够走到的最远距离减去向左走到的最远距离即可。
所以说我们再维护两个状态：$left[i][j]$表示向左走到的最远距离，$right[i][j]$表示享有走到的最远距离。
那么就有如下方程：
$$left\left[i\right]\left[j\right]=max\left(left\left[i\right]\left[j\right],left\left[i-1\right]\left[j\right]\right)$$
$$right\left[i\right]\left[j\right]=min\left(right\left[i\right]\left[j\right],right\left[i-1\right]\left[j\right]\right)$$
注意到这个方程是考虑了上一层的状态的，为什么呢？因为我们最终是需要计算一个矩形的面积，所以说肯定是需要将当前的状态和上一行的比较来更新$left$和$right$还有$up$。（主要是$up$）
然后就可以随便搞这道题目了。
代码如下：

#include<bits/stdc++.h>
using namespace std;
template<typename T> 
void read(T &x){
    T neg=1; x=0;
    static char c=getchar();
    while(!isdigit(c)) c=getchar();
    if(c=='-') neg=-1,c=getchar();
    while(isdigit(c)) x=(x<<1)+(x<<3)+c-'0',c=getchar();
    x*=neg;
}
int mat[2005][2005],Right[2005][2005],Left[2005][2005],up[2005][2005];
int n,m;
int anssq=-0x7fffffff,ansma=-0x7fffffff;
int main(){
    read(n); read(m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            read(mat[i][j]);
            Left[i][j]=Right[i][j]=j;
            up[i][j]=1;
        }
    for(int i=1;i<=n;i++)
        for(int j=2;j<=m;j++)
            if(mat[i][j]!=mat[i][j-1])
                Left[i][j]=Left[i][j-1];
    for(int i=1;i<=n;i++)
        for(int j=m-1;j>=1;j--)
            if(mat[i][j]!=mat[i][j+1])
                Right[i][j]=Right[i][j+1];
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            if(i>1&&mat[i][j]!=mat[i-1][j]){
                Left[i][j]=max(Left[i][j],Left[i-1][j]);
                Right[i][j]=min(Right[i][j],Right[i-1][j]);
                up[i][j]=up[i-1][j]+1;
            }
            int a=Right[i][j]-Left[i][j]+1;
            int b=min(a,up[i][j]);
            anssq=max(anssq,b*b);
            ansma=max(ansma,a*up[i][j]);
        }
    printf("%d\n%d",anssq,ansma);
    return 0;
}