Big Barn_usaco5.3_dp

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

---

两种做法

1. 正解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(n^{2})$
2. 二分答案暴力一下，$O(n^{2}\log_{2}^{n})$的还是复杂度可以

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;
}