[POI2007]POW-The Flood

题目描述

给定一张地势图，所有的点都被水淹没，现在有一些关键点，要求放最少的水泵使所有关键点的水都被抽干

输入输出格式

输入格式：

In the first line of the standard input there are two integers  and , separated by a single space,. The following  lines contain the description of the map. The 'th linedescribes the 'th row ofunitary squares in the map. It contains  integers , separated by single spaces,, .

The number  describes the 'th square of the ![](http://main.edu.pl/images/OI14/pow-en-tex.14.pn…

输出格式：

Your programme should write out one integer to the standard output - the minimum number of pumpsneeded to drain Byteburg.

输入输出样例

输入样例#1：

6 9
-2 -2 -1 -1 -2 -2 -2 -12 -3
-2 1 -1 2 -8 -12 2 -12 -12
-5 3 1 1 -12 4 -6 2 -2
-5 -2 -2 2 -12 -3 4 -3 -1
-5 -6 -2 2 -12 5 6 2 -1
-4 -8 -8 -10 -12 -8 -6 -6 -4

输出样例#1：

2
我们首先考虑如果在格子 a 修建一个抽水机，在什么情况下格子 b 的水也可以被抽干。
我们可以发现当且仅当存在一条从 a 到 b 的路径，中间经过的抽水机（包括 a）的高度都不大于 b 的高度。
即h[b]>=max(h[i])
因此我们可以考虑把所有格子的高度从小到大排序，我们把每一个格子建成一个集合。
然后按照海拔高度从小到大扫描格子，对于当前的格子 i，我们找到所有与 i 相邻并且海拔高度不大于格子 i 的格子，
我们发现如果这些格子中的任意一个洪水要是被解决了，那么格子 i 的洪水也可以被解决，所以我们合并这些格子。
对于当前的格子 i，如果它必须被清理且与它相邻的格子集合中没有任何一个被清理，我们则把这个集合的清理状态标记为真，然后答案加 1。
集合和每个格子是否被清理用并查集来维护就可以了。

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 struct node
 {
     int s,x,y;
 } city[],maze[];
 pair<int,int> set[][];
 int dx[]={,,-,,};
 int dy[]={,,,,-};
 int n,m,a[][],tot,cnt,vis[][],ans;
 bool cmp(node a,node b)
 {
     return (a.s<b.s);
 }
 pair<int,int> find(pair<int,int> x)
 {
     if (set[x.first][x.second]!=x) set[x.first][x.second]=find(set[x.first][x.second]);
     return set[x.first][x.second];
 }
 void union_set(pair<int,int> x,pair<int,int> y)
 {
     x=find(x),y=find(y);
         set[x.first][x.second]=y;
         vis[y.first][y.second]|=vis[x.first][x.second];
 }
 void exam(int x,int y)
 {
     for(int i=;i<=;i++)
     {
         int tx=x+dx[i], ty=y+dy[i];
         if(tx<=||ty<=||tx>m||ty>n) continue;
         if(a[tx][ty]>a[x][y]) continue;
         union_set(make_pair(x, y), make_pair(tx, ty));
     }
 }
 int main()
 {
     int i,j;
     cin>>n>>m;
     for (i=; i<=n; i++)
     {
         for (j=; j<=m; j++)
         {
             scanf("%d",&a[i][j]);
             if (a[i][j]<=)
             {
                 a[i][j]=-a[i][j];
             }
             else
             {
                 city[++tot]=(node){a[i][j],i,j};
             }
             maze[++cnt]=(node){a[i][j],i,j};
             set[i][j]=make_pair(i,j);
         }
     }
     sort(maze+,maze+cnt+,cmp);
     sort(city+,city+tot+,cmp);
     for (i=; i<=tot; i++)
     {
         for (j=; j<=cnt,city[i].s>=maze[j].s; j++)
         {
             exam(maze[j].x,maze[j].y);
             pair<int,int> x=find(make_pair(city[i].x,city[i].y));
             if (vis[x.first][x.second]==)
             {
                 ans++;
                 vis[x.first][x.second]=;
             }
         }
     }
     cout<<ans;
 }