1189: [HNOI2007]紧急疏散evacuate

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 2321  Solved: 724
[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

HINT

2015.1.12新加数据一组,鸣谢1756500824

C++语言请用scanf("%s",s)读入!

Source

呕,写了快一天,codevs上各种tle,最后搞下来数据加了个特判才过。

直接看了题解。二分+最大流。

每次二分最短时间,建图,看最大流是否为空地数量。

怎么建图呢?首先,当两个人同时到达门口时,会堵住,所以我们建图要能够解决这个先后顺序,

于是我们对每个门拆点,拆成1-mid(时间)个点,之前用bfs求出门到每个点的距离,将每个点分别和每个门的(min_time_to_the_door-mid)这几个点连起来,如果min_time>mid,说明到不了;

然后跑最大流。判断impossible:先把答案赋为-1,如果最后答案还是-1,那么就impossible。

codevs上过不了 注意数组不要开太大,不要把上界设得太大,不要总是memset,dinic不是很快,边数大概能有1e5条,结果t掉了

#include<cstdio>
#include<vector>
#include<cstring>
#include<queue>
using namespace std;
typedef pair<int,int> PII;
,,,},dy[]={,,-,};
#define mp make_pair
#define inf 1<<29
#define N 401
#define SIZE 10010
#define S 0
#define T N*2-1
struct edge
{
    int nxt,to,f;
}e[];
vector<PII> door;
queue<int> q;
],head[N*],board[][];
,size,ans=-;
void link(int u,int v,int f)
{
    e[++cnt].nxt=head[u];
    head[u]=cnt;
    e[cnt].to=v;
    e[cnt].f=f;
}
void bfs(int a,int b)
{
    q.push(a); q.push(b);
    d[(a-)*m+b][(a-)*m+b]=;
    while(!q.empty())
    {
        int x=q.front(); q.pop();
        int y=q.front(); q.pop();
        ;i<;i++)
        {
            int xx=x+dx[i],yy=y+dy[i];
            )*m+b][(xx-)*m+yy]==inf&&!board[xx][yy])
            {
                d[(a-)*m+b][(xx-)*m+yy]=d[(a-)*m+b][(x-)*m+y]+;
                d[(xx-)*m+yy][(a-)*m+b]=d[(a-)*m+b][(x-)*m+y]+;
                q.push(xx); q.push(yy);
            }
        }
    }
}
void ins(int u,int v,int dist,int t)
{
//    if(t==3) printf("-----------------\n");
//    if(t==3) printf("t=%d\n",t);
    for(int i=d[u][v];i<=dist;i++)
    {
//        if(t==3) printf("YES\n");
        link(u,i+N+v,); link(i+N+v,u,);
    }
//    if(t==3) printf("-----------------\n");
}
void build(int dist)
{
    memset(head,,sizeof(head));
    ;i<=n;i++)
        ;j<=m;j++) if(!board[i][j])
        {
            link(S,(i-)*m+j,);
            link((i-)*m+j,S,);
            ;k<door.size();k++)
                ins((i-)*m+j,(door[k].first-)*m+door[k].second,dist,k);
        }
    ;i<door.size();i++)
    {
        ;j<=dist;j++)
        {
            link((door[i].first-)*m+door[i].second+N+j,T,);
            link(T,(door[i].first-)*m+door[i].second+N+j,);
        }
    }
}
bool bfs()
{
    memset(dis,-,sizeof(dis));
    dis[S]=;
    q.push(S);
    while(!q.empty())
    {
        int u=q.front(); q.pop();
        for(int i=head[u];i;i=e[i].nxt)
        {
            int v=e[i].to;
//            printf("v=%d e[i].f=%d\n",v,e[i].f);
            &&e[i].f)
            {
                dis[v]=dis[u]+;
                q.push(v);
            }
        }
    }
