最大流——BZOJ1189/Luogu3191 [HNOI2007]紧急疏散evacuate

题面：BZOJ1189 Luogu3191
首先我们预处理出每个空地到每扇门的最短时间
然后我们二分时间t，然后把每扇门拆成t个点
接下来建图：超级原点->每块空地（流1），每扇门的每一个时刻->超级汇点（流1）
然后对于每块空地，枚举能够在t时间内到达的门，连上边（流1，空地->门（对应的最短时间那个时刻点））
然后等待的问题只要每扇门的某一时刻点向下一时刻点连上就行了（流INF）
然后跑最大流判断是否满流（流为所有人）就好了

#include <cstdio>#include <algorithm>#include <cmath>#include <cstring>#include <iostream>#include <ctime>#include <map>#include <queue>#include <cstdlib>#include <string>#include <climits>#include <set>#include <vector>using namespace std;struct ppap{int x,y;};const int dx[5]={0,1,-1,0,0};const int dy[5]={0,0,0,1,-1};bool vis[401][401];int d[201][201][201],n,m,s,t,mn=0,rp[401][401];char S[401];int nedge=-1,p[500001],c[500001],nex[500001],head[500001];int dist[500001],cnt=0,cur[500001];inline void addedge(int x,int y,int z){    p[++nedge]=y;c[nedge]=z;nex[nedge]=head[x];head[x]=nedge;}inline bool bfs(int s,int t){    queue<int>q;q.push(s);    memset(dist,-1,sizeof dist);dist[s]=1;    while(!q.empty()){        int now=q.front();q.pop();        for(int k=head[now];k>-1;k=nex[k])if(c[k]&&dist[p[k]]==-1){            dist[p[k]]=dist[now]+1;q.push(p[k]);        }    }    return dist[t]>-1;}inline int dfs(int x,int low){    if(x==t)return low;    int a,used=0;    for(int k=cur[x];k>-1;k=nex[k])if(c[k]&&dist[p[k]]==dist[x]+1){        a=dfs(p[k],min(c[k],low-used));        if(a)c[k]-=a,c[k^1]+=a,used+=a;        if(c[k])cur[x]=k;        if(used==low)break;    }    if(!used)dist[x]=-1;    return used;}inline int dinic(){    int flow=0;    while(bfs(s,t)){        for(int i=s;i<=t;i++)cur[i]=head[i];        flow+=dfs(s,1e9);    }    return flow;}inline void spfa(int co,int x,int y){    queue<ppap>q;q.push((ppap){x,y});    memset(vis,0,sizeof vis);vis[x][y]=1;d[co][x][y]=0;    while(!q.empty()){        ppap now=q.front();q.pop();        for(int i=1;i<5;i++){            int xx=now.x+dx[i],yy=now.y+dy[i];            if(xx<1||yy<1||xx>n||yy>m||rp[xx][yy]!=0)continue;            if(d[co][xx][yy]>d[co][now.x][now.y]+1){                d[co][xx][yy]=d[co][now.x][now.y]+1;                if(!vis[xx][yy])vis[xx][yy]=1,q.push((ppap){xx,yy});            }        }vis[now.x][now.y]=0;    }}inline bool check(int x){    s=0;t=n*m+mn*x+1;    nedge=-1;memset(nex,-1,sizeof nex);memset(head,-1,sizeof head);    for(int i=1;i<=n;i++)        for(int j=1;j<=m;j++)if(rp[i][j]==0)addedge(s,(i-1)*m+j,1),addedge((i-1)*m+j,s,0);    for(int i=1;i<=mn;i++)        for(int j=1;j<=x;j++)addedge(n*m+(i-1)*x+j,t,1),addedge(t,n*m+(i-1)*x+j,0);    for(int i=1;i<=mn;i++)        for(int j=1;j<x;j++)addedge(n*m+(i-1)*x+j,n*m+(i-1)*x+j+1,1e9),addedge(n*m+(i-1)*x+j+1,n*m+(i-1)*x+j,0);    for(int k=1;k<=mn;k++)        for(int i=1;i<=n;i++)            for(int j=1;j<=m;j++)if(rp[i][j]==0&&d[k][i][j]<=x)addedge((i-1)*m+j,n*m+(k-1)*x+d[k][i][j],1),addedge(n*m+(k-1)*x+d[k][i][j],(i-1)*m+j,0);    return dinic()==cnt;}int main(){    memset(rp,-1,sizeof rp);    scanf("%d%d",&n,&m);    for(int i=1;i<=n;i++){        scanf("%s",S+1);        for(int j=1;j<=m;j++)if(S[j]=='.')rp[i][j]=0,cnt++;        else if(S[j]=='D')rp[i][j]=++mn;    }    for(int i=1;i<=mn;i++)for(int j=1;j<=n;j++)for(int k=1;k<=m;k++)d[i][j][k]=1e9;    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)if(rp[i][j]>0)spfa(rp[i][j],i,j);    int l=0,r=1000,ans=1e9;    while(l<=r){        int mid=l+r>>1;        if(check(mid))ans=min(ans,mid),r=mid-1;        else l=mid+1;    }    if(ans!=1e9)printf("%d",ans);    else puts("impossible");    return 0;}