BZOJ1619: [Usaco2008 Nov]Guarding the Farm 保卫牧场（BFS）

2016-04-03 · 5 min read · bfs dfs

1619: [Usaco2008 Nov]Guarding the Farm 保卫牧场

Time Limit: 5 Sec Memory Limit: 64 MB Submit: 661 Solved: 292 [Submit][Status][Discuss]

Description

The farm has many hills upon which Farmer John would like to place guards to ensure the safety of his valuable milk-cows. He wonders how many guards he will need if he wishes to put one on top of each hill. He has a map supplied as a matrix of integers; the matrix has N (1 < N <= 700) rows and M (1 < M <= 700) columns. Each member of the matrix is an altitude H_ij (0 <= H_ij <= 10,000). Help him determine the number of hilltops on the map. A hilltop is one or more adjacent matrix elements of the same value surrounded exclusively by either the edge of the map or elements with a lower (smaller) altitude. Two different elements are adjacent if the magnitude of difference in their X coordinates is no greater than 1 and the magnitude of differences in their Y coordinates is also no greater than 1.

农夫JOHN的农夫上有很多小山丘，他想要在那里布置一些保镖（……）去保卫他的那些相当值钱的奶牛们。他想知道如果在一座小山丘上布置一名保镖的话，他总共需要招聘多少名保镖。他现在手头有一个用数字矩阵来表示地形的地图。这个矩阵有N行（1 < N < = 100)和M列( 1 < M < = 70) 。矩阵中的每个元素都有一个值H_ij(0 < = H_ij < =10,000)来表示该地区的海拔高度。请你帮助他统计出地图上到底有多少个小山丘。小山丘的定义是：若地图中一个元素所邻接的所有元素都比这个元素高度要小（或它邻接的是地图的边界），则该元素和其周围所有按照这样顺序排列的元素的集合称为一个小山丘。这里邻接的意义是：若一个元素与另一个横坐标纵坐标和它的横纵坐标相差不超过1，则称这两个元素邻接。问题名称：guard 输入格式：第一行：两个由空格隔开的整数N和M 第二行到第N+1行：第I+1行描述了地图上的第I行，有M个由空格隔开的整数：H_ij.

Input

Line 1: Two space-separated integers: N and M
Lines 2..N+1: Line i+1 describes row i of the matrix with M space-separated integers: H_ij

Output

Line 1: A single integer that specifies the number of hilltops

Sample Input

8 7 4 3 2 2 1 0 1 3 3 3 2 1 0 1 2 2 2 2 1 0 0 2 1 1 1 1 0 0 1 1 0 0 0 1 0 0 0 0 1 1 1 0 0 1 2 2 1 1 0 0 1 1 1 2 1 0

Sample Output

HINT

三个山丘分别是：左上角的高度为4的方格，右上角的高度为1的方格，还有最后一行中高度为2的方格．

Source

Silver

[Submit][Status][Discuss] HOME Back\

题目的中文描述的数据范围有问题！

思路：BFS,我的做法是预处理所有“坏点”，也就是如果点（i,j）的相邻点中有比点（i,j）高的，那么点(i,j)就是坏点。

然后bfs的时候，对于存在一个坏点的联通块（只有数字相同的点能构成一个联通块）都标记成坏点。

然后用左上角开始做bfs...注意遇到坏点跳过但是不要直接标记访问，因为可能之后的点会用到之前的坏点标记自身。

然而因为数据范围有问题wa了好多次。。。

于是看了题解：

大概都是另一种做法（而且都加入了输入挂，不懂，怎么感觉这些人是互相抄的2333）

从最高的点开始做dfs，每次把不大于当前点高度的点都访问掉。。。算作一次访问。

访问次数就是答案。

想了下确实比较巧妙，比较好写，然而比我的写法要慢。

然而窝不懂为什么这种写法数组开小就是返回RE,而我的写法数组开销就是返回WA...

我的做法（800ms,bfs）：

/* ***********************************************
Author :111qqz
Created Time :2016年04月03日 星期日 20时01分12秒
File Name :code/bzoj/1619.cpp
************************************************ */

 1#include <cstdio>
 2#include <cstring>
 3#include <iostream>
 4#include <algorithm>
 5#include <vector>
 6#include <queue>
 7#include <set>
 8#include <map>
 9#include <string>
