洛谷-1514 引水入城

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

为了使居民们都尽可能饮用到清澈的湖水，现在要在某些城市建造水利设施。水利设施有两种，分别为蓄水厂和输水站。蓄水厂的功能是利用水泵将湖泊中的水抽取到所在城市的蓄水池中。
因此，只有与湖泊毗邻的第1 行的城市可以建造蓄水厂。而输水站的功能则是通过输水管线利用高度落差，将湖水从高处向低处输送。故一座城市能建造输水站的前提，是存在比它海拔更高且拥有公共边的相邻城市，已经建有水利设施。由于第N 行的城市靠近沙漠，是该国的干旱区，所以要求其中的每座城市都建有水利设施。那么，这个要求能否满足呢？如果能，请计算最少建造几个蓄水厂；如果不能，求干旱区中不可能建有水利设施的城市数目。
输入输出格式
输入格式：
每行两个数，之间用一个空格隔开。输入的第一行是两个正整数N,M，表示矩形的规模。接下来N 行，每行M 个正整数，依次代表每座城市的海拔高度。
输出格式：
两行。如果能满足要求，输出的第一行是整数1，第二行是一个整数，代表最少建造几个蓄水厂；如果不能满足要求，输出的第一行是整数0，第二行是一个整数，代表有几座干旱区中的城市不可能建有水利设施。

输入输出样例
输入样例#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

解释：首先处理联通问题，即第一行和最后一行的联通。对每个点设一个bitset，然后就可以记忆化搜索，最后得到第一行可以覆盖的区域，可以简单证明得到，如果可以全部覆盖最后一行，则必然覆盖区间为连续区间，所以最后就转换成了，用最少的区间覆盖全部区域的问题

#include<iostream>
#include<bitset>
#include<algorithm>
#define N 503
using namespace std;
int n=0,m=0;
bitset<500> dp[N][N];
bool vis[N][N]={0};
int G[N][N]={0};
int dir[4][2]={0,1,0,-1,1,0,-1,0};
struct node{
    int l,r;
}s[N];
bool cmp(node a,node b){
    if(a.l==b.l) return a.r>b.r;
    return a.l<b.l;
}
bool check(int x,int y){
    return x<=n&&x>=1&&y<=m&&y>=1; 
}
void dfs(int x,int y){
    if(vis[x][y]) return;
    vis[x][y]=1;
    if(x==n){
        bitset<500> temp;
        temp.set(y-1);
        dp[x][y]|=temp;
    }
    for(int i=0;i<4;i++){
        int dx=x+dir[i][0];
        int dy=y+dir[i][1];
        if(check(dx,dy)&&G[dx][dy]<G[x][y]){
            dfs(dx,dy);
            dp[x][y]|=dp[dx][dy];
        }
    }
}
int main(){
    ios::sync_with_stdio(false);
    //freopen("C:\\Users\\24867\\Desktop\\in.txt","r",stdin); 
    cin>>n>>m;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            cin>>G[i][j];
        }
    }
    for(int i=1;i<=m;i++) dfs(1,i);
    bitset<500> temp;
    bitset<500> foo;
    for(int i=1;i<=m;i++){
        if(dp[1][i].count()==0) continue;
        temp|=dp[1][i];
        int pos=0;
        while(pos<m&&!(dp[1][i][pos])) pos++;
        s[i].l=pos+1;
        while(pos<m&&(dp[1][i][pos])) pos++;
        s[i].r=pos;
    }
    sort(s+1,s+1+m,cmp);
    if(temp.count()!=m){
        cout<<0<<endl;
        cout<<m-temp.count()<<endl;
    }else{
        cout<<1<<endl;
        int sum=1;
        int index=0;
        while(s[index].l==0) index++;
        int final=s[index].r;
        int Max=index;
        for(int i=index+1;i<=m+1;i++){
            if(s[i].l<=final+1&&i<=m){
                if(s[i].r>s[Max].r) Max=i;
            }else{
				if(final<s[Max].r) sum++; 
            	final=s[Max].r;Max=0;
            	if(i<=m) i--;
            	
            }
        }
        cout<<sum<<endl;
    }
    return 0;
}