题目链接:
题意:
给你一张地图,地图上面有障碍物、广场、和门(若干),门保证在边界上。现在每个广场有一个人,每个有都要通过门出去,且每个人可以每个单位时间可以走相邻的四个位置,且同一时间下,每个广场可以容纳多个人,而门的话可以一个单位时间只能通过一个人,问最少需要多少时间所有人都离开地图?如果不能离开的话输出impossible。
数据范围:
3≤n,m≤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;
}