luogu3191 [HNOI2007]紧急疏散EVACUATE

最新推荐文章于 2024-09-14 09:03:33 发布

elijahqi

最新推荐文章于 2024-09-14 09:03:33 发布

阅读量235

点赞数

分类专栏：网络流

本文链接：https://blog.csdn.net/elijahqi/article/details/78927082

版权

网络流专栏收录该内容

134 篇文章 0 订阅

订阅专栏

博客详细解析了HNOI2007中紧急疏散问题的算法，即在火警情况下，如何通过二维矩阵模拟人员疏散，寻找最短撤离时间。输入输出格式、样例及解决方案被阐述，强调C++的输入方式，并指出应利用最短路径和最大流的方法来确定可能的疏散时间。

摘要由CSDN通过智能技术生成

http://www.elijahqi.win/2017/11/24/luogu3191-hnoi2007%e7%b4%a7%e6%80%a5%e7%96%8f%e6%95%a3evacuate/
题目描述
发生了火警，所有人员需要紧急疏散！假设每个房间是一个N M的矩形区域。每个格子如果是’.’，那么表示这是一块空地；如果是’X’，那么表示这是一面墙，如果是’D’，那么表示这是一扇门，人们可以从这儿撤出房间。已知门一定在房间的边界上，并且边界上不会有空地。最初，每块空地上都有一个人，在疏散的时候，每一秒钟每个人都可以向上下左右四个方向移动一格，当然他也可以站着不动。疏散开始后，每块空地上就没有人数限制了（也就是说每块空地可以同时站无数个人）。但是，由于门很窄，每一秒钟只能有一个人移动到门的位置，一旦移动到门的位置，就表示他已经安全撤离了。现在的问题是：如果希望所有的人安全撤离，最短需要多少时间？或者告知根本不可能。
输入输出格式
输入格式：
输入文件第一行是由空格隔开的一对正整数N与M，3<=N <=20，3<=M<=20，以下N行M列描述一个N M的矩阵。其中的元素可为字符’.’、’X’和’D’，且字符间无空格。
输出格式：
只有一个整数K，表示让所有人安全撤离的最短时间，如果不可能撤离，那么输出’impossible’（不包括引号）。
输入输出样例
输入样例#1：复制
5 5
XXXXX
X…D
XX.XX
X..XX
XXDXX
输出样例#1：复制
3
说明
2015.1.12新加数据一组，鸣谢1756500824
C++语言请用scanf(“%s”,s)读入！
题意描述十分清楚
我们首先把所有人到门的边互相建边然后用宽搜出最短路每次二分答案然后去验证
如何验证每次只把可以在该时间内到达的人向门连边同时每个门需要按照时间拆点因为1s只能出去一个人就每次从最短路到二分枚举的那个时间的门我去连边跑最大流看能否满流

#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define inf 0x3f3f3f3f
using namespace std;
int dx[]={0,1,0,-1},dy[]={-1,0,1,0};
struct node{
    int y,next,z;
}data[5500000],data1[160000];
int dis[440][440],flag[22*22],level[38000],h[38000],T,pe[440],d[440],dd,p,num,h1[22*22];
char s[22][22];int n,m;
inline void insert1(int x,int y,int z){
    data1[++num].y=y;data1[num].z=z;data1[num].next=h1[x];h1[x]=num;
}
inline void insert(int x,int y,int z){
    data[++num].y=y;data[num].z=z;data[num].next=h[x];h[x]=num;
    data[++num].y=x;data[num].z=0;data[num].next=h[y];h[y]=num;
}
inline void bfs1(int u){
    memset(flag,0,sizeof(flag));queue<int>q;flag[u]=1;q.push(u);dis[u][u]=0;
    while (!q.empty()){
        int x=q.front();q.pop();
        for (int i=h1[x];i;i=data1[i].next){
            int y=data1[i].y;if (!flag[y]) q.push(y),flag[y]=1,dis[u][y]=dis[u][x]+1;
        }
    }
}
inline bool bfs(){
    memset(level,0,sizeof(level));queue<int> q;q.push(0);level[0]=1;
    while (!q.empty()){
        int x=q.front();q.pop();
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z;
            if (level[y]||!z) continue;
            level[y]=level[x]+1;q.push(y);if (y==T) return 1;
        }
    }return 0;
}
int dfs(int x,int s){
    if (x==T) return s;int ss=s;
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y,z=data[i].z;
        if (level[x]+1==level[y]&&z){
            int xx=dfs(y,min(z,s));if (!xx) level[y]=0;
            s-=xx;data[i].z-=xx;data[i^1].z+=xx;
            if (!s) return ss;
        }
    }return ss-s;
}
inline int calc(int x,int y){return p+(x-1)*n*m+y;}
inline bool check(int mid){
    memset(h,0,sizeof(h));num=1;
    for (int i=1;i<=p;++i){
        bool flag1=0;
        for (int j=1;j<=dd;++j){
            if (dis[pe[i]][d[j]]<=mid){
                flag1=1;if (dis[pe[i]][d[j]]==inf) continue; 
                for (int k=dis[pe[i]][d[j]];k<=mid;++k) insert(i,calc(j,k),1);
            }
        }
        if (!flag1) return false;
    }int ans=0;T=p+dd*n*m+1;
    for (int i=1;i<=p;++i) insert(0,i,1);
    for (int i=1;i<=dd;++i)
        for (int j=1;j<=n*m;++j) insert(calc(i,j),T,1);
//  for (int i=2;i<=num;++i) printf("%d %d %d\n",data[i].x,data[i].y,data[i].z);
    while (bfs()) ans+=dfs(0,inf); 
    if (ans==p) return 1;else return 0;
}
int main(){
    freopen("3191.in","r",stdin);
    scanf("%d%d",&n,&m);memset(dis,0x3f,sizeof(dis));
    for (int i=1;i<=n;++i) scanf("%s",s[i]+1);
    for (int i=1;i<=n;++i){
        for (int j=1;j<=m;++j){
            if (s[i][j]=='X'||s[i][j]=='D') continue;
            for (int k=0;k<4;++k){
                int x=i+dx[k],y=j+dy[k];
                if (x<1||x>n||y<1||y>m) continue;
                if (s[x][y]=='X') continue;
                insert1((i-1)*m+j,(x-1)*m+y,1);
            }
        }
    }
    for (int i=1;i<=n;++i){
        for (int j=1;j<=m;++j){
            if (s[i][j]=='.') pe[++p]=(i-1)*m+j,bfs1(pe[p]);
            if (s[i][j]=='D') d[++dd]=(i-1)*m+j;
        }
    }

    int l=1,r=n*m;
    while (l<=r){
        int mid=l+r>>1;
        if (check(mid)) r=mid-1;else l=mid+1;
    }
    if (l>n*m) printf("impossible\n");else printf("%d\n",l);
    /*for (int i=1;i<=p;++i){
        for (int j=1;j<=dd;++j) printf("%d %d %d\n",pe[i],d[j],dis[pe[i]][d[j]]);
        printf("\n");
    }*/
    return 0;
}