1189: [HNOI2007]紧急疏散evacuate

题目链接

题目大意：给定一个m*n的地图，每个点有可能是空地、墙或者出口，每个空地初始站着一个人，每一时刻可以向周围走1格，门每一时刻只能通过一个人，求最短多少时间后所有人可以撤离

题解：以每个门为起点bfs，求出其到所有人的最短路，bfs时不能经过门

将门按照时间拆点（限制1s通行1人），然后可以二分答案或者枚举答案

$记每个人为i，每个门为j[x]，其中x为时间$

$连(st,i,1)和(i,j[0],1)，含义显然$

设当前枚举时间为time

$连(j[time-1],j[time],INF)，(j[time],T,1)，显然……$

直至满流

我没有用num[]数组进行GAP优化，但是跑的仍然挺快……

我的收获：新技能get

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
using namespace std;

const int M=234577;
#define T(x) (x)*n*m
#define INF 0x7fffffff
#define pii pair<int,int>
const int dx[]={0,1,-1,0,0};
const int dy[]={0,0,0,1,-1};

int n,m,t,now,tim,cnt,tot,st,ed,door;
int id[25][25],dis[25][25];
int num[M],d[M],head[M],last[M];
bool Exit;
char s[25][25];

queue<pii> q;

struct edge{int to,c,nex;}e[M*10];

void add(int u,int v,int w){e[t]=(edge){v,w,head[u]};last[u]=head[u]=t++;}

void insert(int x,int y,int z){add(x,y,z),add(y,x,0);};

int dfs(int x,int in)
{
    if(x==ed) return in;
    int ans=0,f;
    for(int i=last[x];i!=-1;last[x]=i=e[i].nex)
    {
        int v=e[i].to;
        if(e[i].c&&d[v]==d[x]-1){
            f=dfs(v,min(in-ans,e[i].c));
            ans+=f;e[i].c-=f,e[i^1].c+=f;
            if(Exit||ans==in) return ans;
        }
    }
    if(--num[d[x]]==0) Exit=1;
    d[x]++,num[d[x]]++,last[x]=head[x];
    return ans;
}

int isap(int tim)
{
    memset(d,0,sizeof(d));
    Exit=0;
    int flow=0;
    while(!Exit) flow+=dfs(st,INF);
    return flow;
}

void build()
{
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            if(s[i][j]=='D')
            {
                if(tim) insert(id[i][j]+T(tim-1),id[i][j]+T(tim),INF);
                insert(id[i][j]+T(tim),ed,1);
            }
}

void getdis(int x,int y)
{
    memset(dis,0,sizeof(dis));
    q.push(make_pair(x,y));
    while(!q.empty()){
        int nx=q.front().first,ny=q.front().second;q.pop();
        for(int k=1;k<=4;k++){
            int fx=nx+dx[k],fy=ny+dy[k];
            if(dis[fx][fy]||fx<1||fx>n||fy<1||fy>m||s[fx][fy]!='.') continue;
            dis[fx][fy]=dis[nx][ny]+1;
            insert(id[fx][fy],id[x][y]+T(dis[fx][fy]-1),1);
            q.push(make_pair(fx,fy));
        }
    } 
}

void work()
{
    while(now<tot)
    {
        if(tim>=n*m){puts("impossible");return ;}
        build();
        tim++;now+=isap(tim);
    }
    cout<<tim<<endl;
}

void init()
{
    memset(head,-1,sizeof(head));
    memset(last,-1,sizeof(last));
    cin>>n>>m;
    st=0,ed=T(n*m)+1;
    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++) id[i][j]=++cnt;//！！！！ 
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            if(s[i][j]=='.') tot++,insert(st,id[i][j],1);
            if(s[i][j]=='D') getdis(i,j);
        }
}

int main()
{
    init();
    work();
    return 0;
}