[DP]巨大的牛棚

目录

 

1.先看题

题目描述

输入

输出

样例输入

样例输出

2.解决

方法分析

方法落实

代码实现

---

1.先看题

题目描述

时间限制: 1 Sec  内存限制: 64 MB

农夫约翰想要在他的正方形农场上建造一座正方形大牛棚。他讨厌在他的农场中砍树，想找一个能够让他在空旷无树的地方修建牛棚的地方。我们假定，他的农场划分成 N x N 的方格。输入数据中包括有树的方格的列表。你的任务是计算并输出，在他的农场中，不需要砍树却能够修建的最大正方形牛棚。牛棚的边必须和水平轴或者垂直轴平行。

考虑下面的方格，它表示农夫约翰的农场，‘.'表示没有树的方格，‘#'表示有树的方格

  1 2 3 4 5 6 7 8 
1 . . . . . . . . 
2 . # . . . # . . 
3 . . . . . . . . 
4 . . . . . . . . 
5 . . . . . . . . 
6 . . # . . . . . 
7 . . . . . . . . 
8 . . . . . . . .

最大的牛棚是 5 x 5 的，可以建造在方格右下角的两个位置其中一个。

输入

 Line 1: 两个整数： N （1 <= N <= 1000），农场的大小，和 T （1 <= T <= 10，000）有树的方格的数量 Lines 2..T+1: 两个整数（1 <= 整数 <= N), 有树格子的横纵坐标

输出

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

样例输入

8 3
2 2
2 6
6 3

样例输出

5

---

2.解决

方法分析

这道题的数据范围肯定是不能暴力搜索的。

因为农场最大为 1000*1000 的矩阵，所以把每个元素都遍历一遍的时间复杂度为 o(n*n)，最多会遍历 1000000 次，在这个时间上再增加的话就会有点慢了，这时候就该用DP。

方法落实

定义结构体数组 (itn)f，每个结构体包含 (int)x和(int)y，f[i][j].x表示以 i 行 j 列为矩阵右下角的最大无树矩阵的上下边长，y 则为左右边长。为什么要用结构体？因为这样可以提供两个有用数据：上下边长x及左右边长y。

结构体定义如下

struct itn{
    int x,y;
    itn(){}
    itn(int X,int Y){ //构造函数
        x=X,y=Y;
    }
}f[1005][1005];

状态转移方程为：f[i][j]=itn(min(f[i-1][j].x,f[i-1][j-1].x)+1,min(f[i][j-1].y,f[i-1][j-1].y)+1)

举个例子：

空白表示有树

这样算下来，f[3][3] 就等于 itn(3,2)，正好是题目所求。那么中间为什么要判断 f[i-1][j-1]呢？

比如

如果算f[i][j]=itn(f[i-1][j].x+1,f[i][j-1].y+1)，算粗来就是itn(3,3)，3*3的矩阵。而加了中间的判断，得出来是itn(1,1)。

用 k 来存最大边长，即 k=max(k,min(f[i][j].x,f[i][j].y)) 就是答案。

但是，这个方法只适用于求正方形矩阵。

代码实现

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
struct itn{
    int x,y;
    itn(){}
    itn(int X,int Y){
        x=X,y=Y;
    }
}f[1005][1005];
int n,m,i,j,k,c,d;
bool a[1005][1005];
int main()
{
    scanf("%d%d",&n,&m);
    for(i=1;i<=m;i++)scanf("%d%d",&c,&d),a[c][d]=1;
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++){
            if(!a[i][j])
                f[i][j]=itn(min(f[i-1][j].x,f[i-1][j-1].x)+1,min(f[i][j-1].y,f[i-1][j-1].y)+1);
            k=max(k,min(f[i][j].x,f[i][j].y));
        } //只遍历一次
    printf("%d\n",k);
    return 0;
}