Description
农夫约翰想要在他的正方形农场上建造一座正方形大牛棚。他讨厌在他的农场中砍树,想找一个能够让他在空旷无树的地方修建牛棚的地方。我们假定,他的农场划分成 N x N 的方格。输入数据中包括有树的方格的列表。你的任务是计算并输出,在他的农场中,不需要砍树却能够修建的最大正方形牛棚。牛棚的边必须和水平轴或者垂直轴平行。
格式
EXAMPLE
考虑下面的方格,它表示农夫约翰的农场,‘.’表示没有树的方格,‘#’表示有树的方格
1 2 3 4 5 6 7 8
1 … … . .
2 . # … # . .
3 … … . .
4 … … . .
5 … … . .
6 . . # … . .
7 … … . .
8 … … . .
最大的牛棚是 5 x 5 的,可以建造在方格右下角的两个位置其中一个。
Input
Line 1: 两个整数: N (1 <= N <= 1000),农场的大小,和 T (1 <= T <= 10,000)有树的方格的数量
Lines 2..T+1: 两个整数(1 <= 整数 <= N), 有树格子的横纵坐标
Output
输出文件只由一行组成,约翰的牛棚的最大边长
Analysis
两种做法
- 正解dp, f[i][j]=min(f[i−1][j],f[i][j−1],f[i−1][j−1])+1 ,其中 f[i][j] 表示以 i,j 为右下角的正方形的最大长度, O(n2)
- 二分答案暴力一下, O(n2logn2) 的还是复杂度可以
usaco是没有题目了吗为什么牛棚梗可以玩这么久从1.0一直到5.0
Code
dp解法
/*
ID:wjp13241
PROG:range
LANG:C++
*/
#include <stdio.h>
using namespace std;
int f[1001][1001];
bool map[1001][1001];
int max(int x,int y)
{
return x>y?x:y;
}
int min(int x,int y)
{
return x<y?x:y;
}
int main()
{
int n,m;
scanf("%d%d",&n,&m);
for (int i=1;i<=m;i++)
{
int x,y;
scanf("%d%d",&x,&y);
map[x][y]=true;
}
int ans=0;
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
if (!map[i][j])
{
f[i][j]=min(min(f[i-1][j],f[i][j-1]),f[i-1][j-1])+1;
ans=max(f[i][j],ans);
}
printf("%d\n",ans);
return 0;
}
二分解法
#include <stdio.h>
#include <cmath>
#include <iostream>
using namespace std;
int x[1011001],y[1011001];
int map[1111][1111];
int f[1111][1111];
int acre(int i,int j,int k)
{
return f[i+k-1][j+k-1]-f[i-1][j+k-1]-f[i+k-1][j-1]+f[i-1][j-1];
}
bool check(int i,int j,int k)
{
return acre(i,j,k)==k*k;
}
int main()
{
int n,m,len,ans=0,count=0;
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
map[i][j]=1;
for (int i=1;i<=m;i++)
{
int x,y;
scanf("%d%d",&x,&y);
map[x][y]=0;
}
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
if (map[i][j])
{
x[++count]=i;
y[count]=j;
}
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
f[i][j]=f[i-1][j]+f[i][j-1]-f[i-1][j-1]+map[i][j];
int l=1,r=n;
while (l<=r)
{
int mid=(l+r)>>1;
int k=mid;
bool flag=false;
for (int v=1;v<=count;v++)
{
if ((x[v]+k-1<=n)&&(y[v]+k-1<=n)&&check(x[v],y[v],k))
{
flag=true;
if (acre(x[v],y[v],k)>ans)
{
ans=acre(x[v],y[v],k);
len=k;
}
break;
}
}
if (flag)
l=mid+1;
else
r=mid-1;
}
printf("%d\n",len);
return 0;
}