bzoj1189（的网络流判定+二分）

12 12
DDDDDDDDDDDD
D..........D
D..........D
D..........D
D..........D
D..........D
D..........D
D..........D
D..........D
D..........D
D..........D

DDDDDDDDDDDD

这个数据怎么也过不了。。。挖个坑，先做别的题，不能把时间都用在debug上

果然有的时候硬耗也并不好，我放了一放，终于找到了bug，然而又发现了自己细节上太马虎了，总是考虑不全啊

思路：求最短时间所有的人可以逃脱，那么提前预处理每一个人到各个们的时间，bfs即可

二分时间，如果一个点到一个们的level小于等于时间的话，加上这条边，

原点连向各个点，cap为1，因为要避免出现同一时间人赶到的情况，因为每一秒只能有一个人通过

那么，按时间将们的点拆点，一个点只能向大于等于他level的时间的点连，每一个分的点向汇点连边cap为1，

表示这一时间只能有一个人过去。

网络流判定，如果流量==总人数，就为可行

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<queue>
using namespace std;
const int inf=0x3f3f3f3f;
const int N=100100;
const int qx[]={0,0,1,-1};
const int qy[]={1,-1,0,0};

int n,m,s,t,map[25][25],id[25][25],drx[1000],dry[1000],num,id_dr[1000][1000];
int lev[1000][1000];
bool ok(int i,int x,int y)
{
	if (x<1||x>n||y<1||y>m) return false;
	if (map[x][y]!=1) return false;
	if (lev[i][id[x][y]]!=inf) return false;
	return true;
}
void bfs1(int i)
{
	queue<int> x;queue<int> y;
	x.push(drx[i]);y.push(dry[i]);
	lev[i][id[drx[i]][dry[i]]]=0;
	while (!x.empty())
	{
		int xx=x.front(),yy=y.front();
		x.pop();y.pop();
		for (int l=0;l<4;l++)
		if (ok(i,xx+qx[l],yy+qy[l]))
		{
			lev[i][id[xx+qx[l]][yy+qy[l]]]=lev[i][id[xx][yy]]+1;
			x.push(xx+qx[l]);y.push(yy+qy[l]);
		}
	}
}
int head[N],tot,cur[N];
struct aa
{
	int cap,flow,pre,to;
}edge[N*50];
void addedge(int x,int y,int z)
{
	edge[++tot].to=y;edge[tot].pre=head[x];edge[tot].cap=z;head[x]=tot;
	edge[++tot].to=x;edge[tot].pre=head[y];edge[tot].cap=0;head[y]=tot;
}
int le[N];
bool bfs()
{
	memset(le,0,sizeof(le));
	queue<int> q;q.push(s);
	le[s]=1;
	while (!q.empty())
	{
		int u=q.front();q.pop();
		for (int i=head[u];i;i=edge[i].pre)
		if (edge[i].cap>edge[i].flow&&!le[edge[i].to])
		{
			le[edge[i].to]=le[u]+1;
			q.push(edge[i].to);
			if (edge[i].to==t) return true;
		}
	}
	return false;
}
int dfs(int u,int maxflow)
{
	if (u==t||maxflow==0) return maxflow;
	int ans=0;
	for (int &i=cur[u];i;i=edge[i].pre)
	if (le[edge[i].to]==le[u]+1)
	{
		int flow=dfs(edge[i].to,min(maxflow,edge[i].cap-edge[i].flow));
		ans+=flow;
		maxflow-=flow;
		edge[i].flow+=flow;
		edge[((i-1)^1)+1].flow-=flow;
		if (maxflow==0) return ans;
	}
	return ans;
}
int work(int limit)
{
	int ans=0;
	memset(head,0,sizeof(head));
	memset(edge,0,sizeof(edge));
	tot=0;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++) 
		if (map[i][j]==1)
		{
			addedge(s,id[i][j],1);
			for (int l=1;l<=num;l++)
			if (lev[l][id[i][j]]<=limit) //接下，因为有可能没有被扫到的点还是初始值，若为0，则小于limit，图就错了！！！注意！！
			for (int k=lev[l][id[i][j]];k<=limit;k++)
			addedge(id[i][j],id_dr[l][k],1);
		}
	for (int i=1;i<=num;i++) 
		for (int j=1;j<=limit;j++) addedge(id_dr[i][j],t,1);
	while (bfs()) 
	{
		for (int i=s;i<=t;i++) cur[i]=head[i];
		ans+=dfs(s,inf);
	}
	return ans;
}
int main()
{
	scanf("%d%d",&n,&m);
	char ch[30];
	int bj=0,total=0;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++) id[i][j]=++bj;
	s=0;
	for (int i=1;i<=n;i++)
	{
		scanf("%s",ch);
		for (int j=0;j<strlen(ch);j++)
		switch(ch[j])
		{
			case 'X':map[i][j+1]=0;break;
			case '.':map[i][j+1]=1;total++;break;
			case 'D':map[i][j+1]=2;drx[++num]=i;dry[num]=j+1;break;
		}
	}
	memset(lev,inf,sizeof(lev));//初始时，要把所有点到们的距离置为inf，再更新！！wa了好久。。
	for (int i=1;i<=num;i++) bfs1(i);
	
	for (int i=1;i<=num;i++)
		for (int j=1;j<=total;j++) id_dr[i][j]=++bj;
	
	t=bj+1;
	int mid,l=1,r=total,ans=-1;
	while (l<=r)
	{
		mid=(l+r)>>1;
		if (work(mid)==total) {ans=mid;r=mid-1;}
		else l=mid+1;
	}
	if (ans==-1) printf("impossible");
	else printf("%d",ans);
	return 0;
}