USACO5.3.4Big Barn

最新推荐文章于 2023-04-17 16:02:58 发布

原创最新推荐文章于 2023-04-17 16:02:58 发布 · 510 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#DP #区间前缀和

本博客探讨了如何在包含树木分布的农场中，寻找适合建造最大正方形牛棚的位置，通过分析农场布局，利用二维区间前缀和算法或动态规划方法找到最优解。

题目大意：

农夫约翰想要在他的正方形农场上建造一座正方形大牛棚。他讨厌在他的农场中砍树，想找一个能够让他在空旷无树的地方修建牛棚的地方。我们假定，他的农场划分成 N x N 的方格。输入数据中包括有树的方格的列表。你的任务是计算并输出，在他的农场中，不需要砍树却能够修建的最大正方形牛棚。牛棚的边必须和水平轴或者垂直轴平行。 EXAMPLE 考虑下面的方格，它表示农夫约翰的农场，‘.'表示没有树的方格，‘#'表示有树的方格 1 2 3 4 5 6 7 8 1 . . . . . . . . 2 . # . . . # . . 3 . . . . . . . . 4 . . . . . . . . 5 . . . . . . . . 6 . . # . . . . . 7 . . . . . . . . 8 . . . . . . . . 最大的牛棚是 5 x 5 的，可以建造在方格右下角的两个位置其中一个。

输入：

第一行: 两个整数： N （1 <= N <= 1000），农场的大小，和 T （1 <= T <= 10，000）有树的方格的数量

第2..T+1行: 两个整数（1 <= 整数 <= N), 有树格子的横纵坐标

输出：

只由一行组成，约翰的牛棚的最大边长。

样例输入：

8 3
2 2
2 6
6 3

样例输出：

分析：

本来是一个简单的DP，然而我并没有看出来，在考试的时候乱搞了一个二维区间前缀和，居然过了……下面是我二维区间前缀和的代码：

#include<cstdio>
#define MAXN 1000
using namespace std;
inline int min(int a,int b){return a<b?a:b;}
inline int max(int a,int b){return a>b?a:b;}
int n,t,k;
int maxl=1;
int sum[MAXN+5][MAXN+5];
int main()
{
    int i,j,x,y;
    scanf("%d%d",&n,&t);
    for(i=1;i<=t;i++)
    {
        scanf("%d%d",&x,&y);
        sum[x][y]=1;
    }
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
            sum[i][j]+=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1];
    for(i=2;i<=n;i++)
        for(j=2;j<=n;j++)
        {
            for(k=maxl;k<min(i,j);k++)
                if(sum[i][j]-sum[i-k-1][j]-sum[i][j-k-1]+sum[i-k-1][j-k-1]>0)
                    break;
            maxl=max(k,maxl);
        }
    printf("%d\n",maxl);
}

下面的是DP正解，状态转移方程为f[i][j]=min(f[i-1][j],f[i][j-1],f[i-1][j-1])+1。f[i][j]表示以(i,j)为右下角时能达到的最大边长，比我的区间前缀和简单多了：

#include<cstdio>
#define MAXN 1000
using namespace std;
inline int min(int a,int b){return a<b?a:b;}
inline int max(int a,int b){return a>b?a:b;}
int n,t,maxl;
bool field[MAXN+5][MAXN+5];
int f[MAXN+5][MAXN+5];
int main()
{
	int i,j,x,y;
	scanf("%d%d",&n,&t);
	for(i=1;i<=t;i++)
	{
		scanf("%d%d",&x,&y);
		field[x][y]=true;
	}
	for(i=1;i<=n;i++)
		for(j=1;j<=n;j++)
		{
			if(field[i][j])
				f[i][j]=0;
			else
			{
				f[i][j]=min(f[i-1][j],f[i][j-1]);
				f[i][j]=min(f[i][j],f[i-1][j-1])+1;
			}
			maxl=max(f[i][j],maxl);
		}
	printf("%d\n",maxl);
}