codeforces#375（div.2）723D - Lakes in Berland dfs+bfs

题意：给你一张n×m的图，点代表水，星代表地，当水连在一块时，称为湖，但是，若湖中有水在边界上，则不算它是湖，现在要求你去掉x个湖，并用地将其填上，使原图剩下k个湖，输出你最少需要多少个单位的地，并且输出你填补过的图。

比赛的时候没过，好菜呀。

思路：首先定义两个数组bfs_vis,dfs_vis,分别记录bfs和dfs的访问状况，定义flag，true代表一块区域属于湖，false代表一块区域不属于湖,vector[top]记录第top个湖内的点的情况，top初始化为0。

然后说下dfs的思路，从一个'.'点出发,并且只走'.',若其走到了边界上，则标记flag为false，此时并不结束dfs，而是继续走，将这一整块水域全部走完，并且dfs_vis都标记为true代表已经走过。相反，若走不到边界上，则一定能走完这个水域，flag为true，在dfs的同时，将走过的每一个点都装入vec[top]。当dfs完成后，若flag为false，则清空vec[top],否则top++，进入下个vec[top]记录新的湖。

然后说bfs，就是从(1,1）开始遍历，若找到一个点'.',并且其dfs_vis为false的话，说明它没有在dfs中被访问过，则对其进行dfs，当dfs完成后，若flag为true，说明这个点所在水域为湖，则记录它的大小和序号。

最后只要对这些湖的大小进行升序排序，并且相应的将其中的点填上就行了。

为什么这么强调dfs_vis呢，因为我就坑在这啦，最开始的时候，我在dfs中一旦发现某个点走到了边界就立即停止dfs并返回false，这样出现的问题就是当一个水域只有一个点在边界的时候，假如我们第一次就dfs到了改点，那么就相当于将其变成了陆地。则其剩下的水域在下一次dfs的时候就会被判定为一个湖。

代码：

#include<bits/stdc++.h>
#define MEM(a,x) memset(a,x,sizeof(a));
#define MEMINF(a) memset(a,0x3f,sizeof(a));
using namespace std;
typedef long long LL;
const int MAXN=4000;
const int INF=0x3f3f3f3f;
const int MOD=1000000007;
int n,m,k;
char mp[MAXN][MAXN];
int dx[]={0,0,1,-1};
int dy[]={1,-1,0,0};
int top;
struct node {
  int x,y;
};
bool bfs_vis[MAXN][MAXN];
bool dfs_vis[MAXN][MAXN];
bool cmp(node a,node b) {
  return a.x<b.x;
}

vector<node>que;
queue<node>q;
vector<node> vec[MAXN];
bool flag;
void dfs(node u) {
  vec[top].push_back(u);
  dfs_vis[u.x][u.y]=1;
  if (u.x==1||u.x==n||u.y==1||u.y==m) flag=false;
  for (int i=0; i<4; ++i) {
    node v=u;
    v.x+=dx[i];
    v.y+=dy[i];
    if (dfs_vis[v.x][v.y])continue;
    if (mp[v.x][v.y]!='.') continue;
    //printf("%d %d\n",v.x,v.y);
    if (v.x==1||v.x==n||v.y==1||v.y==m) flag=false;
    dfs(v);
  }
}

void bfs(){
  node u;
  u.x=1,u.y=1;
  q.push(u);
  bfs_vis[1][1]=1;
  MEM(bfs_vis,0);
  while(!q.empty()) {
    node no=q.front();
    q.pop();
    //printf("m %c %d\n",mp[4][3],!dfs_vis[4][3]);
    if (mp[no.x][no.y]=='.'&&!dfs_vis[no.x][no.y]) {
      flag=true;
      //puts("-----");
      dfs(no);
      if (flag) {
       //printf("no:%d %d\n",no.x,no.y);
        //puts("*****");
        node lal;
        lal.x=vec[top].size();
        lal.y=top;
        que.push_back(lal);
        top++;
      }
      else {
        vec[top].clear();
      }
    }
    for (int i=0; i<4; ++i) {
      node v=no;
      v.x+=dx[i];
      v.y+=dy[i];
      if (bfs_vis[v.x][v.y])continue;
      if (v.x<1||v.x>n||v.y<1||v.y>m) continue;
      q.push(v);
      bfs_vis[v.x][v.y]=1;
    }
  }
}
      
int main() {
  scanf("%d %d %d",&n,&m,&k);
  top=0;
  for (int i=1; i<=n; ++i) {
    scanf("%s",mp[i]+1);
  }

  bfs();
  //printf("MAP%ld\n",que.size());
  int sum=0;
  sort(que.begin(),que.end(),cmp);
  for (auto i=0; i<top-k; ++i) {
    sum+=que[i].x;
    int num=que[i].y;
    for (int it=0;it<vec[num].size(); ++it) {
      int x=vec[num][it].x;
      int y=vec[num][it].y;
     mp[x][y]='*';
    }
  }

  printf("%d\n",sum);
  for (int i=1; i<=n; ++i) 
   puts(mp[i]+1);

}