PATA 1030 Travel Plan

题记

这道题我想记录的地方是这道题需要记录最短路径的路线。其他地方跟另一道题思路基本一样，点击查阅PATA 1003 Emergency
思路:
1.首先考虑如何记录;假设当前到原点s最短距离的点是u，u再往外走一步的结点记为v,u的最短距离是d[u],c[u]指到u的最优路径的花费，cost[u][v]是u到v的花费，pre[v]记录到结点v的最优路径的前驱结点;
如果d[u]+G[u][v]<d[v]，那么更新v的前驱结点pre[v]=u;
如果d[u]+G[u][v]==d[v]并且c[u]+cost[u][v]<v[v]那么也更新v的前驱结点pre[v]=u;
2.如何输出最短路径。因为我们是根据终点确定路径的，而且记录的是前驱，所以需要倒着输出，利用递归可以很容易实现，递归边界是当前访问的结点是起点，直接输出，若不是接着dfs然后输出结点v

代码如下

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;

const int INF=0x3f3f3f3f;
const int Maxv=500+10;

//输入
int N,M,S,D;
int G[Maxv][Maxv],cost[Maxv][Maxv];
//最优解的费用,最优解的前缀结点
int c[Maxv],pre[Maxv],d[Maxv];
bool visit[Maxv]={false};

void dijstra(int s){//起点为s
    fill(d,d+Maxv,INF);
    fill(c,c+Maxv,INF);
    for(int i=0;i<N;i++)
        pre[i]=i;
    d[s]=0;
    c[s]=0;
    for(int i=0;i<N;i++){
        int u=-1,MIN=INF;
        //找当前离起点距离最近的点
        for(int j=0;j<N;j++){
            if(visit[j]==false&&d[j]<MIN){
                u=j;
                MIN=d[j];
            }
        }
        if(u==-1)//剩下的结点与起点不连通
            return;
        visit[u]=true;
        for(int v=0;v<N;v++){
            if(visit[v]==false||G[u][v]!=INF){
                if(d[u]+G[u][v]<d[v]){
                    d[v]=d[u]+G[u][v];
                    c[v]=c[u]+cost[u][v];
                    pre[v]=u;
                }
                else if(d[u]+G[u][v]==d[v]){
                    if(c[u]+cost[u][v]<c[v]){
                        c[v]=c[u]+cost[u][v];
                        pre[v]=u;
                    }
                }
            }
        }
    }
}

void dfs(int v){//v是当前结点
    if(v==S){
        printf("%d ",v);
        return;
    }
    dfs(pre[v]);
    printf("%d ",v);
}

int main()
{
    scanf("%d %d %d %d",&N,&M,&S,&D);
    fill(G[0],G[0]+Maxv*Maxv,INF);
    for(int i=0;i<M;i++){
        int u,v;
        scanf("%d %d",&u,&v);
        scanf("%d",&G[u][v]);
        G[v][u]=G[u][v];
        scanf("%d",&cost[u][v]);
        cost[v][u]=cost[u][v];
    }

    dijstra(S);
    dfs(D);
    printf("%d %d\n",d[D],c[D]);
    return 0;
}