算法-搜索专题 Flood Fill算法

Flood Fill算法也叫“泛洪算法”，直观上有种，在宣纸上滴墨水，渐次染开的感觉。在图像处理中常常被用到。

利益相关：题、图来自AcWing
AcWing

1.一句话概括

多重循环遍历整个空间，对满足条件的元素：以该元素为起点，做BFS。
时间复杂度为O(n)，常被用来求解与图的连通集相关的题目，
大致有一下应用：

1. 连通集的大小
2. 连通集的个数
3. 多点之间的连通性

2.伪代码

for item in mp[0]:mp[n-1] {
	if(item未被访问过) {
		bfs(横坐标，纵坐标) 
	} 
}

bfs (横坐标，纵坐标) {
	while(队列不空) {
		...
		标记已访问过
		... 
	} 
}

3.举例

3.1

农夫约翰有一片 N∗M 的矩形土地。

最近，由于降雨的原因，部分土地被水淹没了。

现在用一个字符矩阵来表示他的土地。

每个单元格内，如果包含雨水，则用”W”表示，如果不含雨水，则用”.”表示。

现在，约翰想知道他的土地中形成了多少片池塘。

每组相连的积水单元格集合可以看作是一片池塘。

每个单元格视为与其上、下、左、右、左上、右上、左下、右下八个邻近单元格相连。

请你输出共有多少片池塘，即矩阵中共有多少片相连的”W”块。

输入格式
第一行包含两个整数 N 和 M。

接下来 N 行，每行包含 M 个字符，字符为”W”或”.”，用以表示矩形土地的积水状况，字符之间没有空格。

输出格式
输出一个整数，表示池塘数目。

数据范围
1≤N,M≤1000
输入样例：
10 12

输出样例：
3

3.2

图1是一个城堡的地形图。
请你编写一个程序，计算城堡一共有多少房间，最大的房间有多大。
城堡被分割成 m∗n个方格区域，每个方格区域可以有0~4面墙。
注意：墙体厚度忽略不计。

输入格式
第一行包含两个整数 m 和 n，分别表示城堡南北方向的长度和东西方向的长度。接下来 m 行，每行包含 n 个整数，每个整数都表示平面图对应位置的方块的墙的特征。每个方块中墙的特征由数字 P 来描述，我们用1表示西墙，2表示北墙，4表示东墙，8表示南墙，P 为该方块包含墙的数字之和。例如，如果一个方块的 P 为3，则 3 = 1 + 2，该方块包含西墙和北墙。城堡的内墙被计算两次，方块(1,1)的南墙同时也是方块(2,1)的北墙。

输入的数据保证城堡至少有两个房间。

输出格式
共两行，第一行输出房间总数，第二行输出最大房间的面积（方块数）。

数据范围
1≤m,n≤50,
0≤P≤15
输入样例：
4 7
11 6 11 6 3 10 6
7 9 6 13 5 15 5
1 10 12 7 13 7 5
13 11 10 8 10 12 13
输出样例：
5
9

• Flood Fill的简单应用
• 关键在于处理代表西北东南四面墙的四个数字：1，2，4，8。如何处理呢？用二进制表示为：1，10，100，1000即用1所在数位0，1，2，3表示西北东南四个方向有没有墙。存在1则存在墙。

#include<bits/stdc++.h>

using namespace std;

const int N=55;
typedef pair<int,int> PII;

int m,n，cnt,area;
int mp[N][N];
bool is[N][N];
queue<PII> q;

int bfs(int i,int j) {
	int dx[4] = {0,-1,0,1};
	int dy[4] = {-1,0,1,0};
	
	int res = 0 ;
	q.push({i,j});
	while(q.size()) {
		int x = q.front().first;
		int y = q.front().second;
		int z = mp[x][y];
		q.pop();
		for(int i=0;i<4;i++) {
			int k = z>>i&1; // 看z的二进制的第i个数位是不是1
			int xx = x + dx[i];
			int yy = y + dy[i];
			if(!k && !is[xx][yy] && xx >= 1 && xx <= m && yy >= 1 && yy <= n) {
				q.push({xx,yy});
				is[xx][yy] = 1;
				res ++;
			}
 		} 
	}
	return res;
}

