PAT L2-001. 紧急救援 (最短路变形 Dijkstra + 记录路径)

L2-001. 紧急救援

时间限制：200 ms

内存限制：65536 kB

代码长度限制：8000 B

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

输入格式：

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

题目链接：https://www.patest.cn/contests/gplt/L2-001

题目分析：以前做的最短路变形都是边权有两个值，这题是边权一个值，点权一个值，在最短路的情况下，求最短路上点权和最大的路径，方法类似，记录路径用path数组，path[u]表示u的前驱结点

#include <cstdio>
#include <algorithm>
using namespace std;
int const MAX = 505;
int const INF = 0x3fffffff;
int mp[MAX][MAX], val[MAX], path[MAX], dis[MAX], re[MAX], totval[MAX], pathnum[MAX];
bool vis[MAX];
int n, m, s, d;

void Dijkstra(int v0)
{
    for(int i = 0; i < n; i++)
        dis[i] = INF;
    vis[v0] = true;
    dis[v0] = 0;
    totval[v0] = val[v0];
    pathnum[v0] = 1;
    for(int i = 0; i < n; i++)
    {
        if(mp[v0][i] != INF && i != v0)
        {
            dis[i] = mp[v0][i];
            path[i] = v0;
            totval[i] = val[v0] + val[i];
            pathnum[i] = 1;
        }
    }
    for(int i = 0; i < n - 1; i++)
    {
        int mi = INF, mival = 0, u = v0;
        for(int j = 0; j < n; j++)
        {
            if(!vis[j] && dis[j] < mi)
            {
                mi = dis[j];
                u = j;
            }
        }
        vis[u] = true;
        for(int j = 0; j < n; j++)
        {
            if(!vis[j])
            {
                if(dis[u] + mp[u][j] < dis[j])
                {
                    pathnum[j] = pathnum[u];
                    dis[j] = dis[u] + mp[u][j];
                    totval[j] = totval[u] + val[j];
                    path[j] = u;
                }
                else if(dis[u] + mp[u][j] == dis[j])
                {
                    pathnum[j] += pathnum[u];
                    if(totval[j] < totval[u] + val[j])
                    {
                        totval[j] = totval[u] + val[j];
                        path[j] = u;
                    }
                }
            }
        }
    }
}

int main()
{
    scanf("%d %d %d %d", &n, &m, &s, &d);
    for(int i = 0; i < n; i++)
        scanf("%d", &val[i]);
    for(int i = 0; i < n; i++)
        for(int j = 0; j < n; j++)
            mp[i][j] = INF;
    int x, y, l;
    for(int i = 0; i < m; i++)
    {
        scanf("%d %d %d", &x, &y, &l);
        mp[x][y] = min(mp[x][y], l);
        mp[y][x] = mp[x][y];
    }
    Dijkstra(s);
    int num = 0, cur = d;
    while(cur != s)
    {
        re[num ++] = cur;
        cur = path[cur];
    }
    re[num ++] = s;
    printf("%d %d\n", pathnum[d], totval[d]);
    for(int i = num - 1; i > 0; i--)
        printf("%d ", re[i]);
    printf("%d\n", re[0]);
}