hdu 3143 Speedy Escape 二分+搜索

来源：互联网发布：中国遥感卫星数据现状编辑：程序博客网时间：2024/06/06 06:58

PS. 训练赛的时候看完题就想到做法了，写完之后自觉很对，但是无限WA，实在无解，赛后重写一遍，继续无限WA，换G++交，神奇的过了，（改%lf和%f，C++都过不了）唉，无法理解啊。

英文很长，但很简单，题意不赘述了。

首先求到兄弟到所有点的最短路（不能经过警察局），然后求到警察到所有点的最短路。

然后二分速度，dfs验证可行性，若到达当前点的时间小于等于警察到这个点的时间，则可以扩展。（由于是double，所以小于或者小于等于的无所谓，不必纠结这些细节）。

精度很低，1e-5都能过=。=

G++交

#include <cstdio>#include <iostream>#include <algorithm>#include <cstring>#include<queue>#define lson num<<1,s,mid#define rson num<<1|1,mid+1,e#define maxn 105using namespace std;#define eps 1e-6#define maxn 105#define maxm 10005struct node{    int v,w,next;}edge[maxm];int head[maxn];bool in[maxn];int d[maxn];int ed[maxn];int a,b,c,n,m,e,id;int bro,pli;void add(int a,int b,int c){    edge[id].v=b;    edge[id].w=c;    edge[id].next=head[a];    head[a]=id++;}void init(){    memset(head,-1,sizeof(int)*(n+1));    memset(ed,0,sizeof(ed));    id=0;}void SPFA(int st){    queue<int> Q;    memset(in,0,sizeof(int)*(n+1));    memset(d,0x3f,sizeof(int)*(n+1));    Q.push(st);    in[st]=1;    d[st]=0;    int tmp;    while(!Q.empty())    {        tmp=Q.front();Q.pop();        in[tmp]=0;        for(int i=head[tmp];i!=-1;i=edge[i].next)        {            if(edge[i].v==pli) continue;            if(d[edge[i].v]>d[tmp]+edge[i].w)            {                d[edge[i].v]=d[tmp]+edge[i].w;                if(!in[edge[i].v])                {                    in[edge[i].v]=1;                    Q.push(edge[i].v);                }            }        }    }}int vis[maxn];int db[maxn];double mid;bool isok(int now){    vis[now]=1;    if(d[now]*mid<db[now]*160.0) return false;    if(ed[now]) return true;    for(int i=head[now];i!=-1;i=edge[i].next)    {        int to=edge[i].v;        if(vis[to]) continue;        if(isok(to)) return true;    }    return false;}int main(){    while(~scanf("%d%d%d",&n,&m,&e))    {        init();        for(int i=1;i<=m;i++)        {            scanf("%d%d%d",&a,&b,&c);            add(a,b,c);            add(b,a,c);        }        for(int i=1;i<=e;i++)        {            scanf("%d",&a);            ed[a]=1;        }        scanf("%d%d",&bro,&pli);        SPFA(bro);        for(int i=1;i<=n;i++) db[i]=d[i];        SPFA(pli);        int tzf=0;        double l=0,r=10000000;        while(r-l>=eps)        {            memset(vis,0,sizeof(vis));            mid=(l+r)/2;            if(isok(bro)) {tzf=1;r=mid;}            else l=mid;        }        if(!tzf) puts("IMPOSSIBLE");        else printf("%.10lf\n",mid);    }    return 0;}

4 0