Flood Fill算法也叫“泛洪算法”,直观上有种,在宣纸上滴墨水,渐次染开的感觉。在图像处理中常常被用到。
利益相关:题、图来自AcWing
AcWing
1.一句话概括
多重循环遍历整个空间,对满足条件的元素:以该元素为起点,做BFS。
时间复杂度为O(n),常被用来求解与图的连通集相关的题目,
大致有一下应用:
- 连通集的大小
- 连通集的个数
- 多点之间的连通性
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;
}