[洛谷]P1514 引水入城 (#搜索 线性dp)

题目描述

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

为了使居民们都尽可能饮用到清澈的湖水，现在要在某些城市建造水利设施。水利设施有两种，分别为蓄水厂和输水站。蓄水厂的功能是利用水泵将湖泊中的水抽取到所在城市的蓄水池中。

因此，只有与湖泊毗邻的第11 行的城市可以建造蓄水厂。而输水站的功能则是通过输水管线利用高度落差，将湖水从高处向低处输送。故一座城市能建造输水站的前提，是存在比它海拔更高且拥有公共边的相邻城市，已经建有水利设施。由于第NN 行的城市靠近沙漠，是该国的干旱区，所以要求其中的每座城市都建有水利设施。那么，这个要求能否满足呢？如果能，请计算最少建造几个蓄水厂；如果不能，求干旱区中不可能建有水利设施的城市数目。

输入格式

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

输出格式

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

输入输出样例

输入 #1复制

2 5
9 1 5 4 3
8 7 6 1 2

输出 #1复制

1
1

输入 #2复制

3 6
8 4 5 6 4 4
7 3 4 3 3 3
3 2 2 1 1 2

输出 #2复制

1
3

说明/提示

【样例1 说明】

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

【样例2 说明】

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

【数据范围】

---

思路

先将地图第1排每个点用dfs搜一遍，然后标记它能到达的点，判断是否能覆盖完。

如果能覆盖全部点，上方的每一个点所对应覆盖的将会是一个区间，于是我们将问题转化为区间覆盖问题，这里用dp来实现。

#include <stdio.h>
#include <iostream>
#include <memory.h>
#define N 501
#define inf 2e9+7
using namespace std;
int n,m,s,a[N][N],dp[N],cnt;
bool vis[N][N],check[N];
struct lxydl
{
	int l,r;
}cover[N];
void dfs(int x,int y,int from)
{
	if(x==n)
	{
		check[y]=1;
		cover[from].l=min(cover[from].l,y);
		cover[from].r=max(cover[from].r,y);
	}
	vis[x][y]=1;
	if(x!=n && vis[x+1][y]==0 && a[x+1][y]<a[x][y]) dfs(x+1,y,from);
	if(x!=1 && vis[x-1][y]==0 && a[x-1][y]<a[x][y]) dfs(x-1,y,from);
	if(y!=m && vis[x][y+1]==0 && a[x][y+1]<a[x][y]) dfs(x,y+1,from);
	if(y!=1 && vis[x][y-1]==0 && a[x][y-1]<a[x][y]) dfs(x,y-1,from);
}
signed main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	register int i,j;
	cin>>n>>m;
	for(i=1;i<=m;i++)
	{
		cover[i].l=dp[i]=inf;
	}
	for(i=1;i<=n;i++)
	{
		for(j=1;j<=m;j++)
		{
			cin>>a[i][j];
		}
	}
	for(i=1;i<=m;i++)
	{
		dfs(1,i,i);
		memset(vis,0,sizeof(vis));
	}
	for(i=1;i<=m;i++)
	{
		if(!check[i])
		{
			cnt++;
		}
	}
	if(cnt)
	{
		cout<<0<<endl<<cnt<<endl;
		return 0;
	}
	/*for(i=1;i<=m;i++)
	{
		cout<<cover[i].l<<' '<<cover[i].r<<endl;
	}*/
	cout<<1<<endl;
	for(i=1;i<=m;i++)//求最少点数 
	{
		for(j=1;j<=m;j++)
		{
			if(i>=cover[j].l && i<=cover[j].r)
			{
				dp[i]=min(dp[i],dp[cover[j].l-1]+1);
			}
		}
	}
	cout<<dp[m]<<endl;
	return 0;
}