L2-001. 紧急救援（多级最短路 spfa）

L2-001. 紧急救援

时间限制

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>
using namespace std;
#define inf 0x3f3f3f3f
struct node
{
    int to;
    int w;
    int next;
} s[550*550];
int head[550*550];
int vis[550];
int dis[550];
int pre[550];
int val[550];
int sum[550];
//int cn[550];
int tmp[550];
int n, m;
int cnt;
int ss,e;
void add(int u, int v, int w)
{
    s[cnt].to = v;
    s[cnt].w = w;
    s[cnt].next = head[u];
    head[u] = cnt++;
}
void spfa(int sc)
{
    memset(vis, 0, sizeof(vis));
    memset(sum, 0, sizeof(sum));
    memset(cn, 0, sizeof(cn));
    memset(dis, inf, sizeof(dis));
    memset(pre, -1, sizeof(pre));
    dis[sc] = 0;
    vis[sc] = 1;
    sum[sc] = val[sc];
    //cn[sc] = 1;
    queue<int>q;
    q.push(sc);
    while(!q.empty())
    {
        int u = q.front();
        q.pop();
        vis[u] = 0;
        for(int i = head[u]; i != -1; i = s[i].next)
        {
            int v = s[i].to;
            int w = s[i].w;
            if(dis[v] > dis[u] + w)
            {
                dis[v] = dis[u] + w;
                sum[v] = sum[u] + val[v];
                //cn[v] = cn[u];
                pre[v] = u;
                if(!vis[v])
                {
                    vis[v] = 1;
                    q.push(v);
                }
            }
            else if(dis[v] == dis[u] + w)
            {
                dis[v] = dis[u] + w;
                //cn[v] += cn[u];
                if(sum[v] < sum[u] + val[v])
                {
                    sum[v] = sum[u] + val[v];
                    pre[v] = u;
                    if(!vis[v])
                    {
                        vis[v] = 1;
                        q.push(v);
                    }

                }

            }

        }
    }
}
int ra;
void dfs(int now)
{
    if(now == e)
    {
        ra++;return;
    }
    for(int i = head[now]; i != -1; i = s[i].next)
    {
        int v =s[i].to;
        if(!vis[v] && dis[v] == dis[now] + s[i].w)
        {
            vis[v] = 1;
            dfs(v);
            vis[v] = 0;
        }
    }
}
int main()
{
    int i,j;

    scanf("%d%d%d%d", &n, &m, &ss, &e);
    for( i = 0; i < n; i++)
        scanf("%d", &val[i]);
    cnt = 0;
    memset(head, -1, sizeof(head));
    while(m--)
    {
        int u,v,w;
        scanf("%d%d%d", &u, &v, &w);
        add(u, v, w);
        add(v, u, w);
    }
    spfa(ss);
    ra = 0;
    memset(vis, 0, sizeof(vis));
    dfs(ss);
    printf("%d %d\n",ra,sum[e]);
    int pos = e;
    int cc = 0;
    while(pos != -1)
    {
        tmp[cc++] = pos;
        pos = pre[pos];
    }
    for(int i = cc - 1; i > 0; i--)
        printf("%d ",tmp[i]);
    printf("%d\n",tmp[0]);
    return 0;
}