码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj 1189 二分+最大流

时间:2018-11-26 20:58:49      阅读:192      评论:0      收藏:0      [点我收藏+]

标签:题目   else   add   type   lin   head   div   com   href   

题目传送门

思路:

  先预处理出每个人到每扇门的时间,用门作为起点进行bfs处理。

  然后二分时间,假设时间为x,将每扇门拆成1到x,x个时间点,表示这扇门有几个时间点是可以出去的。对于一扇门,每个时间点都向后一个时间点建边,表示人在当前时间点到达,可以在下一时间点出去。

  先将s连上所有的空地,流量为1,建立每个空地每个门的对应的时间点流量为1的边,表示这个空地的人会再某一时间点到达这扇门。然后每个门流向t,流量为inf。只要最大流为空地的数量,则代表该时间是可以的,继续向下二分。

#include<bits/stdc++.h>
#define CLR(a,b) memset(a,b,sizeof(a))
using namespace std;
typedef long long ll;
int fx[4][2]={{0,1},{0,-1},{1,0},{-1,0}};
int dis[100][30][30],n,m;
char s[30][30];
int ed[1000],g[30][30],cnt,sum,st,t;
const int maxn=40010;
const int inf=0x3f3f3f3f;
int head[maxn],tot=0;
struct edge{
    int to,f,Next;
}a[maxn<<1];
void addv(int u,int v,int w){
    a[++tot].to=v,a[tot].f=w,a[tot].Next=head[u],head[u]=tot;
    a[++tot].to=u,a[tot].f=0,a[tot].Next=head[v],head[v]=tot;
}
inline void bfs(int id){
    queue<int >q;
    dis[id][ed[id]/m][ed[id]%m]=0;
    q.push(ed[id]);
    while(!q.empty()){
        int u=q.front();
        int x=u/m,y=u%m;
        q.pop();
        for(int i=0;i<4;i++){
            int xx=x+fx[i][0],yy=y+fx[i][1];
            if(xx<0||xx>=n||yy<0||yy>=m||s[xx][yy]!=.)continue;
            if (dis[id][xx][yy]>dis[id][x][y]+1)
                dis[id][xx][yy]=dis[id][x][y]+1,q.push(xx*m+yy);

        }
    }
}
int deep[maxn];
bool vis[maxn];
inline void find(){
    CLR(deep,0);
    queue<int >q;
    vis[st]=1;
    q.push(st);
    while(!q.empty()){
        int u=q.front();
        q.pop();
        for(int i=head[u];i!=-1;i=a[i].Next)
        {
            if(a[i].f && !vis[a[i].to]){
                q.push(a[i].to);
                deep[a[i].to]=deep[u]+1;
                vis[a[i].to]=1;
            }
        }
    }
}
int dfs(int u,int delta){
    if(u==t)return delta;
    int ret=0;
    for(int i=head[u];delta&&i!=-1;i=a[i].Next){
        if(a[i].f&&deep[a[i].to]==deep[u]+1){
            int dd=dfs(a[i].to,min(a[i].f,delta));
            a[i].f-=dd;
            a[i^1].f+=dd;
            delta-=dd;
            ret+=dd;
        }
    }
    return ret;
}
inline bool check(int x){
    st=0,t=maxn-2,tot=1;
    CLR(head,-1);
    for(int i=1;i<=sum;i++)addv(st,i,1);
    for(int k=1;k<=cnt;k++)
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
                if(s[i][j]==.&&dis[k][i][j]<=x)
                    addv(g[i][j],sum+(k-1)*x+dis[k][i][j],1);
    
    for (int i=1; i<=cnt; i++)
        for (int j=1; j<=x; j++) {
            int tmp=(i-1)*x+sum;
            addv(tmp+j,t,1);
            if (j<x) addv(tmp+j,tmp+j+1,inf);
        }
    int ret=0;
    while(1){
        CLR(vis,0);
        find();
        if(!vis[t]){
            return ret==sum;
        }
        ret+=dfs(st,inf);
    }
}
int main(){
    cin>>n>>m;
    for(int i=0;i<n;i++)
    {
        scanf("%s",s[i]);
    }
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<m;j++)
        {
            if(s[i][j]==D){
                ed[++cnt]=i*m+j;
            }
            if(s[i][j]==.)g[i][j]=++sum;
        }
    }
    CLR(dis,0x3f);
    for(int i=1;i<=cnt;i++)
    {
        bfs(i);
    }
    int l=1,r=n*m,mid;
    int ans;
    if(!check(r)){
        puts("impossible");
        return 0;
    }
    while(l<=r){
        
        mid=(l+r)>>1;//printf("l:%d  r:%d  mid:%d\n",l,r,mid);
        if(check(mid)){
        ans=mid,r=mid-1;
    //    printf("ans:%d\n",ans);
    }
        else l=mid+1;
    }
    cout<<ans<<endl;
    
}

1189: [HNOI2007]紧急疏散evacuate

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 3825  Solved: 1118
[Submit][Status][Discuss]

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

bzoj 1189 二分+最大流

标签:题目   else   add   type   lin   head   div   com   href   

原文地址:https://www.cnblogs.com/mountaink/p/10022786.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!