Subway -铁路-步行

最新推荐文章于 2024-03-15 09:30:18 发布

ID_BePosit

最新推荐文章于 2024-03-15 09:30:18 发布

阅读量326

点赞数

CC 4.0 BY-SA版权

分类专栏：最短路

本文链接：https://blog.youkuaiyun.com/BePosit/article/details/82805248

最短路专栏收录该内容

28 篇文章

订阅专栏

本文深入解析了POJ-2502地铁路径规划问题的算法实现，通过Dijkstra算法计算最短路径，特别关注地铁线路的特殊约束条件，如不能直接跳过站点。文章详细介绍了如何构建图模型，设置速度参数，并使用C++代码实现了算法，最终输出经过四舍五入处理后的最短时间。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Subway
POJ - 2502
注意建边就好了
round四舍五入
地铁呢是不能从第一站蹦到除第二站外的站的，只能一站一站走。。。。

#include<iostream>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<vector>
#include<math.h>
using namespace std;
#define maxn 100000
#define inf 0x3f3f3f3f
double v1=10000.0/60;
double v2=40000.0/60;
struct node
{
    double x;
    double y;
} vetx[maxn];
vector<node>mmp[maxn];
int sx,sy,ex,ey,cnt;
double dis[maxn],x,y;
bool vis[maxn];
double cal(double x1,double y1,double x2,double y2)
{
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
void dij()
{
    for(int i=0; i<=cnt; i++)
    {
        dis[i]=inf;
        vis[i]=0;
    }
    dis[0]=0;
    for(int i=0; i<=cnt; i++)
    {
        double temp=inf;
        int k;
        for(int j=0; j<=cnt; j++)
            if(dis[j]<temp&&!vis[j])
                temp=dis[k=j];
        vis[k]=1;
        int len=mmp[k].size();
        for(int j=0; j<len; j++)
        {
            if(dis[int (mmp[k][j].x)]>dis[k]+mmp[k][j].y&&!vis[int (mmp[k][j].x)])
                dis[int (mmp[k][j].x)]=dis[k]+mmp[k][j].y;
        }
    }
}
int main()
{
    ios::sync_with_stdio(false);
    cin>>sx>>sy>>ex>>ey;
    vetx[0].x=sx;
    vetx[0].y=sy;
    while(cin>>x>>y)
    {
        vetx[++cnt].x=x;
        vetx[cnt].y=y;
        while(cin>>x>>y)
        {
            if(x==-1&&y==-1)
                break;
            vetx[++cnt].x=x;
            vetx[cnt].y=y;
            mmp[cnt].push_back((node)
            {
                cnt-1,cal(vetx[cnt].x,vetx[cnt].y,vetx[cnt-1].x,vetx[cnt-1].y)/v2
            });
            mmp[cnt-1].push_back((node)
            {
                cnt,cal(vetx[cnt].x,vetx[cnt].y,vetx[cnt-1].x,vetx[cnt-1].y)/v2
            });
//            cout<<cnt<<" "<<cnt-1<<endl;
        }
    }
    vetx[++cnt].x=ex;
    vetx[cnt].y=ey;
    for(int i=0; i<cnt; i++)
        for(int j=i+1; j<=cnt; j++)
        {
            mmp[i].push_back((node)
            {
                j,cal(vetx[j].x,vetx[j].y,vetx[i].x,vetx[i].y)/v1
            });
            mmp[j].push_back((node)
            {
                i,cal(vetx[j].x,vetx[j].y,vetx[i].x,vetx[i].y)/v1
            });
        }
    dij();
//    for(int i=0; i<=cnt; i++)
    cout<<round(dis[cnt])<<endl;
    return 0;
}