首页 > 代码库 > 【BZOJ】1189: [HNOI2007]紧急疏散evacuate(二分+bfs+网络流)

【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

Sample Output

3

HINT

 

Source

 

 

【BZOJ】1189: [HNOI2007]紧急疏散evacuate(二分+bfs+网络流)