猫猫把鱼全部转移到了她那长方形大池子中,然后开始思考:到底要以何种方法吃鱼呢(猫猫就是这么可爱,吃鱼也要想好吃法 ^_*)。她发现,把大池子视为01矩阵(0表示对应位置无鱼,1表示对应位置有鱼)有助于决定吃鱼策略。在代表池子的01矩阵中,有很多的正方形子矩阵,如果某个正方形子矩阵的某条对角线上都有鱼,且此正方形子矩阵的其他地方无鱼,猫猫就可以从这个正方形子矩阵“对角线的一端”下口,只一吸,就能把对角线上的那一队鲜鱼吸入口中。
猫猫是个贪婪的家伙,所以她想一口吃掉尽量多的鱼。请你帮猫猫计算一下,她一口下去,最多可以吃掉多少条鱼?
Input
有多组输入数据,每组数据:
第一行有两个整数n和m(n,m≥1),描述池塘规模。接下来的n行,每行有m个数字(非“0”即“1”)每两个数字之间用空格隔开。
Output
只有一个整数——猫猫一口下去可以吃掉的鱼的数量,占一行,行末有回车。
Sample Input
4 6
0 1 0 1 0 0
0 0 1 0 1 0
1 1 0 0 0 1
0 1 1 0 1 0
Sample Output
3
Hint
对于30%的数据,有n,m≤100
对于60%的数据,有n,m≤1000
对于100%的数据,有n,m≤2500
分析:
F[i][j]表示以(i,j)为对角顶点(map[i][j]=1)的一个正方形的对角线的最小值.
up[i][j]表示(i,j)上面(包含(i,j))连续的一段0的个数。
left[i][j]表示(i,j)左边(包含(i,j))连续的一段0的个数。
right[i][j]表示(i,j)右边(包含(i,j))连续的一段0的个数。
因为对角线有两条,所以要分左对角线和右对角线分别DP,answer=max(F[i][j])
【参考程序】:
#include<cstring>
#include<cstdio>
using namespace std;
typedef int array[2510][2510];
array map,up,left,right,F;
int n,m,ans;
void yuchuli()
{
memset(up,0,sizeof(up));
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
if (map[i][j]==0) up[i][j]=up[i-1][j]+1;
memset(left,0,sizeof(left));
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
if (map[i][j]==0) left[i][j]=left[i][j-1]+1;
memset(right,0,sizeof(right));
for (int i=1;i<=n;i++)
for (int j=m;j>=1;j--)
if (map[i][j]==0) right[i][j]=right[i][j+1]+1;
}
int Min(int x,int y,int z)
{
int t=x;
if (y<x) return y;
if (z<x) return z;
return t;
}
void solve()
{
ans=0;
memset(F,0,sizeof(F));
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
if (map[i][j]==1)
{
F[i][j]=Min(F[i-1][j-1],left[i][j-1],up[i-1][j])+1;
if (ans<F[i][j]) ans=F[i][j];
}
memset(F,0,sizeof(F));
for (int i=1;i<=n;i++)
for (int j=m;j>=1;j--)
if (map[i][j]==1)
{
F[i][j]=Min(F[i-1][j+1],up[i-1][j],right[i][j+1])+1;
if (ans<F[i][j]) ans=F[i][j];
}
}
int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
scanf("%d",&map[i][j]);
yuchuli();
solve();
printf("%d\n",ans);
return 0;
}