NOIP2010 引水入城

来源：互联网发布：手机短信软件哪个好编辑：程序博客网时间：2024/06/08 03:11

题目来源：https://www.luogu.org/problem/show?pid=1514

水厂只能在湖泊旁建，如果能够满足要求，那么每一个水厂其输水可达的干旱区的城市一定是连续的。

对每一个靠近湖泊的城市bfs，求出其可达的干旱区的城市的区间左右端点。

如果bfs过程中有没被遍历到的城市，则为无解。

如果有解，则问题转化为了区间覆盖问题，运用贪心或dp解决。

代码：

#include <iostream>#include <cstdio>#include <cstdlib>#include <cstring>#include <algorithm>#include <queue>using namespace std;int a[501][501];int du[501];int vis[501][501];int start[501],end[501];int f[521];int q[1000001];int main(){    int n,m;scanf("%d%d",&n,&m);    for(int i=1;i<=n;i++)    {        for(int j=1;j<=m;j++)        {            scanf("%d",&a[i][j]);        }    }    memset(start,63,sizeof(start));    memset(end,0,sizeof(end));    for(int j=1;j<=m;j++)    {        memset(vis,0,sizeof(vis));        int front=1,back=0;        vis[1][j]=1;        q[++back]=1;q[++back]=j;        while(front<=back)        {            int x=q[front++];            int y=q[front++];            if(x==n)            {                du[y]=1;                start[j]=min(start[j],y);                end[j]=max(end[j],y);            }            if(x-1>=1&&a[x-1][y]<a[x][y]&&vis[x-1][y]==0)            {                q[++back]=x-1;q[++back]=y;                vis[x-1][y]=1;            }            if(x+1<=n&&a[x+1][y]<a[x][y]&&vis[x+1][y]==0)            {                q[++back]=x+1;q[++back]=y;                vis[x+1][y]=1;            }            if(y-1>=1&&a[x][y-1]<a[x][y]&&vis[x][y-1]==0)            {                q[++back]=x;q[++back]=y-1;                vis[x][y-1]=1;            }            if(y+1<=m&&a[x][y+1]<a[x][y]&&vis[x][y+1]==0)            {                q[++back]=x;q[++back]=y+1;                vis[x][y+1]=1;            }        }    }    int cnt=0;    for(int i=1;i<=m;i++)if(du[i]==0)cnt++;    if(cnt)    {        cout<<"0"<<endl;        cout<<cnt;        return 0;    }    else cout<<"1"<<endl;    memset(f,63,sizeof(f));    f[0]=0;    for(int i=1;i<=m;i++)    {        for(int j=1;j<=m;j++)        {            if(i>=start[j]&&i<=end[j])            {                f[i]=min(f[i],f[start[j]-1]+1);            }        }    }    cout<<f[m];    return 0;}

0 0