CSU 1808 - 地铁 - [最短路变形]

Posted tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了CSU 1808 - 地铁 - [最短路变形]相关的知识,希望对你有一定的参考价值。

题目链接:http://acm.csu.edu.cn/csuoj/problemset/problem?pid=1808

Time limit: 5000 ms Memory limit: 131072 kB

Bobo 居住在大城市 ICPCCamp。

ICPCCamp 有 n 个地铁站,用 1,2,…,n 编号。 m 段双向的地铁线路连接 n 个地铁站,其中第 i 段地铁属于 c i 号线,位于站 a i,b i 之间,往返均需要花费 t i 分钟(即从 a i 到 b i需要 t i 分钟,从 b i 到 a i 也需要 t i 分钟)。
众所周知,换乘线路很麻烦。如果乘坐第 i 段地铁来到地铁站 s,又乘坐第 j 段地铁离开地铁站 s,那么需要额外花费 |c i-c j | 分钟。注意,换乘只能在地铁站内进行。
Bobo 想知道从地铁站 1 到地铁站 n 所需要花费的最小时间。

Input

输入包含不超过 20 组数据。
每组数据的第一行包含两个整数 n,m (2≤n≤10 5,1≤m≤10 5).
接下来 m 行的第 i 行包含四个整数 a i,b i,c i,t i (1≤a i,b i,c i≤n,1≤t i≤10 9).
保证存在从地铁站 1 到 n 的地铁线路(不一定直达)。

Output

对于每组数据,输出一个整数表示要求的值。

Sample Input

3 3
1 2 1 1
2 3 2 1
1 3 1 1
3 3
1 2 1 1
2 3 2 1
1 3 1 10
3 2
1 2 1 1
2 3 1 1

Sample Output

1
3
2

 

题解:

如果只记录到某个节点x的最短路长度d[x],并且记录对应于d[x],是坐哪号线来到节点x的,这样显然是错误的。

原因比如这样的样例:

3 3
1 2 1 2
1 2 3 3
2 3 3 5

可以看出,d[x]要扩展到d[x][c],即这题的状态有两个量决定:到了节点x,最后乘坐的是c号线;

那么,如果我们把节点x用若干条边Edge(u1→x)…Edge(uk→x)来代替,那么我们就相当于把d[x]要扩展到d[x][c]了;

所以我们可以直接把边当成点,对边做最短路(用堆优化dijkstra)。

 

AC代码:

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
using namespace std;
typedef long long LL;
typedef pair<LL,int> Pair;

const LL INF=1e18;
const int maxn=1e5+10;
const int maxm=2e5+10; //无向边拆成两条有向边

int n,m;

struct Edge{
    int u,v,c;
    int next;
    LL t;
};
Edge E[maxm];
int head[maxn],ne;
void init()
{
    ne=0;
    memset(head,-1,sizeof(head));
}
void addedge(int u,int v,int c,LL t)
{
    E[ne].u=u, E[ne].v=v, E[ne].c=c, E[ne].t=t;
    E[ne].next=head[u];
    head[u]=ne++;
}

LL ans;
LL d[maxm];
bool vis[maxm];
void dijkstra(int st)
{
    priority_queue< Pair, vector<Pair>, greater<Pair> > Q;

    memset(vis,0,sizeof(vis));
    for(int i=0;i<ne;i++) d[i]=INF;
    ans=INF;

    for(int i=head[st];i!=-1;i=E[i].next)
    {
        d[i]=E[i].t;
        Q.push(Pair(d[i],i));
    }
    while(!Q.empty())
    {
        int x=Q.top().second; Q.pop();

        if(vis[x]) continue;
        vis[x]=1;
        if(E[x].v==n) ans=min(ans,d[x]);

        for(int y=head[E[x].v];y!=-1;y=E[y].next)
        {
            if(vis[y]) continue;
            if(d[y]>d[x]+E[y].t+abs(E[y].c-E[x].c))
            {
                d[y]=d[x]+E[y].t+abs(E[y].c-E[x].c);
                Q.push(Pair(d[y],y));
            }
        }
    }
}
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        init();
        for(int i=1;i<=m;i++)
        {
            int u,v,c; LL t;
            scanf("%d%d%d%lld",&u,&v,&c,&t);
            addedge(u,v,c,t);
            addedge(v,u,c,t);
        }

        dijkstra(1);
        printf("%lld\n",ans);
    }
}
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
using namespace std;
typedef long long LL;
typedef pair<LL,int> Pair;

const LL INF=1e18;
const int maxn=1e5+10;
const int maxm=2e5+10; //无向边拆成两条有向边

int n,m;

struct Edge{
    int u,v,c;
    LL t;
};
vector<Edge> E;
vector<int> G[maxn];
void init(int l,int r)
{
    E.clear();
    for(int i=l;i<=r;i++) G[i].clear();
}
void addedge(int u,int v,int c,LL t)
{
    E.push_back((Edge){u,v,c,t});
    G[u].push_back(E.size()-1);
}

LL ans;
LL d[maxm];
bool vis[maxm];
void dijkstra(int st)
{
    priority_queue< Pair, vector<Pair>, greater<Pair> > Q;

    memset(vis,0,sizeof(vis));
    for(int i=0;i<E.size();i++) d[i]=INF;
    ans=INF;

    for(int i=0;i<G[st].size();i++)
    {
        int x=G[st][i];
        d[x]=E[x].t;
        Q.push(Pair(d[x],x));
    }
    while(!Q.empty())
    {
        int x=Q.top().second; Q.pop();

        if(vis[x]) continue;
        vis[x]=1;
        if(E[x].v==n) ans=min(ans,d[x]);

        for(int i=0;i<G[E[x].v].size();i++)
        {
            int y=G[E[x].v][i];
            if(vis[y]) continue;
            if(d[y]>d[x]+E[y].t+abs(E[y].c-E[x].c))
            {
                d[y]=d[x]+E[y].t+abs(E[y].c-E[x].c);
                Q.push(Pair(d[y],y));
            }
        }
    }
}
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        init(1,n);
        for(int i=1;i<=m;i++)
        {
            int u,v,c; LL t;
            scanf("%d%d%d%lld",&u,&v,&c,&t);
            addedge(u,v,c,t);
            addedge(v,u,c,t);
        }

        dijkstra(1);
        printf("%lld\n",ans);
    }
}

注:两份代码的区别是分别用链式前向星和vector邻接表存图。

以上是关于CSU 1808 - 地铁 - [最短路变形]的主要内容,如果未能解决你的问题,请参考以下文章

[CSUOJ1808] 地铁(dijkstra,堆,边最短路)

最短路STLCSU 1808 地铁 (2016湖南省第十二届大学生计算机程序设计竞赛)

CSU1808地铁

CSU 1808 地铁

CSU 1808:地铁(Dijkstra)

CSU 1333 & Uva 12661 Funny Car Racing最短路变形+spfa算法,链式前向星建图