L2-001. 紧急救援 (Dijkstra算法打印路径)

作为一个城市的应急救援队伍的负责人，你有一张特殊的全国地图。在地图上显示有多个分散的城市和一些连接城市的快速道路。每个城市的救援队数量和每一条连接两个城市的快速道路长度都标在地图上。当其他城市有紧急求助电话给你的时候，你的任务是带领你的救援队尽快赶往事发地，同时，一路上召集尽可能多的救援队。

输入格式：

输入第一行给出4个正整数N、M、S、D，其中N（2<=N<=500）是城市的个数，顺便假设城市的编号为0~(N-1)；M是快速道路的条数；S是出发地的城市编号；D是目的地的城市编号。第二行给出N个正整数，其中第i个数是第i个城市的救援队的数目，数字间以空格分隔。随后的M行中，每行给出一条快速道路的信息，分别是：城市1、城市2、快速道路的长度，中间用空格分开，数字均为整数且不超过500。输入保证救援可行且最优解唯一。

输出格式：

第一行输出不同的最短路径的条数和能够召集的最多的救援队数量。第二行输出从S到D的路径中经过的城市编号。数字间以空格分隔，输出首尾不能有多余空格。

输入样例：

4 5 0 3
20 30 40 10
0 1 1
1 3 2
0 3 3
0 2 2
2 3 2

输出样例：

2 60
0 1 3

#include<stdio.h>
#include<math.h>
#include<stdlib.h>

#define MAX 502
#define INFI 9999999

int G[MAX][MAX],nv,ne,c1,c2; //图，结点数，边数，起点，终点
int know[MAX];              //访问标记
int distance[MAX];          //距离
int num[MAX];               //最短路径条数
int weight[MAX];           //每个点的权值
int w[MAX];                //最短路径上的点权值总和
int pre[MAX];              //存放每个结点前一个结点
void printPath( int v);
void CreateGraph();
void Dijkstra( );
int main()
{
    int i;
    scanf("%d%d%d%d",&nv,&ne,&c1,&c2);
    for( i=0; i<nv; i++)
        scanf("%d",&weight[i]);
    CreateGraph();
    Dijkstra();
    printf("%d %d
",num[c2],w[c2]);
    printPath(c2);
    return 0;
}

void CreateGraph()
{
    int i,j;
    int v1,v2;
    int dn;
    for( i=0; i<nv; i++)
        for( j=0; j<nv; j++)
            G[i][j]=INFI;
    for( i=0; i<ne; i++)
    {
        scanf("%d%d%d",&v1,&v2,&dn);
        G[v1][v2]=G[v2][v1]=dn;  //无向图
    }
}
void Dijkstra( )
{
    int i,j,k;
    int min;
    for( i=0; i<nv; i++)
    {
        know[i]=0;
        distance[i]=G[c1][i];  //更新起点到所有顶点的距离
    }
    distance[c1]=0;           //起点到起点的距离为0
    w[c1]=weight[c1];
    num[c1]=1;
    for( i=1; i<nv; i++)
    {
        k=-1;
        min = INFI;
        for( j=0; j<nv; j++)
        {
            if( !know[j] && distance[j]<min)
            {
                k=j;
                min = distance[j];
            }
        }
        if( k==-1) break;
        know[k]=1;  //寻找到最短的距离，标记该点
        for( j=0; j<nv; j++)
        {
            if( !know[j] && min+G[k][j]<distance[j])  //更新距离
            {
                distance[j] = min+G[k][j];
                num[j] = num[k];
                w[j] = w[k] + weight[j];
                pre[j]=k;
            }
            else if( !know[j] && min+G[k][j]==distance[j])
            {
                num[j]=num[j]+num[k];
                if( w[k]+weight[j]>w[j])
                {
                    w[j]=w[k]+weight[j];
                    pre[j]=k;
                }
            }
        }
    }
}

void printPath( int v)  //递归打印路径
{
    if( v==c1){
        printf("%d",v);
        return;
    }
    printPath( pre[v]);
    printf(" %d",v);
}