天梯赛 L2 城市间紧急救援

题目:

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

输入输出格式:

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

本题考查Dijkstra算法，且需要维护额外信息，详见代码

#include<bits/stdc++.h>

using namespace std;

const int N=510, M=N*N+100;
int n,m,s,d;
int dist[N];   //dist用于存储起始点到目标点的当前最大距离
int cnt[N];    //cnt表示有几条路经过该点（多条最短路径）
int sum[N];    //sum表示经过该点的所有路径中，其中救援队最多的一条的数量和
int num[N];    //num表示该城市初始时有多少救援队
int f[N];      //表示当前最优情况下该点前一点是谁
//两种更新最优情况：
//当dist[j]>dis+w[i]，最短路更新为当前最优
//当dist[j]==dis+w[i]，sum更新为当前最优
bool st[N];    //最短路是否已经确定
int h[N], e[M], w[M], ne[M], idx;     //邻接表存储
typedef pair<int,int> PII;

void add (int a, int b, int c)
{
    e[idx]=b, w[idx]=c, ne[idx]=h[a], h[a]=idx++;
}

void ff (int u)    //递归输出最短路径
{
    if(u==s)
    {
        cout<<u;
    }
    else
    {
        ff(f[u]);
        cout<<' '<<u;
    }
}

void dijkstra ()
{
    memset(dist,0x3f,sizeof dist);
    dist[s]=0;
    cnt[s]=1;
    sum[s]=num[s];
    //初始化起始点，最短路为0，经过路径为1，最大救援队数量为初始数量

    priority_queue<PII,vector<PII>,greater<PII>> heap;
    heap.push({0,s});
    //优先队列模拟堆

    while(heap.size())
    {
        auto t=heap.top();
        heap.pop();
        int ver=t.second, dis=t.first;

        if(st[ver]==true) continue;
        else st[ver]=true;
        //取出当前最短路长度最小的点ver，该点最短路大小即为dist[ver]
        //若ver最短路已经确定，continue

        for(int i=h[ver];i!=-1;i=ne[i])
        {
            int j=e[i];

            //更新情况1：最短路更新
            if(dist[j]>dis+w[i])
            {
                dist[j]=dis+w[i];
                cnt[j]=cnt[ver];
                sum[j]=num[j]+sum[ver];
                f[j]=ver;
                heap.push({dist[j],j});
            }
            //更新情况2：救援队最大值更新
            else if(dist[j]==dis+w[i])
            {
                if(sum[ver]+num[j]>sum[j])
                {
                    sum[j]=sum[ver]+num[j];
                    f[j]=ver;
                }
                cnt[j]+=cnt[ver];
            }
        }
    }

    cout<<cnt[d]<<' '<<sum[d]<<endl;
}

int main ()
{
    memset(h,-1,sizeof h);
    cin>>n>>m>>s>>d;
    for(int i=0;i<n;i++)
    {
        int x;
        cin>>x;
        num[i]=x;
        //存入救援队初始值
    }
    for(int i=0;i<m;i++)
    {
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
        add(b,a,c);
        //存入双向图
    }
    dijkstra();

    ff(d);

    return 0;
}