标题：全球变暖
你有一张某海域NxN像素的照片，".“表示海洋、”#"表示陆地，如下所示：
.......
.##....
.##....
....##.
..####.
...###.
.......
其中"上下左右"四个方向上连在一起的一片陆地组成一座岛屿。例如上图就有2座岛屿。
由于全球变暖导致了海面上升，科学家预测未来几十年，岛屿边缘一个像素的范围会被海水淹没。具体来说如果一块陆地像素与海洋相邻(上下左右四个相邻像素中有海洋)，它就会被淹没。
例如上图中的海域未来会变成如下样子：
.......
.......
.......
.......
....#..
.......
.......
请你计算：依照科学家的预测，照片中有多少岛屿会被完全淹没。
【输入格式】
第一行包含一个整数N。 (1 <= N <= 1000)
以下N行N列代表一张海域照片。
照片保证第1行、第1列、第N行、第N列的像素都是海洋。
【输出格式】
一个整数表示答案。
【输入样例】
7
.......
.##....
.##....
....##.
..####.
...###.
.......
【输出样例】
1
资源约定：
峰值内存消耗（含虚拟机） < 256M
CPU消耗 < 1000ms
请严格按要求输出，不要画蛇添足地打印类似：“请您输入…” 的多余内容。
注意：
main函数需要返回0;
只使用ANSI C/ANSI C++ 标准;
不要调用依赖于编译环境或操作系统的特殊函数。
所有依赖的函数必须明确地在源文件中 #include
不能通过工程设置而省略常用头文件。
提交程序时，注意选择所期望的语言类型和编译器类型。
//深度优先搜索找出所有的没有淹没之前的岛屿数量count 更新岛屿，计算出淹没之后的岛屿数量num
//使用了标记数组
// if(count > num) 输出count - num
//else 输出0
//遍历了两次，但是有一个bug就是如果正好有一块岛屿被淹没，又出现了另外一块岛屿被分成了两座岛屿，
//实际上应该输出1，但是程序输出结果为0；测试数据为第三组。
第一组测试数据
/*
7
.......
.##....
.##....
....##.
..####.
...###.
.......
*/
第二组测试数据
/*
8
.#.....#
####.###
########
###..###
###....#
.##..###
.#...#..
........
*/
第三组测试数据
/*
8
.#.....#
...#.###
########
###..###
###....#
.##..###
.#...#..
........
*/
#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int maxn = 1000+10;
int id[4][2] = {1, 0, 0, 1, -1, 0, 0, -1};
int idx[maxn][maxn];
int idx1[maxn][maxn];
char map[maxn][maxn];
int n;
int count = 0; //统计淹没之前的岛屿数量
int num = 0; //统计淹没之后还存在的岛屿数量
bool check(int x, int y){//为统计淹没之后的岛屿做准备
for(int i = 0; i < 4; i++){
int nx = x + id[i][0];
int ny = y + id[i][1];
if(map[nx][ny] != '#' || nx < 0 || nx >= n || ny < 0 || ny >= n){
return false;
}
}
return true;
}
void dfs(int x, int y){//统计更新之前的岛屿数量
if(x >= 0 && x < n && y >= 0 && y < n && map[x][y] == '#' && idx[x][y] == 0){
idx[x][y] = 1;
if(check(x, y)) idx[x][y] = 2;
for(int i = 0; i < 4; i++){
int nx = x + id[i][0];
int ny = y + id[i][1];
dfs(nx, ny);
}
}else{
return;
}
}
void dfs1(int x, int y){//统计更新之后的岛屿数量
if(x >= 0 && x < n && y >= 0 && y < n && idx[x][y] == 2 && idx1[x][y] == 0){
idx1[x][y] = 1;
for(int i = 0; i < 4; i++){
int nx = x + id[i][0];
int ny = y + id[i][1];
dfs1(nx, ny);
}
}else{
return;
}
}
int main(){
memset(idx, 0, sizeof(idx));
memset(idx1, 0, sizeof(idx1));
scanf("%d", &n);
for(int i = 0; i < n; i++){
cin >> map[i];
}
for(int i = 0; i < n; i++){
for(int j = 0; j < n; j++){
if(map[i][j] == '#' && idx[i][j] == 0){
count++;
dfs(i,j);
}
}
}
for(int i = 0;i < n; i++){
for(int j = 0; j < n; j++){
if(idx[i][j] == 2 && idx1[i][j] == 0){
num++;
dfs1(i,j);
}
}
}
if(count > num) cout << count-num << endl;
else cout << 0 << endl;
return 0;
}