L2-001 紧急救援

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

代码长度限制 16 KB 时间限制 200 ms 内存限制 64 MB

输入格式:

输入第一行给出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 <bits/stdc++.h>
using namespace std;
const int N=510;
int dis[N],g[N][N];
int n,m,s,d;
int num[N],xjy[N],jy[N];
int path[N];
bool st[N];
void dist()
{
    memset(dis,0x3f,sizeof dis);//路径长度初始化
    memset(num,0,sizeof num);//路径个数初始化，表示到每个点的路径个数都是0
    memset(xjy,0,sizeof xjy);//新的救援队的数量
    memset(path,-1,sizeof path);//路径地址
    dis[s]=0;
    num[s]=1;
    xjy[s]=jy[s];
    for(int i=0;i<n;i++)
    {
        int t=-1;
        for(int j=0;j<n;j++)
        {
            if(!st[j]&&(t==-1||dis[j]<dis[t]))
                t=j;
        }
        st[t]=true;
        for(int j=0;j<n;j++)
        {
            if(dis[j]>dis[t]+g[t][j])
            {
                dis[j]=dis[t]+g[t][j];
                num[j]=num[t];  //使用上一条路径的个数
                xjy[j]=xjy[t]+jy[j];//对救援数量及及进行更新
                path[j]=t;//路径更新
            }
            else if(dis[j]==dis[t]+g[t][j])//这里添加判断条件，进行修改救援数量，并判断是否要修改路径
            {
                num[j]+=num[t];//要加上上一个点的全部路径个数
                if(xjy[t]+jy[j]>xjy[j])
                {
                    xjy[j]=xjy[t]+jy[j];
                    path[j]=t;//进行路径更新
                }
            }
        }
    }

    
}
int main()
{
    cin>>n>>m>>s>>d;
    memset(g,0x3f,sizeof g);
    for(int i=0;i<n;i++) cin>>jy[i];
    while(m--)
    {
        int a,b,c;
        cin>>a>>b>>c;
        //这里是无向图，需要同时存入两条边，否则会被卡
        g[a][b]=min(g[a][b],c);
        g[b][a]=min(g[b][a],c);
    }
    dist();
    cout<<num[d]<<" "<<xjy[d]<<endl;
    vector<int> a;
    for(int i=d;i!=-1;i=path[i])
        a.push_back(i);
    while(a.size())
    {
        if(a.size()>1) cout<<a.back()<<" ";
        else cout<<a.back();
        a.pop_back();
    }
    cout<<endl;


    return 0;
}