Codevs 1066 引水入城

题目描述 Description

 

在一个遥远的国度，一侧是风景秀美的湖泊，另一侧则是漫无边际的沙漠。该国的行政 区划十分特殊，刚好构成一个N行M列的矩形，如上图所示，其中每个格子都代表一座城 市，每座城市都有一个海拔高度。 为了使居民们都尽可能饮用到清澈的湖水，现在要在某些城市建造水利设施。水利设施 有两种，分别为蓄水厂和输水站。蓄水厂的功能是利用水泵将湖泊中的水抽取到所在城市的 蓄水池中。因此，只有与湖泊毗邻的第1行的城市可以建造蓄水厂。而输水站的功能则是通 过输水管线利用高度落差，将湖水从高处向低处输送。故一座城市能建造输水站的前提，是 存在比它海拔更高且拥有公共边的相邻城市，已经建有水利设施。 由于第N行的城市靠近沙漠，是该国的干旱区，所以要求其中的每座城市都建有水利 设施。那么，这个要求能否满足呢？如果能，请计算最少建造几个蓄水厂；如果不能，求干 旱区中不可能建有水利设施的城市数目。

输入描述 Input Description

输入的每行中两个数之间用一个空格隔开。 输入的第一行是两个正整数N和M，表示矩形的规模。 接下来N行，每行M个正整数，依次代表每座城市的海拔高度。

输出描述 Output Description

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

这道题六个月之前第一次听讲，只得了30，六个月后写，还是30...寻找原因，原来是自己不会写线段覆盖...

第一步呢，我的想法是，对于每一个河边的点进行DFS，求出其能控制的沙漠城市的范围，用一个visall[i]来反过来表示每一个沙漠城市i能否被浇灌，只要某一次dfs浇灌到了城市i，就将visall[i]改为1。

第二步，判断能否满足要求，即是否有visall=0，如果有即不能满足要求，输出visall=0的个数。第三步，就会开始忧虑了啊，每一个储水站可能控制的区域不是连续的，那该怎么办呢？其实，就是连续的一段。如果不连续，那么则会与另一个储水站的储水路线产生交点，则另一个储水站也可以控制那不连续的一段，最后可看做每一个储水站控制了一段连续的区间。

第四部就是做线段覆盖了，求取最少多少线段能够覆盖一个连续区间。我最开始的思路是，当现在选取的右端点能覆盖到现在考虑线段的左端而覆盖不到下一个线段的左端，就选取现在考虑的这一条线段，其实会少很多情况。例如有三个线段[l1,r1],[l2,r2],[l3,r3] , l1<l2<l3; r1>r2; l1,l2能被这时选取的右端覆盖到，l3却覆盖不到。这时我们会选取r2来作为新的右端点，但是实际上选取r1更长，一般会使得答案更好，所以这种贪心方法是不对的。因此我们需要存下来当前右端点能覆盖到的所有左端点的右端点，当答案需要更新时，选出最大的右端点，才为正确的线段覆盖问题解法。

下附AC代码。

#include<iostream>
#include<string.h>
#include<queue>
#include<algorithm>
#define maxn 505
using namespace std;
int height[maxn][maxn];
int cont[maxn][maxn];
int visall[maxn];
int vis[maxn][maxn];
int m,n;
int fx[maxn]={0,-1,1,0,0};
int fy[maxn]={0,0,0,-1,1};
int times=0;
struct line
{
	int x,y;
}c[maxn];
void dfs(int fa,int x,int y)
{
	if(x==m)
	{
		cont[fa][y]=1;
		visall[y]=1;
	}
	vis[x][y]=1;
	for(int i=1;i<=4;i++)
	if(vis[x+fx[i]][y+fy[i]]==0 && height[x+fx[i]][y+fy[i]]<height[x][y])
	{
		if(1<=x+fx[i] && x+fx[i]<=m &&1<=y+fy[i] && y+fy[i]<=n)
		{
			dfs(fa,x+fx[i],y+fy[i]);
		}
	}
	return;
}
bool cmp(line a,line b)
{
	if(a.x!=b.x)
	return a.x<b.x;
	return a.y>b.y;
}
priority_queue<int> ans;

int main()
{
	cin>>m>>n;
	for(int i=1;i<=m;i++)
	for(int j=1;j<=n;j++)
	{
		cin>>height[i][j];
	}
	for(int i=1;i<=n;i++)
	{
		memset(vis,0,sizeof(vis));
		dfs(i,1,i);
	}
	int cnt=0;
	for(int i=1;i<=n;i++)
	{
		if(visall[i]==0)
		{
			cnt++;
		}
	}
	
	if(cnt)
	{
		cout<<0<<endl;
		cout<<cnt<<endl;
		return 0;
	}
	
	for(int i=1;i<=n;i++)
	{
		bool flag=false;
		for(int j=1;j<=n;j++)
		{
			if(cont[i][j]==1 && flag==false)
			{
				c[i].x=j;
				flag=true;
			}
			if(cont[i][j]==1 && flag==true)
			{
				c[i].y=j;
			}
			if(cont[i][j]==0 && flag==true)
			{
				break;
			}
		}
	}	
	sort(c+1,c+1+n,cmp);
	int ans=0;
	int pos=1;	
	while (c[pos].y==0)	
	pos++;	
	int maxr=0;	
	int j=0;
	while (pos<=n)
	{	
		ans++;		
		int i=pos+1;		
		j=0;		
		maxr=c[pos].y;		
		while (c[i].x<=c[pos].y+1&&i<=n)
		{		
			if (c[i].y>maxr) 
			{
				maxr=c[i].y;
				j=i;
			}			
			i++;		
		}		
		if (j) 
		pos=j;		
		else pos=i;	
	}
	cout<<"1"<<endl;
	cout<<ans<<endl;
}