P3191 [HNOI2007]紧急疏散EVACUATE（费用流）

P3191 [HNOI2007]紧急疏散EVACUATE

费用流+卡常优化

我们只关心一个人通过门时的时间，在空地的行走时间可以分层维护

于是根据时间分层，到门的时候再计算代价，即代价$=$层数

每经过$1$单位时间就向下走一层

然后就是优化：

1. 删去多余点（层）：只要开和人数相同的层，因为一个人的等待时间一定小等于人数

2. 能预处理的尽量预处理

3. register，快读，inline等

4. 能不用stl尽量不用（使用stl:queue=TLE），建议手写队列

复杂度$O($卡常能过$)$，详情可看code

#include<iostream>
#include<cstring>
#include<cstdio>
#define re register
using namespace std;
#define N 200005
#define M 3000005
const int d1[]={,,,-};
const int d2[]={,,-,};
char a[][];
int n,m,nm,mx,ans,tot,S,T,v[N],p[N],id[][];
bool vis[N]; int h[M],L,R;
int Cnt=,hd[N],nxt[M],ed[N],poi[M],con[M],val[M];
void adde(int x,int y,int v1,int v2){
    nxt[ed[x]]=++Cnt; hd[x]=hd[x]?hd[x]:Cnt;
    ed[x]=Cnt; poi[Cnt]=y; con[Cnt]=v1; val[Cnt]=v2;
}
inline void link(int x,int y,int v1,int v2){adde(x,y,v1,v2),adde(y,x,,-v2);}
#define to poi[i]
bool bfs(){//普通费用流
    memset(v,,sizeof(v)); int inf=v[];
    h[L=]=S; R=; vis[S]=; v[S]=;
    while(L!=R){
        re int x=h[L++]; vis[x]=;
        if(L>=M) L=;
        for(re int i=hd[x];i;i=nxt[i])
            if(con[i]&&v[to]>v[x]+val[i]){
                v[to]=v[x]+val[i]; p[to]=i;
                if(!vis[to]){
                    vis[to]=,h[R++]=to;
                    if(R>=M) R=;
                }
            }
    }if(v[T]==inf) return ;
    tot-=; ans=max(ans,v[T]);
    for(re int u=T;u!=S;u=poi[p[u]^])
        con[p[u]]-=,con[p[u]^]+=;
    return ;
}//因为每次流量均为1，可以省去流量数组
void draw(int x,int y){
    re int p=id[x][y];
    if(a[x][y]=='X') return ;
    if(a[x][y]=='D') for(re int i=;i<=mx;++i) link(p+i*nm,T,,i);//分层与终点连边，代价=层数
    if(a[x][y]=='.'){
        link(S,p,,); ++tot;
        for(re int i=;i<mx;++i){
            link(p+i*nm,p+(i+)*nm,1e9,);
            for(re int j=;j<;++j){
                int rx=x+d1[j],ry=y+d2[j];
                if(a[rx][ry]=='X') continue;
                link(p+i*nm,id[rx][ry]+(i+)*nm,1e9,);//与四周连边
            }
        }
    }
}
int main(){
    scanf("%d%d",&n,&m);  S=N-; T=N-; nm=n*m;
    for(re int i=;i<=n;++i) scanf("%s",a[i]+);
    for(re int i=;i<=n;++i)
        for(re int j=;j<=m;++j)
            mx+=(a[i][j]=='.'),id[i][j]=(i-)*m+j;//mx：人数，建图需要的层数
    for(re int i=;i<=n;++i) for(re int j=;j<=m;++j) draw(i,j);
    while(bfs()) ;
    if(tot) puts("impossible");
    else printf("%d",ans);
    return ;
}