bzoj1189&洛谷P3191: [HNOI2007]紧急疏散evacuate

1189: [HNOI2007]紧急疏散evacuate

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 3635   Solved: 1063
[ Submit][ Status][ Discuss]

Description

发生了火警，所有人员需要紧急疏散！假设每个房间是一个N M的矩形区域。每个格子如果是'.'，那么表示这是一

块空地；如果是'X'，那么表示这是一面墙，如果是'D'，那么表示这是一扇门，人们可以从这儿撤出房间。已知门

一定在房间的边界上，并且边界上不会有空地。最初，每块空地上都有一个人，在疏散的时候，每一秒钟每个人都

可以向上下左右四个方向移动一格，当然他也可以站着不动。疏散开始后，每块空地上就没有人数限制了（也就是

说每块空地可以同时站无数个人）。但是，由于门很窄，每一秒钟只能有一个人移动到门的位置，一旦移动到门的

位置，就表示他已经安全撤离了。现在的问题是：如果希望所有的人安全撤离，最短需要多少时间？或者告知根本

不可能。

Input

第一行是由空格隔开的一对正整数N与M，3<=N <=20，3<=M<=20，

以下N行M列描述一个N M的矩阵。其中的元素可为字符'.'、'X'和'D'，且字符间无空格。

Output

只有一个整数K，表示让所有人安全撤离的最短时间，

如果不可能撤离，那么输出'impossible'（不包括引号）。

Sample Input

5 5
XXXXX
X...D
XX.XX
X..XX
XXDXX

Sample Output

3

HINT

2015.1.12新加数据一组，鸣谢1756500824

题解：我都笑死，百度第一个题解建边是错的！！！

这题当然先二分答案mid，把每个门拆成mid个点。每个点都与汇点之间连一条容量为1的边（保证每秒只能有一个人到终点），然后对于每个门的每个i到i+1连一条（mid-i+1）的点，因为当前时间后只剩下（mid-i+1）个人可以走了，如果i到终点容量已满（即有人这一秒出去了），那他就得等下一秒或更之后了，所以说往后连边，加权值保证了不能有更多人走过。

然后每个空地到每个门之间的最短路预处理出来，算第i个空地到第j个门之间的距离为x，那么我们就对于i，j的第x个点（即第x秒达到）连一条容量为1的边，最后把每个空地与源点相连，跑最大流即可。

代码（原谅我copy）：

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
using namespace std;
const int inf=100000000,N=161000;
struct ee{int to,next,f;}e[5000001];
struct bb{int x,y;}b[1010];
int head[N],q[N*2],dis[N],dis1[410][21][21],map[21][21];
int S,T,n,m,cnt=1,ans1,w,ans,door=1,num;
char s[30];
int fx[4]={0,0,1,-1},fy[4]={1,-1,0,0}; 
void ins(int u,int v,int f){
    e[++cnt].to=v;e[cnt].f=f;e[cnt].next=head[u];head[u]=cnt;
    e[++cnt].to=u;e[cnt].f=0;e[cnt].next=head[v];head[v]=cnt;
}
  
bool bfs(){
    for (int i=1;i<=T;i++) dis[i]=inf;
    int h=0,t=1,now;
    q[1]=S;dis[S]=0;
    while(h!=t){
        now=q[++h];
        for (int i=head[now];i;i=e[i].next){
            int v=e[i].to;
            if (e[i].f&&dis[now]+1<dis[v]){
                dis[v]=dis[now]+1;
                if (v==T)return 1;
                q[++t]=v;
            }
        }
    }
    if (dis[T]==inf) return 0; return 1;
}
  
int dinic(int now,int f){
    if (now==T) return f;
    int rest=f;
    for (int i=head[now];i;i=e[i].next){
        int v=e[i].to;
        if (e[i].f&&dis[v]==dis[now]+1&&rest){
            int t=dinic(v,min(rest,e[i].f));
            if (!t) dis[v]=0;
            e[i].f-=t;
            e[i^1].f+=t;
            rest-=t;
            //if(t) printf("%d %d %d\n",now,v,e[i].f);
        }
    }
    return f-rest;
}
void bfs1(int k,int x1,int y1){
    int h,t;h=t=0;
    b[++t].x=x1,b[t].y=y1;dis1[k][x1][y1]=0;
    while(h<t){
        int x=b[++h].x,y=b[h].y;
        for(int i=0;i<4;i++){
            int nx=x+fx[i],ny=y+fy[i];
            if(map[nx][ny]==1&&dis1[k][nx][ny]>dis1[k][x][y]+1){
                dis1[k][nx][ny]=dis1[k][x][y]+1;
                b[++t].x=x+fx[i],b[t].y=y+fy[i];
            }
        }
    }
}
  
void build(int x){
    memset(head,0,sizeof(head));
    memset(e,0,sizeof(e));
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            if (map[i][j]==1) {
                int u=(i-1)*m+j;
                ins(S,u,1);
                for(int k=2;k<=door;k++)
                    if(dis1[k][i][j]<=x) ins(u,dis1[k][i][j]+(k-2)*x+m*n,x);
            }
            if(map[i][j]>1){
                for(int k=1;k<=x;k++){
                    int u=((map[i][j]-2)*x+k)+m*n;
                    if(k!=x)ins(u,u+1,x-k+1);
                    ins(u,T,1);
                }
            }
        }
} 
  
bool judge(int x){
    build(x);
    ans=0;
    while(bfs()) 
    ans+=dinic(S,inf);
    if(ans>=num) return 1;
    return 0;
         
}  
int main(){
    scanf("%d%d",&n,&m);
    S=0,T=160010;
    for(int i=1;i<=n;i++){
            scanf("%s",s+1);
            for(int j=1;j<=m;j++) {
                if (s[j]=='.')map[i][j]=1,num++;
                else if(s[j]=='D') map[i][j]=++door;
            }
        }
    for(int i=2;i<=door;i++)
        for(int j=1;j<=n;j++)
            for(int k=1;k<=m;k++) dis1[i][j][k]=inf; 
        for(int j=1;j<=n;j++)
            for(int k=1;k<=m;k++) 
                if(map[j][k]>1) bfs1(map[j][k],j,k);
    int l=0,r=400;
    ans1=-1;
    while(l<=r){
        int mid=(l+r)>>1;
        if(judge(mid)) ans1=mid,r=mid-1;
        else l=mid+1;
    }
    if(ans1==-1) printf("impossible");
    else printf("%d",ans1);
}