poj3057 Evacuation（二分+最大匹配）

题目链接：

poj3057

题意：

给你一张地图，地图上面有障碍物、广场、和门（若干），门保证在边界上。现在每个广场有一个人，每个有都要通过门出去，且每个人可以每个单位时间可以走相邻的四个位置，且同一时间下，每个广场可以容纳多个人，而门的话可以一个单位时间只能通过一个人，问最少需要多少时间所有人都离开地图?如果不能离开的话输出impossible。

数据范围：

$3\leq n,m \leq 12$

题解：

用网络流调了一天不知道哪里错了。。。最后改成匈牙利算法才过。。到现在也不知道哪里错了。
1、我们可以利用bfs求出每个人每个门的最短距离。
2、我们知道每个时间点只能通过一个人，我们通过枚举从小到大时间ti，然后把时间拆点，如果i点能在ti之前到达j号门，我们就把i与j号门对应的能到的时间点连边，然后利用匈牙利算法判断匹配能否到达所有人都出去。
3、其实这里我们直接二分时间就可以了，因为这是个单调的。。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<sstream>
#include<algorithm>
#include<vector>
#include<bitset>
#include<set>
#include<queue>
#include<stack>
#include<map>
#include<cstdlib>
#include<cmath>
#define PI 2*asin(1.0)
#define LL long long
#define pb push_back
#define pa pair<int,int>
#define clr(a,b) memset(a,b,sizeof(a))
#define lson lr<<1,l,mid
#define rson lr<<1|1,mid+1,r
#define bug(x) printf("%d++++++++++++++++++++%d\n",x,x)
#define key_value ch[ch[root][1]][0]
const int  MOD = 1000000007;
const int N = 25;
const int maxn = 200+ 14;
const int mm=10010;
const int letter = 130;
const int INF = 1e9;
const double pi=acos(-1.0);
const double eps=1e-8;
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int dir[4][2]={{1,0},{0,1},{-1,0},{0,-1}};
int n,m,ps[N][N],vis[N][N];
char mp[N][N];
int num,cnt,s,t;
int d[maxn][maxn];
int cx[mm],cy[mm],mark[mm],tot,nx,ny,head[mm];
struct node{
    int to,next;
}e[1000100];
void add(int u,int v){

    e[tot].to=v,e[tot].next=head[u],head[u]=tot++;
    e[tot].to=u,e[tot].next=head[v],head[v]=tot++;
    //printf("tot = %d\n",tot);
}
bool ok(int x,int y){
    return x>=0&&x<n&&y>=0&&y<m;
}
int path(int x){
    for(int i=head[x];i!=-1;i=e[i].next){
        int to=e[i].to;
        if(!mark[to]){
            mark[to]=1;
            if(cy[to]==-1||path(cy[to])){
                cx[x]=to,cy[to]=x;
                return 1;
            }
        }
    }
    return 0;
}
bool check(int ti){
    tot=0,clr(head,-1);
    int res=0;
    clr(cx,-1),clr(cy,-1);
    nx=num;
    for(int j=1;j<=num;j++)
        for(int i=1;i<=cnt;i++){
            for(int p=d[i][j];p<=ti;p++){
                add(j,num+(i-1)*ti+p);
            }
        }

    for(int i=1;i<=nx;i++){
        if(cx[i]==-1){
            clr(mark,0);
            res+=path(i);
        }
    }
    return res==num;
}
void solve(){
    int l=0,r=120;
    int ans=-1;
    while(l<=r){
        int mid=(l+r)>>1;
        if(check(mid)) ans=mid,r=mid-1;
        else l=mid+1;
    }
    if(ans==-1) puts("impossible");
    else printf("%d\n",ans);
}
void bfs(int p,int xx,int yy){
    clr(vis,-1);
    queue<pa>pq;
    while(!pq.empty()) pq.pop();
    vis[xx][yy]=0;
    pq.push(make_pair(xx,yy));
    while(!pq.empty()){
        pa now=pq.front();
        pq.pop();
        for(int k=0;k<4;k++){
            int x=now.first+dir[k][0],y=now.second+dir[k][1];
            if(ok(x,y)&&mp[x][y]=='.'&&vis[x][y]==-1){
                vis[x][y]=vis[now.first][now.second]+1;
                d[p][ps[x][y]]=vis[x][y];
                pq.push(make_pair(x,y));
            }
        }
    }
}
void init(){
    num=cnt=0;
    clr(ps,0);
    for(int i=0;i<maxn;i++)
        for(int j=0;j<maxn;j++) d[i][j]=INF;
    for(int i=0;i<n;i++)
        for(int j=0;j<m;j++) if(mp[i][j]=='.') ps[i][j]=++num;
    for(int i=0;i<n;i++)
        for(int j=0;j<m;j++) if(mp[i][j]=='D') bfs(++cnt,i,j);
}
int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        for(int i=0;i<n;i++) scanf("%s",mp[i]);
        init();
        solve();
    }
    return 0;
}