題目大意:給定一個矩陣,'.'表示空位,'X'表示人,'#'表示牆,'@'表示門,每個位置至多隻能站一個人,人不能穿越牆,人能從門中出去。每個人每分鐘只能上下左右移動一步,問最少需要多少時間讓所有的人出去。
思路:
(1):把每個點按照天數拆成d個點。
(2):添加源匯點S和T。
(3):源點向第0天地圖上人所在位置的點連一天容量爲1的邊。
(4):枚舉時間Ti
(5):每次枚舉只需向殘留網絡裏重新添加一些新的點和邊。然後再從源點向匯點跑最大流。
(6):直到前Ti天的最大流總和=總人數時,Ti就是所有的最少天數。
添加新的點和邊:
如果當前枚舉的天數爲d,那麼地圖中每個位置第d-1天的點向四周和本身第d天的點(新添加的點)連一條容量爲1 的邊。
如果當前點爲出口,則連接在第d天的這個點和出口一條容量爲1的邊
PS.一開始先用BFS判一下可達性。
My Code:
#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
#include <queue>
using namespace std;
#define maxn 40010 //點數
#define INF 0x7fffffff
int N,M,S,T,step,sum,ans;
int dir[5][2]={0,1,0,-1,1,0,-1,0,0,0};
struct Edge
{
int from,to,cap,flow;
};
vector<Edge> edges; //邊數 edges[e]和edges[e^1]互爲反向弧
vector<int> G[maxn]; //G[i][j]表示節點i的第j條邊在edges數組中的序號
bool vis[maxn]; //BFS使用
int d[maxn]; //從起點到i的距離
int cur[maxn]; //當前弧的下標
void init() //每組數據運行前記得初始化
{
edges.clear();
for(int i = 0; i < maxn; i++)G[i].clear();
}
void AddEdge(int from, int to, int cap) //加邊,依次是起點、終點、容量
{
edges.push_back((Edge){from,to,cap,0});
edges.push_back((Edge){to,from,0,0});
int tmp=edges.size();
G[from].push_back(tmp-2);
G[to].push_back(tmp-1);
}
bool BFS(int s, int t)
{
memset(vis,0,sizeof(vis));
queue<int> Q;
Q.push(s);
d[s]=0;
vis[s]=1;
while(!Q.empty())
{
int x=Q.front();Q.pop();
for(int i = 0; i < G[x].size(); i++)
{
Edge& e=edges[G[x][i]];
if(!vis[e.to]&&e.cap>e.flow) //只考慮殘量網絡中的弧
{
vis[e.to]=1;
d[e.to]=d[x]+1;
Q.push(e.to);
}
}
}
return vis[t];
}
int DFS(int x, int a, int t)
{
if(x==t||a==0)return a;
int flow=0,f;
for(int& i = cur[x]; i < G[x].size(); i++) //從上次考慮的弧
{
Edge& e = edges[G[x][i]];
if(d[x]+1==d[e.to]&&(f=DFS(e.to,min(a,e.cap-e.flow), t))>0)
{
e.flow+=f;
edges[G[x][i]^1].flow-=f;
flow+=f;
a-=f;
if(a==0)break;
}
}
return flow;
}
int Maxflow(int s, int t) //求最大流,傳參:源點,匯點
{
int flow=0;
while(BFS(s, t))
{
memset(cur,0,sizeof(cur));
flow+=DFS(s,INF, t);
}
return flow;
}
char data[20][20];
int to(int x, int y){return x*M+y;}
void Read()
{
sum=0;ans=0;
S=0;T=maxn-1;step=M*N;
for(int i = 1; i <= N; i++)scanf("%s",&data[i][1]);
for(int i = 1; i <= N; i++)
for(int j = 1; j <= M; j++)
if(data[i][j]=='X'){AddEdge(S,to(i,j),1);sum++;}
}
int out(int x, int y){return (x<1||x>N||y<1||y>M);}
int vis2[20][20];
struct node
{
int x,y;
};
int judge(int x, int y)
{
memset(vis2,0,sizeof(vis2));
queue<node> hehe;
node tmp;tmp.x=x;tmp.y=y;
hehe.push(tmp);
vis2[x][y]=1;
while(!hehe.empty())
{
tmp=hehe.front();
hehe.pop();
if(data[tmp.x][tmp.y]=='@')return 1;
for(int i = 0; i < 4; i++)
{
int dx=tmp.x+dir[i][0],dy=tmp.y+dir[i][1];
if(!out(dx,dy)&&!vis2[dx][dy]&&data[dx][dy]!='#')
{
vis2[dx][dy]=1;
hehe.push((node){dx,dy});
}
}
}
return 0;
}
int OK()
{
for(int i = 1; i <= N; i++)
for(int j = 1; j <= M; j++)
{
if(data[i][j]=='X')
{
if(!judge(i,j))return 0;
}
}
return 1;
}
int can(int t)
{
int Pre=(t-1)*step;
int Now=t*step;
for(int i = 1; i <= N; i++)
for(int j = 1; j <= M; j++)
{
if(data[i][j]=='#')continue;
int u=Pre+to(i,j);
for(int k = 0; k < 5; k++)
{
int dx=i+dir[k][0];
int dy=j+dir[k][1];
if(!out(dx,dy)&&data[dx][dy]!='#')
{
int v=Now+to(dx,dy);
AddEdge(u,v,1);
}
}
if(data[i][j]=='@')
{
int v=Now+to(i,j);
AddEdge(v,T,1);
}
}
int ch=Maxflow(S,T);
ans+=ch;
return ans==sum;
}
void solve()
{
if(sum==0){printf("0\n");return ;}
if(!OK()){printf("-1\n");return ;}
int Ti=0;
while(1)
{
Ti++;
if(can(Ti))break;
}
printf("%d\n",Ti);
}
int main()
{
while(~scanf("%d%d",&N,&M))
{
init();
Read();
solve();
}
return 0;
}