JDOJ-1239: VIJOS-P1057 盖房子

[文章目录]

Description

永恒の灵魂最近得到了面积为n*m的一大块土地(高兴ING^_^),他想在这块土地上建造一所房子,这个房子必须是正方形的。 但是,这块土地并非十全十美,上面有很多不平坦的地方(也可以叫瑕疵)。这些瑕疵十分恶心,以至于根本不能在上面盖一砖一瓦。 他希望找到一块最大的正方形无瑕疵土地来盖房子。 不过,这并不是什么难题,永恒の灵魂在10分钟内就轻松解决了这个问题。 现在,您也来试试吧。

又是一发神奇的动态规划。

dp[i][j]表示以i,j为右下角的正方形最长的边长。

得到转移方程:

如果map[i][j]==1,dp[i][j]=min(dp[i-1][j],dp[i-1][j-1],dp[i][j-1])+1;

若map[i][j]==0,dp[i][j]=0;进行dp,初值将第一行和第一列设成map

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
int n,m;
int map[1100][1100],dp[1100][1100],ans=-1;
int read()
{
	char ch=getchar();int re=0;
	while(ch<'0'||ch>'9')ch=getchar();
	while(ch>='0'&&ch<='9') re=(re<<3)+(re<<1)+ch-'0',ch=getchar();
	return re;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			map[i][j]=read();
	for(int i=1;i<=n;i++) dp[i][1]=map[i][1],ans=max(ans,dp[i][1]);
	for(int j=1;j<=m;j++) dp[1][j]=map[1][j],ans=max(ans,dp[1][j]);
	for(int i=2;i<=n;i++)
		for(int j=2;j<=m;j++)
		{
			if(map[i][j])
				dp[i][j]=min(min(dp[i-1][j],dp[i-1][j-1]),dp[i][j-1])+1;
			ans=max(dp[i][j],ans);
		}
	printf("%d",ans);
	return 0;
}

 

发表评论

邮箱地址不会被公开。 必填项已用*标注