10#include <cmath>
11#include <cstdlib>
12#include <ctime>
13#define fst first
14#define sec second
15#define lson l,m,rt<<1
16#define rson m+1,r,rt<<1|1
17#define ms(a,x) memset(a,x,sizeof(a))
18typedef long long LL;
19#define pi pair < int ,int >
20#define MP make_pair

 1using namespace std;
 2const double eps = 1E-8;
 3const int dx4[4]={1,0,0,-1};
 4const int dy4[4]={0,-1,1,0};
 5const int dx8[8]={1,1,1,0,0,-1,-1,-1};
 6const int dy8[8]={1,0,-1,-1,1,1,0,-1};
 7const int inf = 0x3f3f3f3f;
 8const int N=705;
 9int h[N][N];
10bool bad[N][N];
11int n,m;
12bool vis[N][N];
13void checkvis();

1struct node
2{
3    int x,y;

1    bool inmaze ()
2    {
3	if (x>=1&&y>=1&&x<=n&&y<=m)  return true;
4	return false;
5    }

 1    void look()
 2    {
 3	cout<<x<<" "<<y<<endl;
 4    }
 5}s;
 6bool ok ( int x,int y)
 7{
 8 //   if (x-1>=1&&h[x-1][y]>h[x][y])  return false;
 9 //   if (x+1<=n&&h[x+1][y]>h[x][y])  return false;
10 //   if (y-1>=1&&h[x][y-1]>h[x][y])  return false;
11 //   if (y+1<=m&&h[x][y+1]>h[x][y])  return false;

1    for ( int i = 0 ; i < 8 ; i++)
2    {
3	int nx = x + dx8[i];
4	int ny = y + dy8[i];
5	if (nx<1||nx>n||ny<1||ny>m) continue;
6	if (h[nx][ny]>h[x][y]) return false;
7    }
8    return true;
9}

 1bool bfs( int x,int y)
 2{
 3    queue<node>q;
 4    s.x = x;
 5    s.y = y;
 6    q.push(s);
 7    vis[x][y] = true;
 8    bool res = true;
 9    while (!q.empty())
10    {
11	node pre = q.front() ;q.pop();
12//	pre.look();
13//	checkvis();
14	if (!res) bad[pre.x][pre.y] = true;

 1	for ( int i = 0 ; i < 8 ; i++)
 2	{
 3	    node nxt;
 4	    nxt.x = pre.x + dx8[i];
 5	    nxt.y = pre.y + dy8[i];
 6	    if (!nxt.inmaze()) continue;
 7	    if (vis[nxt.x][nxt.y]) continue;
 8	  //  vis[nxt.x][nxt.y] = true;
 9	    if (h[nxt.x][nxt.y]!=h[pre.x][pre.y]) continue;
10	    if (bad[nxt.x][nxt.y])
11	    {
12		res = false;
13	    }
14	    if (!ok(nxt.x,nxt.y))
15	    {
16		res = false;
17	    }
18//	    nxt.look();
19	    q.push(nxt);
20	    vis[nxt.x][nxt.y] = true;
21	}
22    }

    return res;

}

 1void checkbad()
 2{
 3    for ( int i = 1 ; i <= n ; i++)
 4    {
 5	for ( int j = 1 ; j <= m ; j++)
 6	{
 7	    cout<<bad[i][j]<<" ";
 8	}
 9	cout<<endl;
10    }
11}

 1void checkvis()
 2{
 3    for ( int i = 1 ; i <= n ; i++)
 4    {
 5	for ( int j = 1 ; j <= m ; j++)
 6	{
 7	    cout<<vis[i][j]<<" "; 
 8	}
 9	cout<<endl;
10    }
11}
12int main()
13{
14	#ifndef  ONLINE_JUDGE 
15	freopen("code/in.txt","r",stdin);
16  #endif

1	ios::sync_with_stdio(false);
2	cin>>n>>m;
3	for ( int i = 1 ; i <= n ; i++)
4	    for ( int j = 1 ; j <= m ; j++) cin>>h[i][j];
5	ms(bad,false);

 1	for ( int i = 1 ; i <= n ; i++)
 2	{
 3	    for ( int j = 1 ; j <= m ; j++)
 4	    {
 5		if (!ok(i,j)) bad[i][j] = true;
 6	    }
 7	}
 8	ms(vis,false);
 9//	checkbad();
10	int ans = 0 ;
11	for ( int i = 1 ; i <= n ; i++)
12	{
13	    for ( int j = 1 ; j  <= m ; j++)
14	    {
15		if (vis[i][j]) continue;
16		if (bad[i][j]) continue;
17		if (bfs(i,j))
18		{
19		    ans++;
20		  //  cout<<i<<" "<<j<<endl;
21		}
22	//	else cout<<"bad:"<<endl;
23	    }
24	}

	cout<<ans<<endl;

1  #ifndef ONLINE_JUDGE  
2  fclose(stdin);
3  #endif
4    return 0;
5}

