【BZOJ】1189: [HNOI2007]紧急疏散evacuate（二分+bfs+网络流）
http://www.lydsy.com/JudgeOnline/problem.php?id=1189
表示完全不会QAQ。。。。
于是膜拜题解orz
二分时间。。。。。。。。。。。
于是转换成判定性问题：即如何在有限时间内通过。
假设当前有t时间可供通过。。。那么每一个门最多能通过t个人。。。。。。。。
然后将所有能够到达门的点连边，容量为无限。。。。
然后源向每个可行点连边。。容量1。。
然后每个门向汇连边。。容量为t。。
然后判断即可。。。。
（一开始bfs写错了啊QAQ。。
#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long ll;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define error(x) (!(x)?puts("error"):0)
#define rdm(x, i) for(int i=ihead[x]; i; i=e[i].next)
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<'0'||c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0'&&c<='9'; c=getchar()) r=r*10+c-'0'; return k*r; }
const int N=415, oo=~0u>>1;
const int dx[4]={1, -1, 0, 0}, dy[4]={0, 0, -1, 1};
int ihead[N], cnt=1;
struct dat { int next, to, from, cap; }e[N*N*2];
void add(int u, int v, int c) {
e[++cnt].next=ihead[u]; ihead[u]=cnt; e[cnt].from=u; e[cnt].to=v; e[cnt].cap=c;
e[++cnt].next=ihead[v]; ihead[v]=cnt; e[cnt].from=v; e[cnt].to=u; e[cnt].cap=0;
}
int p[N], d[N], gap[N], cur[N];
int isap(int s, int t, int n) {
for1(i, 0, n) p[i]=0, d[i]=0, gap[i]=0, cur[i]=ihead[i];
gap[0]=n; int ret=0, f, u=s, i;
while(d[s]<n) {
for(i=cur[u]; i; i=e[i].next) if(e[i].cap && d[e[i].to]+1==d[u]) break;
if(i) {
p[e[i].to]=cur[u]=i; u=e[i].to;
if(u==t) {
for(f=oo; u!=s; u=e[p[u]].from) f=min(f, e[p[u]].cap);
for(u=t; u!=s; u=e[p[u]].from) e[p[u]].cap-=f, e[p[u]^1].cap+=f;
ret+=f;
}
}
else {
if(!(--gap[d[u]])) break;
d[u]=n; cur[u]=ihead[u];
for(i=ihead[u]; i; i=e[i].next) if(e[i].cap && d[e[i].to]+1<d[u]) d[u]=d[e[i].to]+1;
++gap[d[u]];
if(u!=s) u=e[p[u]].from;
}
}
return ret;
}
int dis[N][23][23], n, tot, m, a[30][30], S, T;
struct Q { int x, y; }q[N];
inline int ID(int x, int y) { return (x-1)*m+y; }
void bfs(int x, int y, int d[23][23]) {
int front=0, tail=0;
for1(i, 1, n) for1(j, 1, m) d[i][j]=oo;
q[tail].x=x, q[tail].y=y; ++tail;
d[x][y]=0;
Q tp;
while(front!=tail) {
tp=q[front++];
x=tp.x, y=tp.y;
rep(i, 4) {
int fx=dx[i]+x, fy=dy[i]+y;
if(fx<1 || fy<1 || fx>n || fy>m || a[fx][fy]!=1 || d[fx][fy]!=oo) continue;
d[fx][fy]=d[x][y]+1;
q[tail].x=fx; q[tail].y=fy; ++tail;
}
}
//for1(i, 1, n) { for1(j, 1, m) printf("%d ", d[i][j]==oo?0:d[i][j]); puts(""); } puts("");
}
void build() {
for1(i, 1, n) for1(j, 1, m) if(a[i][j]==2)
bfs(i, j, dis[ID(i, j)]);
}
bool check(int t) {
for1(i, 1, T) ihead[i]=0;
cnt=1;
for1(i, 1, n) for1(j, 1, m) if(a[i][j]==2) {
int id=ID(i, j);
for1(x, 1, n) for1(y, 1, m) if(a[x][y]==1 && dis[id][x][y]<=t)
add(ID(x, y), id, oo);
add(id, T, t);
}
for1(i, 1, n) for1(j, 1, m) if(a[i][j]==1) add(S, ID(i, j), 1);
return isap(S, T, T)==tot;
}
int main() {
read(n); read(m); S=n*m+1, T=S+1;
for1(i, 1, n) for1(j, 1, m) { char c=getchar(); while(c!='X'&&c!='.'&&c!='D') c=getchar(); if(c=='.') a[i][j]=1, ++tot; if(c=='D') a[i][j]=2; }
build();
int l=0, r=n*m+1, mid;
while(l<=r) {
mid=(l+r)>>1;
if(check(mid)) r=mid-1;
else l=mid+1;
}
if(!check(r+1)) puts("impossible");
else printf("%d\n", r+1);
return 0;
}
Description
发生了火警，所有人员需要紧急疏散！假设每个房间是一个N M的矩形区域。每个格子如果是'.'，那么表示这是一块空地；如果是'X'，那么表示这是一面墙，如果是'D'，那么表示这是一扇门，人们可以从这儿撤出房间。已知门一定在房间的边界上，并且边界上不会有空地。最初，每块空地上都有一个人，在疏散的时候，每一秒钟每个人都可以向上下左右四个方向移动一格，当然他也可以站着不动。疏散开始后，每块空地上就没有人数限制了（也就是说每块空地可以同时站无数个人）。但是，由于门很窄，每一秒钟只能有一个人移动到门的位置，一旦移动到门的位置，就表示他已经安全撤离了。现在的问题是：如果希望所有的人安全撤离，最短需要多少时间？或者告知根本不可能。
Input
输入文件第一行是由空格隔开的一对正整数N与M，3<=N <=20，3<=M<=20，以下N行M列描述一个N M的矩阵。其中的元素可为字符'.'、'X'和'D'，且字符间无空格。
Output
只有一个整数K，表示让所有人安全撤离的最短时间，如果不可能撤离，那么输出'impossible'（不包括引号）。
Sample Input
5 5
XXXXX
X...D
XX.XX
X..XX
XXDXX
XXXXX
X...D
XX.XX
X..XX
XXDXX
Sample Output
3
HINT
Source