PTA 拯救007(Floyd、SPFA)

最新推荐文章于 2024-07-04 09:00:00 发布

m0_51864047

最新推荐文章于 2024-07-04 09:00:00 发布

阅读量330

点赞数

分类专栏：图论天梯赛文章标签：图论算法

本文链接：https://blog.csdn.net/m0_51864047/article/details/115291932

版权

天梯赛同时被 2 个专栏收录

14 篇文章 0 订阅

订阅专栏

图论

9 篇文章 0 订阅

订阅专栏

该博客介绍了如何通过编程解决一个类似'最小瓶颈路径'的问题，即在电影《007之生死关头》中的场景，007如何从鳄鱼池中逃脱。作者提供了两种算法实现，一是Floyd算法，二是SPFA算法，来计算007能否跳过鳄鱼到达岸边。通过输入鳄鱼的数量和007的跳跃最大距离，程序会判断007是否有逃生可能。博客内容涉及到距离计算、路径瓶颈和图论算法的应用。

摘要由CSDN通过智能技术生成

PTA 拯救007 (25 分)

在老电影“007之生死关头”（Live and Let Die）中有一个情节，007被毒贩抓到一个鳄鱼池中心的小岛上，他用了一种极为大胆的方法逃脱 —— 直接踩着池子里一系列鳄鱼的大脑袋跳上岸去！（据说当年替身演员被最后一条鳄鱼咬住了脚，幸好穿的是特别加厚的靴子才逃过一劫。）

设鳄鱼池是长宽为100米的方形，中心坐标为 (0, 0)，且东北角坐标为 (50, 50)。池心岛是以 (0, 0) 为圆心、直径15米的圆。给定池中分布的鳄鱼的坐标、以及007一次能跳跃的最大距离，你需要告诉他是否有可能逃出生天。

输入格式：

首先第一行给出两个正整数：鳄鱼数量 N（≤100）和007一次能跳跃的最大距离 D。随后 N 行，每行给出一条鳄鱼的 (x,y) 坐标。注意：不会有两条鳄鱼待在同一个点上。

输出格式：

如果007有可能逃脱，就在一行中输出"Yes"，否则输出"No"。

输入样例 1：

输出样例 1：

Yes

输入样例 2：

输出样例 2：

No

这是一道类似于“最小瓶颈路径”的题。总体思路是：先算出小岛到每只鳄鱼的“最小瓶颈路径”。然后判断每一只鳄鱼，如果可以到达（瓶颈距离小于D），并且可以从这里离开“水池”，那么就说明是可以逃离的。
小岛中心在（0,0）点，小岛到每只鳄鱼的距离要减去7.5（小岛半径）。之后可以用Floyd或者SPFA来计算起点到任意点的瓶颈距离，还应特判一下是否可以从小岛上一步逃离（D>=42.5）。~~虽然不特判也能过。~~

方法一（Floyd）：

#include<cstdio>
#include<cmath>
#include<algorithm>
#define dis(a,b,c,d)    sqrt((a-c)*(a-c)+(b-d)*(b-d))
using namespace std;
const int maxn=105;
int n;
double p[maxn][2],d,mp[maxn][maxn];

int main(){
    scanf("%d%lf",&n,&d);
    for(int i=1;i<=n;i++){
        scanf("%lf%lf",&p[i][0],&p[i][1]);
        mp[0][i]=mp[i][0]=max(0.0,dis(0,0,p[i][0],p[i][1])-7.5);     //单独计算小岛到鳄鱼的距离（减去7.5）
    }
    for(int i=1;i<=n;i++)   //鳄鱼之间的距离
        for(int j=i+1;j<=n;j++)
            mp[i][j]=mp[j][i]=dis(p[i][0],p[i][1],p[j][0],p[j][1]);
    for(int k=0;k<=n;k++)
        for(int i=0;i<=n;i++)
            for(int j=0;j<=n;j++)
                mp[i][j]=min(mp[i][j],max(mp[i][k],mp[k][j]));
    for(int i=1;i<=n;i++)
        if(abs(p[i][0])+d>=50||abs(p[i][1]+d>=50))
            if(mp[0][i]<=d){ puts("Yes");return 0;}
    puts(d>=42.5?"Yes\n":"No\n");
}

方法二（SPFA）:

#include<cstdio>
#include<cmath>
#include<queue>
#define dis(a,b,c,d) sqrt((a-c)*(a-c)+(b-d)*(b-d))
using namespace std;
const int maxn=105;
int n,m,v[maxn];
double p[maxn][2],d[maxn];
queue<int> que;

void spfa(){
    for(int i=0;i<maxn;i++) d[i]=1e9;   d[0]=0;
    for(int i=1;i<=n;i++)   que.push(i),v[i]=1,d[i]=dis(0,0,p[i][0],p[i][1])-7.5;
    while(que.size()){
        int t=que.front();que.pop();
        for(int i=1;i<=n;i++){
            double z=dis(p[t][0],p[t][1],p[i][0],p[i][1]);
            if(d[i]>max(d[t],z)){
                d[i]=max(d[t],z);
                if(!v[i])   que.push(i),v[i]=1;
            }
        }
        v[t]=0;
    }
}

int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)   scanf("%lf%lf",&p[i][0],&p[i][1]);
    spfa();
    for(int i=1;i<=n;i++)
        if((abs(p[i][0])+m>=50||abs(p[i][1])+m>=50)&&d[i]<=m){
            puts("Yes");    return 0;
        }
    puts(m>=42.5?"Yes":"No");
}