网上题解的做法（940ms,dfs）：

/* ***********************************************
Author :111qqz
Created Time :2016年04月03日 星期日 21时33分17秒
File Name :code/bzoj/r1619.cpp
************************************************ */

 1#include <cstdio>
 2#include <cstring>
 3#include <iostream>
 4#include <algorithm>
 5#include <vector>
 6#include <queue>
 7#include <set>
 8#include <map>
 9#include <string>
10#include <cmath>
11#include <cstdlib>
12#include <ctime>
13#include <assert.h>
14#define fst first
15#define sec second
16#define lson l,m,rt<<1
17#define rson m+1,r,rt<<1|1
18#define ms(a,x) memset(a,x,sizeof(a))
19typedef long long LL;
20#define pi pair < int ,int >
21#define MP make_pair

 1using namespace std;
 2const double eps = 1E-8;
 3const int dx4[4]={1,0,0,-1};
 4const int dy4[4]={0,-1,1,0};
 5const int dx8[8]={1,1,1,0,0,-1,-1,-1};
 6const int dy8[8]={-1,0,1,1,-1,-1,0,1};
 7const int inf = 0x3f3f3f3f;
 8const int N=705;
 9int n,m;
10bool vis[N][N];
11int he[N][N];

1struct node
2{
3    int x,y;
4    int val;

1    bool operator < (node b)const
2    {
3	return val>b.val;
4    }
5}h[N*N];

 1inline int read()
 2{
 3    int x=0,f=1;char ch=getchar();
 4    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
 5    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
 6    return x*f;
 7}
 8void dfs ( int x,int y)
 9{
10    vis[x][y] = true;
11    for ( int i = 0 ; i < 8 ; i++)
12    {
13	int nx = x + dx8[i];
14	int ny = y + dy8[i];

 1	if (nx<1||ny<1||nx>n||ny>m) continue;
 2	if (he[x][y]<he[nx][ny]) continue;
 3	if (vis[nx][ny]) continue;
 4	dfs(nx,ny);
 5    }
 6}
 7int main()
 8{
 9	#ifndef  ONLINE_JUDGE 
10	freopen("code/in.txt","r",stdin);
11  #endif

 1	cin>>n>>m;
 2	int cnt = 0 ;
 3	for ( int i = 1 ; i <= n ; i++)
 4	    for ( int j = 1 ; j <= m ; j++)
 5	    {	
 6		    cin>>he[i][j];
 7		    cnt++;
 8		    h[cnt].x = i;
 9		    h[cnt].y = j;
10		    h[cnt].val = he[i][j];
11	    }

	sort(h+1,h+cnt+1);
	ms(vis,false);

 1	int ans = 0 ;
 2	for ( int i = 1 ; i <= cnt ; i++)
 3	{
 4	    int x = h[i].x;
 5	    int y = h[i].y;
 6	    assert(x>=1&&x<N);
 7	    assert(y>=1&&y<N);
 8	    if (!vis[x][y])
 9	    {
10		dfs(x,y);
11		ans++;
12	    }

	}

	cout<<ans<<endl;

1  #ifndef ONLINE_JUDGE  
2  fclose(stdin);
3  #endif
4    return 0;
5}

BZOJ1619: [Usaco2008 Nov]Guarding the Farm 保卫牧场 （BFS）