Luogu3191 HNOI2007 紧急疏散 二分答案、最大流

传送门

题意：自己去看

---

考虑二分答案。$BFS$预处理出每一个人到每一扇门的最短时间，设二分的值为$mid$，那么把门拆成$mid$个点，每一个点代表第$1,2,...,mid$秒时的状态。$i-1$时刻的门向$i$时刻的门连一条流量为$INF$的边，表示有无限多的人可以在门口等待。每一个门拆出来的$mid$个点向汇点连流量为$1$的边表示逃出的一个人；源点向每一个人连一条流量为$1$的边，每一个人对应的点向每一扇门的最早到达时刻对应的点连一条流量为$1$的边，最后计算最大流是否等于人数即可。

 #include<bits/stdc++.h>
 #define MAXN 7010
 #define XX now.x + dir[i][0]
 #define YY now.y + dir[i][1]
 #define INF 0x7fffffff
 using namespace std;

 struct Edge{
     int end , upEd , flow;
 }Ed[MAXN << ];
 struct node{
     int x , y , flo;
 }now;
 const int dir[][] = {,,,-,,,-,};
 int head[MAXN] , flo[MAXN] , cur[MAXN] , cntEd , cntBlock , N , M , cntDoor , dis[][] , mid;
 map < pair < int , int > , int > door;
 char c[][];
 bool vis[][] , be[MAXN];
 queue < int > q1;
 queue < node > q2;

 void bfs(int x , int y){
     cntBlock++;
     memset(dis[cntBlock] , 0x3f , sizeof(dis[cntBlock]));
     q2.push((node){x , y , });
     memset(vis ,  , sizeof(vis));
     vis[x][y] = ;
     while(!q2.empty()){
         now = q2.front();
         q2.pop();
         for(int i =  ; i <  ; i++)
             if(c[XX][YY] != 'X' && !vis[XX][YY]){
                 vis[XX][YY] = ;
                 if(c[XX][YY] == 'D')
                     dis[cntBlock][door.find(make_pair(XX , YY))->second] = now.flo + ;
                 else
                     q2.push((node){XX , YY , now.flo + });
             }
     }
 }

 bool div(){
     while(!q1.empty())
         q1.pop();
     memset(be ,  , sizeof(be));
     q1.push();
     be[] = flo[] = ;
     while(!q1.empty()){
         int t = q1.front();
         q1.pop();
         for(int i = head[t] ; i ; i = Ed[i].upEd){
             if(!be[Ed[i].end] && Ed[i].flow){
                 flo[Ed[i].end] = flo[t] + ;
                 be[Ed[i].end] = ;
                 if(Ed[i].end == cntBlock + cntDoor * mid + ){
                     memcpy(cur , head , sizeof(head));
                     return ;
                 }
                 q1.push(Ed[i].end);
             }
         }
     }
     return ;
 }

 bool dinic(int now){
     if(now == cntBlock + cntDoor * mid + )
         return ;
     for(int &i = cur[now] ; i ; i = Ed[i].upEd)
         if(flo[Ed[i].end] == flo[now] +  && Ed[i].flow)
             if(dinic(Ed[i].end)){
                 Ed[i].flow--;
                 Ed[i ^ ].flow++;
                 return ;
             }
     return ;
 }

 inline void addEd(int a , int b , int c){
     Ed[++cntEd].end = b;
     Ed[cntEd].flow = c;
     Ed[cntEd].upEd = head[a];
     head[a] = cntEd;
 }

 bool check(){
     memset(head ,  , sizeof(head));
     cntEd = ;
     int T = cntDoor * mid + cntBlock + ;
     for(int i =  ; i < cntDoor ; i++)
         for(int j =  ; j < mid ; j++){
             addEd(mid * i + j + cntBlock , mid * i + j + cntBlock +  , INF);
             addEd(mid * i + j + cntBlock +  , mid * i + j + cntBlock , );
         }
     for(int i =  ; i < cntDoor ; i++)
         for(int j =  ; j <= mid ; j++){
             addEd(mid * i + j + cntBlock , T , );
             addEd(T , mid * i + j + cntBlock , );
         }
     int cnt = ;
     for(int i =  ; i < N ; i++)
         for(int j =  ; j < M ; j++)
             if(c[i][j] == '.'){
                 cnt++;
                 for(int k =  ; k <= cntDoor ; k++)
                     if(dis[cnt][k] <= mid){
                         addEd(cnt , mid * (k - ) + dis[cnt][k] + cntBlock , );
                         addEd(mid * (k - ) + dis[cnt][k] + cntBlock , cnt , );
                     }
                 addEd( , cnt , );
                 addEd(cnt ,  , );
             }
     int ans = ;
     while(div())
         while(dinic())
             ans++;
     return ans == cntBlock;
 }

 int main(){
     cin >> N >> M;
     for(int i =  ; i <= N ; i++)
         for(int j =  ; j <= M ; j++){
             cin >> c[i][j];
             if(c[i][j] == 'D')
                 door.insert(make_pair(make_pair(i , j) , ++cntDoor));
         }
     for(int i =  ; i < N ; i++)
         for(int j =  ; j < M ; j++)
             if(c[i][j] == '.')
                 bfs(i , j);
     int L =  , R = N * M;
     while(L < R){
         mid = L + R >> ;
         check() ? R = mid : L = mid + ;
     }
     if(R == N * M)
         cout << "impossible";
     else
         cout << R;
     return ;
 }