洛谷 P1736 创意吃鱼法

最新推荐文章于 2020-09-09 12:37:19 发布

Alex Su (*^▽^*)

最新推荐文章于 2020-09-09 12:37:19 发布

阅读量222

点赞数

CC 4.0 BY-SA版权

分类专栏：动态规划

本文链接：https://blog.youkuaiyun.com/qq_34416123/article/details/82528973

动态规划专栏收录该内容

86 篇文章

订阅专栏

本文深入探讨了寻找矩阵中最大正方形的算法实现，通过动态规划的方法，利用两个辅助矩阵来记录非零元素的最长连续长度，从而找到最大的正方形。文章详细介绍了算法流程，并附带了完整的代码实现。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

和最大正方形很像。
记得对角线有两条
还有2500不要开成3000，否则会mle

#include<cstdio>
#include<algorithm>
#include<cstring>
#define REP(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

const int MAXN = 2512;
int f[MAXN][MAXN], x[MAXN][MAXN], y[MAXN][MAXN];
int a[MAXN][MAXN], n, m, ans, now;

int main()
{
	scanf("%d%d", &n, &m);
	_for(i, 1, n)
		_for(j, 1, m)
			scanf("%d", &a[i][j]);
	 
	_for(i, 1, n)
		_for(j, 1, m)	
		{
			if(a[i][j])
			{	
				f[i][j] = min(f[i-1][j-1], min(y[i][j-1], x[i-1][j])) + 1;
				ans = max(ans, f[i][j]);
			}
			else 
			{
				x[i][j] = x[i-1][j] + 1;
				y[i][j] = y[i][j-1] + 1;
			}
	}
	
	memset(f, 0, sizeof(f));
	memset(x, 0, sizeof(x));
	memset(y, 0, sizeof(y));
	
	_for(i, 1, n)
		for(int j = m; j >= 1; j--)	
		{
			if(a[i][j])
			{	
				f[i][j] = min(f[i-1][j+1], min(y[i][j+1], x[i-1][j])) + 1;
				ans = max(ans, f[i][j]);
			}
			else 
			{
				x[i][j] = x[i-1][j] + 1;
				y[i][j] = y[i][j+1] + 1;
			}
	}

	printf("%d\n", ans);
	
	return 0;
}