动态规划：最大矩阵寻找

例题选自洛谷P1387 最大正方形寻找
题目要求在一个矩阵中寻找一个全为1的正方形，要求可能的最大的正方形的边长。

这个题如果用暴力搜索既不直观实现起来也不方便。我们用动态规划来求解。作为一个很典型的题型，看到正方形、矩阵时，应该要想到右下角坐标的记录作用。那么，对于右下角点的状态可以直接影响的是它左边、上边和左上边的三个点。

我们可以设 dp[i][j] 为以坐标（i，j）为右下方边界点的正方形的最大边长。那么，对于（i，j）这个点而言，我们要把以（i-1，j-1）为右下角的正方形的部分尽量取完，最理想的情况就是dp[i][j]=dp[i-1][j-1]+1。但是由于与（i，j）同行和同列的边不一定和dp[i-1][j-1]的值相同，所以不一定能够满足最好的情况。

我们只需要取同行最大边长、同列最大边长以及dp[i-1][j-1]里的最小值，再加上1，就是该点可以取到的最大边长。同行最大边长即为dp[i][j-1]，同列最大边长即为dp[i-1][j]。

那么，状态转移方程即为：

if(num==1)  f[i][j]=min(min(f[i][j-1],f[i-1][j]),f[i-1][j-1]);

由于我们只需要某一个点左上方的状态，所以可以一边读取数据一边转状态。

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
int f[105][105];
int main()
{
    int n,m;
    cin>>n>>m;
    int ans=-1;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            int num;
            cin>>num;
            if(num==1)f[i][j]=min(min(f[i][j-1],f[i-1][j]),f[i-1][j-1])+1;
            ans=max(ans,f[i][j]);
        }
    }
    cout<<ans;
}

这里还有一道变式题洛谷P1681
和本题不同的是，对于每一个点，可以转移状态的标准是，这一点的值必须和其上方以及左方的值不同，和其左上方的值相同，这样才确保了正方形内黑白相间。否则更新dp[i][j]=1;

代码如下：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
int f[1505][1505];
int num[1505][1505];
int main()
{
    int n,m;
    cin>>n>>m;
    int ans=-1;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            cin>>num[i][j];
            if(i==1){f[i][j]=1;ans=max(ans,f[i][j]);continue;}
            if(j==1){f[i][j]=1;ans=max(ans,f[i][j]);continue;}
            int a=num[i][j],b=num[i][j-1],c=num[i-1][j],d=num[i-1][j-1];            
            if(a==d&&a!=c&&a!=b)f[i][j]=min(min(f[i][j-1],f[i-1][j]),f[i-1][j-1])+1;
            else f[i][j]=1;
            ans=max(ans,f[i][j]);
        }
    }
    cout<<ans;
}

正方形的右下角在类似的题型中是问题的关键，在某些题中也会有类似于矩阵内是随机的1-9的数，找数字相同的最大矩形，那么就可以如上建立状态转移方程，只需要添加一些限制条件即可。或者也可建立一个三维数组，那么记录每一位上1-9的数字分别组成的正方形的最大边长，就可以直接像第一个例题一样写出代码了。