int main()
{
	cin >> m >> n;
	for(int i=1;i<=m;i++) 
		for(int j=1;j<=n;j++) 
			cin >> mp[i][j];

	for(int i=1;i<=m;i++) 
		for(int j=1;j<=n;j++)
			if(!is[i][j]) {
				area = max(area,bfs(i,j));
				cnt ++;
			}
	cout << cnt << endl;
	cout << area << endl;
	
	return 0;
}

3.3

FGD小朋友特别喜欢爬山，在爬山的时候他就在研究山峰和山谷。
为了能够对旅程有一个安排，他想知道山峰和山谷的数量。
给定一个地图，为FGD想要旅行的区域，地图被分为 n×n 的网格，每个格子 (i,j) 的高度 w(i,j) 是给定的。

若两个格子有公共顶点，那么它们就是相邻的格子，如与 (i,j) 相邻的格子有(i−1,j−1),(i−1,j),(i−1,j+1),(i,j−1),(i,j+1),(i+1,j−1),(i+1,j),(i+1,j+1)。
我们定义一个格子的集合 S 为山峰（山谷）当且仅当：

S 的所有格子都有相同的高度。S 的所有格子都连通。对于 s 属于 S，与 s 相邻的 s′ 不属于 S，都有 ws>ws′（山峰），或者 ws<ws′（山谷）。如果周围不存在相邻区域，则同时将其视为山峰和山谷。你的任务是，对于给定的地图，求出山峰和山谷的数量，如果所有格子都有相同的高度，那么整个地图即是山峰，又是山谷。

输入格式
第一行包含一个正整数 n，表示地图的大小。

接下来一个 n×n 的矩阵，表示地图上每个格子的高度 w。

输出格式
共一行，包含两个整数，表示山峰和山谷的数量。

数据范围
1≤n≤1000,
0≤w≤10^9
输入样例1：
5
8 8 8 7 7
7 7 8 8 7
7 7 7 7 7
7 8 8 7 8
7 8 8 8 8
输出样例1：
2 1
输入样例2：
5
5 7 8 3 1
5 5 7 6 6
6 6 6 2 8
5 7 2 5 8
7 1 0 1 7
输出样例2：
3 3
样例解释
样例1：

样例2：

#include<bits/stdc++.h>

using namespace std;

typedef pair<int,int> PII;
const int N=1010;
int n,peak,vag;
int mp[N][N];
bool is[N][N];

queue<PII> q;
set<int> s;

int bfs(int i,int j) {
	q.push({i,j});
	bool f = 0,g = 0;
	while(q.size()) {
		int x = q.front().first;
		int y = q.front().second;
		int z = mp[x][y];
		for(int i=x-1;i <= x+1; i++) {
			for(int j=y-1;j <= y+1; j++) {
				if(i == x && j == y) continue;
				if(i<1 || i>n || j<1 || j>n) continue;
				
				if(mp[i][j] == z ) {
					if(!is[i][j]) {
						q.push({i,j});
						is[i][j] = 1;
					}
				}else {
					if(z > mp[i][j]) f = 1;
					if(z < mp[i][j]) g = 1;
				}
			}
		}		
		q.pop();
	}
	if(f && g) return 0;
	if(!f && !g) {
	    peak ++;
	    vag ++;
	    return 0;
	}
	if(f) peak++;
	if(g) vag ++;
	return 0; 
}

int main()
{
    cin.tie(0);
	cin >> n;
	bool ip = 0;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++) 
			cin >> mp[i][j];
	for(int i=1;i<=n;i++) 
		for(int j=1;j<=n;j++) 
			if(!is[i][j]) 
				bfs(i,j);
	cout << peak <<" "<< vag <<endl;
	
	return 0;
}