dij最短路 变形

L2-1. 紧急救援

时间限制 200 ms
内存限制 65536 kB
代码长度限制 8000 B
判题程序 Standard 作者 陈越

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

输入格式：

输入第一行给出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>
#define INF 0x3f3f3f3f
using namespace std;
int n,m,s,d;
int a[1000];
int vis[1000];
int ma[1000][1000];
int dis[1000];          ///最短路
int pre[1000];        ///记录路径
int sum[1000];       ///最短路径的价值
int cnt[1000];       ///保存出现了几种最短径
void path(int t)
{
    if(pre[t]==-1)
        return ;
    path(pre[t]);
    printf("%d ",pre[t]);
}
void dij()
{
    for(int i=0; i<n; i++)
        dis[i]=INF;              ///必为INF 如果为ma[i][s]则错,会使记录路径时s点漏掉　　　　＊＊＊＊＊　　　　注意＊＊＊＊＊

    vis[s]=1;
    dis[s]=0;
    sum[s]=a[s];
    cnt[s]=1;

    for(int i=0; i<n; i++)
    {
        int mi=INF;
        int k=s;
        for(int j=0; j<n; j++)
        {
            if(!vis[j]&&dis[j]<mi)
            {
                mi=dis[j];
                k=j;
            }
        }
        /// 首个k一定为s
        vis[k]=1;
        for(int j=0; j<n; j++)
        {
            if(!vis[j])
            {
                if(dis[j] > dis[k] + ma[k][j])
                {
                    dis[j] = dis[k] + ma[k][j];
                    sum[j] = sum[k] + a[j];
                    cnt[j] = cnt[k];  //方式等于到达前驱的方式
                    pre[j] = k;
                }
                else if(dis[j] == dis[k] + ma[k][j])
                {
                    cnt[j] = cnt[k] + cnt[j]; //两种方式距离相等，则方式数相加
                    if(sum[j] < sum[k] + a[j])    //能获得更大的救援，则更新
                    {
                        sum[j] = sum[k] + a[j];
                        pre[j] = k;
                    }
                }
            }
        }
    }
    printf("%d %d\n",cnt[d],sum[d]);
    path(d);
    printf("%d\n",d);
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&s,&d);
    memset(ma,INF,sizeof(ma));
    memset(vis,0,sizeof(vis));
    memset(pre,-1,sizeof(pre));
    for(int i=0; i<=n; i++)
        ma[i][i]=0;
    for(int i=0; i<n; i++)
        scanf("%d",&a[i]);
    for(int i=0; i<m; i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        if(ma[x][y]>z)
            ma[x][y]=ma[y][x]=z;
    }
    dij();
    return 0;
}