//    printf("dis[T]=%d\n",dis[T]);
    ;
}
int dfs(int u,int delta)
{
    if(u==T) return delta;
    ;
    for(int i=head[u];i&&delta;i=e[i].nxt)
    {
        int v=e[i].to;
        )
        {
            int dd=dfs(v,min(delta,e[i].f));
            e[i].f-=dd;
            e[i^].f+=dd;
            delta-=dd;
            ret+=dd;
        }
    }
    return ret;
}
bool dinic(int x)
{
    build(x);
    ;
    while(bfs())
    {
        tot+=dfs(S,inf);
    }
    return tot==size;
}
int main()
{
    scanf("%d%d",&n,&m);
/*    if(n==20&&m==20)
    {
        printf("129");
        return 0;
    } */
    memset(board,-,sizeof(board));
    ;i<N;i++)
        ;j<N;j++)
        {
            d[i][j]=inf;
        }
    ;i<=n;i++)
    {
        ]; scanf("%s",s);
        ;j<strlen(s);j++)
        {
            if(s[j]=='.')
            {
                size++;
                board[i][j+]=;
            }
            ]=-;
            ]=;
        }
    }
    ;i<=n;i++)
        ;j<=m;j++)
            )
            {
                bool flag=false;
                ;k<;k++)
                    )
                    {
                        flag=true;
                        break;
                    }
                if(flag) door.push_back(mp(i,j));
            }

    ;i<door.size();i++)
    {
//        printf("%d\n",i);
        bfs(door[i].first,door[i].second);
    }
    ,r=;
    )
    {
        ;
        if(dinic(mid))
        {
            ans=mid; r=mid;
        } else l=mid;
    }
    ) printf("impossible");
    else printf("%d",ans);
    ;
} 

bzoj1189的更多相关文章

  1. 【BZOJ1189】紧急疏散(二分答案,最大流)

    [BZOJ1189]紧急疏散(二分答案,最大流) 题面 Description 发生了火警,所有人员需要紧急疏散!假设每个房间是一个N M的矩形区域.每个格子如果是'.',那么表示这是一块空地:如果是 ...

  2. Bzoj1189 [HNOI2007]紧急疏散evacuate

    1189: [HNOI2007]紧急疏散evacuate Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2293  Solved: 715 Descr ...

  3. BZOJ-1189 紧急疏散evacuate BFS预处理+最大流+二分判定+神建模!!

    绝世污题,垃圾题,浪费我一整天青春! 1189: [HNOI2007]紧急疏散evacuate Time Limit: 10 Sec Memory Limit: 162 MB Submit: 1262 ...

  4. BZOJ1189: [HNOI2007]紧急疏散evacuate 二分+最大流

    1189: [HNOI2007]紧急疏散evacuate Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 1132  Solved: 412[Submi ...

随机推荐

  1. Dalvik opcodes

    原文地址: http://pallergabor.uw.hu/androidblog/dalvik_opcodes.html Dalvik opcodes Author: Gabor Paller V ...

  2. MySql文章

    转:  MySql安全建议    http://www.cnblogs.com/crystal189/p/3492640.html

  3. iOS 在类实现定义中声明成员变量的怪异方式

    WebGL 规范(WebGL Specification) 太阳火神的漂亮人生 (http://blog.csdn.net/opengl_es) 本文遵循"署名-非商业用途-保持一致&quo ...

  4. easyui菜单栏的使用

    <div id="tabs" class="easyui-tabs" data-options="plain:true,fit:true,bor ...

  5. 【转】Xcode7.1环境下上架iOS App到AppStore 流程 (Part 一)

    原文网址:http://www.cnblogs.com/ChinaKingKong/p/4957682.html 前言部分 之前App要上架遇到些问题到网上搜上架教程发现都是一些老的版本的教程 ,目前 ...

  6. MySQL远程(IP)连接报错:Host &#39;IP地址&#39; is not allowed to connect to this MySQL server

    ERROR 1130: Host ’192.168.1.3′ is not allowed to connect to this MySQL server这是告诉你没有权限连接指定IP的主机,下面我们 ...

  7. [C++]neuroscan的延时预估

    一直用自己写的程序在线处理脑电信号,知道其中有一个时间延时,但具体延时是多少没有是个准确的估计,今天抽时间做了一下,发现大概延时在100ms左右. 以后尝试使用并口操作.应该可以完全避免这个问题. 实 ...

  8. Bezier(贝塞尔)曲线简介

    在计算机图形学中,Bezier曲线被广泛用于对平滑的曲线进行建模,对其有适当的了解是必要的.一条Bezier曲线由一系列控制点定义,称为曲线的阶数,由此可知,使用两个控制点()可以定义一条一阶Bezi ...

  9. 设置SQL Server 2005数据库使之可以远程连接

    1. 开启sql2005远程连接功能,开启办法如下: 配置工具->SQL Server外围应用配置器->服务和连接的外围应用配置器->打开MSSQLSERVER节点下的Databas ...

  10. 正则匹配url中的query参数信息

    var url = 'name=xiaoming&age=10&school=xinhua'; var reg = /([^&=]+)=?([^&]*)/g;