[NOIP2010] 引水入城

★★★   输入文件:flow.in   输出文件:flow.out   简单对比
时间限制:1 s   内存限制:128 MB

在一个遥远的国度,一侧是风景秀美的湖泊,另一侧则是漫无边际的沙漠。该国的行政区划十分特殊,刚好构成一个N行M列的矩形,如上图所示,其中每个格子都代表一座城市,每座城市都有一个海拔高度。

为了使居民们都尽可能饮用到清澈的湖水,现在要在某些城市建造水利设施。水利设施有两种,分别为蓄水厂和输水站。蓄水厂的功能是利用水泵将湖泊中的水抽取到所在城市的蓄水池中。因此,只有与湖泊毗邻的第1行的城市可以建造蓄水厂。而输水站的功能则是通过输水管线利用高度落差,将湖水从高处向低处输送。故一座城市能建造输水站的前提,是存在比它海拔更高且拥有公共边的相邻城市,已经建有水利设施。

由于第N行的城市靠近沙漠,是该国的干旱区,所以要求其中的每座城市都建有水利设施。那么,这个要求能否满足呢?如果能,请计算最少建造几个蓄水厂;如果不能,求干旱区中不可能建有水利设施的城市数目。

【输入】

输入文件名为flow.in。输入文件的每行中两个数之间用一个空格隔开。

输入的第一行是两个正整数N和M,表示矩形的规模。

接下来N行,每行M个正整数,依次代表每座城市的海拔高度。

【输出】

输出文件名为flow.out。

输出有两行。如果能满足要求,输出的第一行是整数1,第二行是一个整数,代表最少建造几个蓄水厂;如果不能满足要求,输出的第一行是整数0,第二行是一个整数,代表有几座干旱区中的城市不可能建有水利设施。

【输入输出样例1】

flow.in

2 5

9 1 5 4 3

8 7 6 1 2

flow.out

1

1

【样例1说明】

只需要在海拔为9的那座城市中建造蓄水厂,即可满足要求。

【输入输出样例2】

flow.in

3 6

8 4 5 6 4 4

7 3 4 3 3 3

3 2 2 1 1 2

flow.out

1

3

【样例2说明】

湖泊

8 4 5 6 4 4

7 3 4 3 3 3

3 2 2 1 1 2

沙漠

上图中,在3个粗线框出的城市中建造蓄水厂,可以满足要求。以这3个蓄水厂为源头在干旱区中建造的输水站分别用3种颜色标出。当然,建造方法可能不唯一。

【数据范围】

本题共有10个测试数据,每个数据的范围如下表所示:

测试数据编号能否满足要求N M

1不能 N≤10 M ≤ 10

2不能 N≤100M≤ 100

3不能 N≤500 M≤ 500

4能 N= 1 M≤ 10

5能 N≤10 M ≤ 10

6能 N≤100 M≤ 20

7能 N≤100 M≤ 50

8能 N≤100 M≤100

9能 N≤200 M≤ 200

10能N≤500 M≤ 500

对于所有的10个数据,每座城市的海拔高度都不超过10^6。

思路

bfs||dfs+线段覆盖

仔细想一下可以发现,每一个临海点对应的沙漠点集一定为一条线段,然后就可以用线段覆盖的贪心解决。

代码实现

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 using namespace std;
 5 inline int max_(int x,int y){return x>y?x:y;}
 6 const int maxn=1e3;
 7 int n,m,ok,now,ans;
 8 int a,b,head,tail;
 9 int map[maxn][maxn];
10 struct nate{int l,r;}s[maxn];
11 bool v[maxn][maxn],d[maxn],f[maxn];
12 int ox,oy,nx,ny;
13 int bx[]={1,0,-1,0};
14 int by[]={0,1,0,-1};
15 bool comp(const nate&x,const nate&y){return x.l<y.l;}
16 struct nata{int x,y;}q[maxn*maxn];
17 void bfs(int x,int y){
18     memset(v,0,sizeof(v));
19     head=tail=0;
20     q[head++]=(nata){x,y};
21     v[x][y]=1;
22     while(head>tail){
23         ox=q[tail].x;
24         oy=q[tail].y;
25         tail++;
26         for(int i=0;i<4;i++){
27             nx=ox+bx[i];
28             ny=oy+by[i];
29             if(!v[nx][ny]&&map[nx][ny]<map[ox][oy]){
30                 q[head++]=(nata){nx,ny};
31                 v[nx][ny]=1;
32             }
33         }
34     }
35 }
36 int main(){
37        freopen("flow.in","r",stdin);
38        freopen("flow.out","w",stdout);
39     memset(map,0x7f,sizeof(map));
40        scanf("%d%d",&n,&m);
41        for(int i=1;i<=n;i++)
42        for(int j=1;j<=m;j++)
43        scanf("%d",&map[i][j]);
44        for(int i=1;i<=m;i++){
45            if(!d[i]){
46                bfs(1,i);
47                for(int j=1;j<=m;j++){
48                    if(v[1][j]) d[j]=1;
49                    if(!s[now].l&&v[n][j]) s[now].l=j;
50                    if(v[n][j]) s[now].r=j,f[j]=1;
51             }
52             if(s[now].l) now++;
53            }
54        }
55        for(int i=1;i<=m;i++) if(f[i]) ++ok;
56        if(ok==m){
57            sort(s,s+now,comp);
58            for(int i=0;i<now;i++){
59                if(s[i].l<=a||s[i].r<=b) b=max_(s[i].r,b);
60                else a=max_(s[i].l,b+1),b=s[i].r,ans++;
61                
62         }
63            printf("1
%d
",ans);
64        }
65        else printf("0
%d
",m-ok);
66     return 0;
67 }

代码极丑,求轻喷。

贪心挂了好多好多次,我真的弱的不行。。。

原文地址:https://www.cnblogs.com/J-william/p/6880392.html