多条最短路径、点权和最大问题

7-4 城市间紧急救援 (25分)

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

输入格式:

输入第一行给出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

思路：首先用二维数组来储存各个点之间的权值，用一个数组来表示每个顶点的最短路径有多少条，用一个数组来表示前一个节点。当出现两条路路径相等时，根据点权和找到最优路径，比如本题的例子，两条路径相等，判断两条路径中哪条路径的救援队多就选择哪条路径。
代码出处：http://www.manongjc.com/article/31534.html
ac代码：

#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
int n,m,s,d;
int Map[1001][1001];
int pace[1001],vis[1001],num[1001],sum[1001],pre[1001];
//pace[]数组记录每个点最短路径的条数
//vis[]标记函数
//num[]数组表每个点救援队的数量
//sum[]数组表示到每个点最多召集的救援队（包括自身点的救援队）
//pre[]数组储存前一个点
int path[501];
void init() //初始化
{
    memset(vis,0,sizeof vis);
    for(int i=0; i<n; i++)
        for(int j=0; j<n; j++)
            if(i!=j) Map[i][j]=Map[j][i]=INF;
}
void Dijkstra()
{
    pace[s]=vis[s]=1;//最短路径也要相应初始化
    for(int i=0; i<n; i++)
    {
        path[i]=Map[s][i];
    }
    for(int i=0; i<n-1; i++) //仅遍历n-1次即可；
    {
        int minn=INF,f=-1;//f赋初值-1
        for(int j=0; j<n; j++) //寻找当前未被标记的且里出发地最近的点
        {
            if(!vis[j]&&path[j]<minn)
            {
                minn=path[j];
                f=j;
            }
        }
        if(f==-1) break;
        vis[f]=1;//该点标记
        for(int j=0; j<n; j++) //以该点位中转站，对未标记的的点松弛
        {
            if(!vis[j]&&path[j]>path[f]+Map[f][j])
            {
                path[j]=path[f]+Map[f][j];//路径长度更新
                pre[j]=f;//记录上一节点位置
                pace[j]=pace[f];//当前最短路径的个数
                sum[j]=sum[f]+num[j];//召集救援队的总人数
            }
            else if(!vis[j]&&path[j]==path[f]+Map[f][j])
            {
                pace[j]+=pace[f];//当前也为最短路径，数目加一
                if(sum[j]<sum[f]+num[j]) //实时更新可召集救援队总数目
                {
                    pre[j]=f;
                    sum[j]=sum[f]+num[j];
                }
            }
        }
    }
}
int main()
{
    int u,v,l;
    cin>>n>>m>>s>>d;
    init();
    for(int i=0; i<n; i++)
    {
        scanf("%d",&num[i]);
        sum[i]=num[i];//分别为该点最多召集救援队数目和该点含救援队数目
        pace[i]=1;//最短路径数目初始为1
    }
    for(int i=0; i<m; i++)
    {
        scanf("%d%d%d",&u,&v,&l);
        if(l<Map[u][v])
        {
            Map[u][v]=Map[v][u]=l;
        }
    }
    Dijkstra();
    cout<<pace[d]<<" "<<sum[d]+num[s]<<endl;//加上出发点所含队伍
    int road[1001];
    int cur=d,cnt=0;
    road[cnt++]=d;
    while(pre[cur]!=0) //由终点向起点寻找路径
    {
        road[cnt++]=pre[cur];
        cur=pre[cur];
    }
    road[cnt++]=s;
    for(int i=cnt-1; i>=0; i--)
    {
        if(i!=cnt-1)
            cout<<" ";
        cout<<road[i];
    }
    return 0;
}