1689: [HNOI2007]紧急疏散evacuate

时间限制: 1 Sec  内存限制: 128 MB

题目描述

输入

输出

样例输入

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

样例输出

3

   因为每一秒门只能出一个人，所以人可能堵住。因此要把门按时间拆点，

   而答案为二分时间，只要求出空地到每扇门的最短路（最短时间），然后把此空地与门连接（时间为最短路径到二分的时间），边权为一。

    注意拆出的新点标号别重复。。身败名裂。。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<queue>
#define inf 100000000
using namespace std;
int n,m,a[404][404],cnt=0,door[404],peo[404];
int tim[404][704],adj[5000000],s=0,dep[5000000],e;
int S=0,T;
char tu[25][25];
struct node
{int v,l,next;
} lu[5000100];
void add(int u,int v,int l)
{lu[++e].v=v;lu[e].next=adj[u];adj[u]=e;lu[e].l=l;}
int bfs()
{memset(dep,0,sizeof(dep));queue<int> q;dep[S]=1;q.push(S);while(!q.empty()){int x=q.front();q.pop();for(int i=adj[x];i!=-1;i=lu[i].next){int to=lu[i].v;if(!dep[to]&&lu[i].l){dep[to]=dep[x]+1;if(to==T)return 1;q.push(to);}}}return 0;
}
int dfs(int x,int fw){if(x==T) return fw;int tmp=fw,k;for(int i=adj[x];i!=-1;i=lu[i].next){int v=lu[i].v;if(lu[i].l && tmp && dep[v]==dep[x]+1){k=dfs(v,min(tmp,lu[i].l));if(!k){dep[v]=0;continue;}lu[i].l-=k; lu[i^1].l+=k; tmp-=k;}}return fw-tmp;
}
int check(int len)
{memset(adj,-1,sizeof(adj));e=0;int sum=cnt;for(int i=1;i<=cnt;i++)if(door[i])for(int j=1;j<=len;j++)tim[i][j]=(i-1)*len+cnt;T=300000;for(int i=1;i<=cnt;i++)if(peo[i]){add(S,i,1),add(i,S,0);for(int j=1;j<=cnt;j++)if(door[j]){int k=a[i][j];for(k=a[i][j];k<=len;k++)add(i,tim[j][k],1),add(tim[j][k],i,0);}}for(int i=1;i<=cnt;i++)if(door[i]){for(int j=1;j<=len;j++)add(tim[i][j],T,1),add(T,tim[i][j],0);}int ans=0,tt;while(bfs())ans+=dfs(S,inf);if(ans>=s)return 1;return 0;
}
int main()
{
//  freopen("data.in","r",stdin);
//  freopen("data.out","w",stdout);scanf("%d%d",&n,&m);for(int i=1;i<=n;i++)scanf("%s",tu[i]+1);memset(a,30,sizeof(a));for(int i=1;i<=n;i++)for(int j=1;j<=m;j++){cnt++;if(tu[i][j]=='.'){s++;peo[cnt]=1;if(i!=1&&(tu[i-1][j]=='.'||tu[i-1][j]=='D'))a[cnt][cnt-m]=1,a[cnt-m][cnt]=1;if(i!=n&&(tu[i+1][j]=='.'||tu[i+1][j]=='D'))a[cnt][cnt+m]=1,a[cnt+m][cnt]=1;if(j!=1&&(tu[i][j-1]=='.'||tu[i][j-1]=='D'))a[cnt][cnt-1]=1,a[cnt-1][cnt]=1;if(j!=m&&(tu[i][j+1]=='.'||tu[i][j+1]=='D'))a[cnt][cnt+1]=1,a[cnt+1][cnt]=1;}if(tu[i][j]=='D'){door[cnt]=1;}}for(int k=1;k<=cnt;k++)for(int i=1;i<=cnt;i++)for(int j=1;j<=cnt;j++) if(a[i][j]>a[i][k]+a[k][j])a[i][j]=a[i][k]+a[k][j];for(int i=1;i<=cnt;i++)if(peo[i]){int p=1;for(int j=1;j<=cnt;j++)if(door[j]&&a[i][j]<1000){p=0;break;}if(p==1){printf("impossible");exit(0);}       }int l=0,r=600,mid,ans=600;while(l<=r){mid=(l+r)/2;if(check(mid))r=mid-1,ans=mid;elsel=mid+1;}printf("%d",ans);
}