（dp）洛谷P1736创意吃鱼法

洛谷P1736创意吃鱼法

思路：

他要求一个矩阵只能并且必须在对角线上有1才能被吃掉。
如下：

1 0 0
0 1 0
0 0 1

我们用dp[i][j]表示这点满足矩阵的最大值。
那么看下图：

1 0 0 0
0 1 0 1
0 0 1 0
0 0 0 1

可以看出dp[3][3]=3，那怎么求dp[4][4]=2。他显然不是dp[3][3]+1=4≠2。我们观察可以发现，他是被（2，4）的1给限制了，（4，4）的上方只有一个连续的0。
我们可以看出在判断一个点的时候需要关注他上方和左（右）边连续的0。
用row表示左（右）边连续的0的个数，col表示上面连续0的个数。
可以得到左上——右下的转移方程为：
dp[i][j]=min(dp[i-1][j-1],min（row[i][j-1],col[i-1][j]))+1。
右上——左下的转移方程为：
dp[i][j]=min(dp[i-1][j-1],min(dow[i][j+1],col[i-1][j]))+1。

代码：

#include<bits/stdc++.h>
#define pii pair<int,int>
#define ll long long
#define cl(x,y) memset(x,y,sizeof(x))
#define ct cerr<<"Time elapsed:"<<1.0*clock()/CLOCKS_PER_SEC<<"s.\n";
const int N=2550;
const int mod=1e7+9;
const int maxn=0x3f3f3f3f;
const int minn=0xc0c0c0c0;
const int inf=99999999;
using namespace std;
int a[N][N],dp[N][N]={0},row[N][N]={0},col[N][N]={0};
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	int n,m,i,j,res=0;
	cin>>n>>m;
	for(i=1;i<=n;i++)
		for(j=1;j<=m;j++)
			cin>>a[i][j];
	for(i=1;i<=n;i++)
		for(j=1;j<=m;j++)
		{
			if(a[i][j])
				dp[i][j]=min(dp[i-1][j-1],min(row[i][j-1],col[i-1][j]))+1;
			else
			{
				row[i][j]=row[i][j-1]+1;
				col[i][j]=col[i-1][j]+1;
			}
			res=max(res,dp[i][j]);
		}
	cl(dp,0);
	cl(row,0);
	cl(col,0);
	for(i=1;i<=n;i++)
		for(j=m;j>=1;j--)
		{
			if(a[i][j])
				dp[i][j]=min(dp[i-1][j+1],min(row[i][j+1],col[i-1][j]))+1;
			else
			{
				row[i][j]=row[i][j+1]+1;
				col[i][j]=col[i-1][j]+1;
			}
			res=max(res,dp[i][j]);
		}
	cout<<res<<endl;
	return 0;
}