P736 创意吃鱼法

题目描述
回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中，然后开始思考：到底要以何种方法吃鱼呢（猫猫就是这么可爱，吃鱼也要想好吃法 ^_*）。她发现，把大池子视为01矩阵（0表示对应位置无鱼，1表示对应位置有鱼）有助于决定吃鱼策略。

在代表池子的01矩阵中，有很多的正方形子矩阵，如果某个正方形子矩阵的某条对角线上都有鱼，且此正方形子矩阵的其他地方无鱼，猫猫就可以从这个正方形子矩阵“对角线的一端”下口，只一吸，就能把对角线上的那一队鲜鱼吸入口中。

猫猫是个贪婪的家伙，所以她想一口吃掉尽量多的鱼。请你帮猫猫计算一下，她一口下去，最多可以吃掉多少条鱼？

输入格式
有多组输入数据，每组数据：

第一行有两个整数n和m（n,m≥1），描述池塘规模。接下来的n行，每行有m个数字（非“0”即“1”）。每两个数字之间用空格隔开。

对于30%的数据，有n,m≤100

对于60%的数据，有n,m≤1000

对于100%的数据，有n,m≤2500

输出格式
只有一个整数——猫猫一口下去可以吃掉的鱼的数量，占一行，行末有回车。

输入输出样例
输入 #1 复制
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
输出 #1 复制
3
说明/提示
右上角的

1 0 0
0 1 0
0 0 1

思路：此题类似于最大1矩阵
不过需要预处理
需统计连续的横，纵的零的个数
再与连续对角1取最小值，看能构成的正方形矩阵大小
状态转移方程： 	dp[i][j] = min(dp[i - 1][j - 1], min(heng[i][j - 1], zhong[i - 1][j]))+1 和 dp[i][j] = min(dp[i + 1][j - 1], min(heng[i][j - 1], zhong[i + 1][j]))+1;
前一个是主对角线的正方形，第二个是副对角线的正方形 

#include  <iostream>
#include <algorithm>
#include <cstring>

using namespace std;
#define ms(i,j) memset(i,j,sizeof i);
const int  N = 2510;

int g[N][N];
int dp[N][N];
int heng[N][N];
int zhong[N][N];
int n, m;

int main()
{
	cin >> n >> m;
	int ans=0;
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= m; j++)
		{
			scanf("%d", &g[i][j]);
			if (!g[i][j])//看连续的横的零，纵的零有多少
			{
				heng[i][j]=heng[i][j - 1]+1;
				zhong[i][j] = zhong[i - 1][j]+1;
			}
			else//看01构成的正方形有多大
			{
				dp[i][j] = min(dp[i - 1][j - 1], min(heng[i][j - 1], zhong[i - 1][j]))+1;
				ans = max(ans, dp[i][j]);
			}
		}
	for (int i = 1; i <= n; i++)
	{
		for (int j = 1; j <= m; j++)
			cout << heng[i][j] << " ";
		cout << endl;
	}
	ms(dp, 0);
	ms(heng, 0);
	ms(zhong, 0);
	for(int i=n;i>0;i--)
		for (int j = 1; j <= m; j++)
		{
			if (!g[i][j])
			{
				heng[i][j]= heng[i][j - 1]+1;
				zhong[i][j] = zhong[i +1][j]+1;
			}
			else
			{
				dp[i][j] = min(dp[i + 1][j - 1], min(heng[i][j - 1], zhong[i + 1][j]))+1;
				ans = max(ans, dp[i][j]);
			}
		}
	cout << ans << endl;
	system("pause");
	return 0;
}

传送门