计算几何+SPFA

我已经不想看我的提交记录了。。。。

HDU 我起码WA了2页。。。。

都是浮点数惹的祸。

const double eps=1e-4;

abs(a-b)<=eps;

总共 n 条鳄鱼，最多有 n*(n+1)/2 条路。

抽象化处理。

把 中心的起点当作 起点0 ； 最多有 n+1 条路。

把鳄鱼和周围的边界的终点都当作 n+1 ; 最多有 n+1 条

总共就只存在 n+2个点。

就是计算0 和 n+1 的最短距离。

有个小优化，就是当 跳跃距离能够直接跳到岸上的时候就直接输出 42.5 1；

G++ ,C++ 都过了。

#include<cstdio>
#include<cstring>
#include<string>
#include<queue>
#include<cmath>
#include<algorithm>
#include<queue>
#include<map>
#include<iostream>
using namespace std;
const double INF= 0xfffffff;
bool vis[103];
double dis[103];
int path[103];
int n;
double m;
const double eps=1e-4;
struct lx
{
    int v;
    double len;
};
vector<lx>g[103];
struct node
{
    double x,y;
} point[101];
void build()
{
    int u,v;
    for(int i=1; i<=n; i++)
    {
        double x=point[i].x;
        double y=point[i].y;
        double len;
        u=i;
        double xx,yy;
        lx now;
        xx=min(50.00-x,50.00-y);
        yy=min(50.00+x,50.00+y);
        len=min(xx,yy);
        if(len<=m)
        {
            now.v=n+1,now.len=len;
            g[u].push_back(now);
        }
        len=sqrt(x*x+y*y);
        if(len-7.5<=m)
        {
            now.len=len-7.5;
            now.v=u;
            g[0].push_back(now);
        }
        for(int j=i+1; j<=n; j++)
        {
            len=sqrt((x-point[j].x)*(x-point[j].x)+(y-point[j].y)*(y-point[j].y));
            if(len<=m)
            {
                now.len=len;
                now.v=j;
                g[u].push_back(now);
                now.v=u;
                g[j].push_back(now);
            }
        }
    }
}
void SPFA()
{
    for(int i=0; i<103; i++)
        dis[i]=INF,vis[i]=0,path[i]=INF;
    queue<int>q;
    dis[0]=0,vis[0]=1;
    q.push(0);path[0]=0;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int j=0; j<g[u].size(); j++)
        {
            int v=g[u][j].v;
            double len=g[u][j].len;
            if(abs(dis[v]-dis[u]-len)<=eps)
            {
                if(path[v]>path[u]+1)
                path[v]=path[u]+1;
            }
            else if(dis[v]>=dis[u]+len)
            {
                dis[v]=dis[u]+len;
                path[v]=path[u]+1;
                if(!vis[v])
                {
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
    if(abs(dis[n+1]-INF)<=eps)
        puts("can't be saved");
    else
        printf("%.2f %d\n",dis[n+1],path[n+1]);
}
int main()
{
    while(cin>>n>>m)
    {
        for(int i=0; i<103; i++)
            g[i].clear();
        for(int i=1; i<=n; i++)
            scanf("%lf%lf",&point[i].x,&point[i].y);
        if(m>=42.50)
        {
            puts("42.50 1");
            continue;
        }
        build();
        SPFA();
    }
}