PAT-L2-001. 紧急救援（最短路 dijkstar）

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

分析：
其实就是最短路路径模板修改一下
求多条路，用一个数组记录多少条路通过这个点
如果最小距离改变那么 lu[i]=lu[v];
否则lu[i]=lu[i]+lu[v];

AC代码

#include<stdio.h>
#include<string.h>
#include<algorithm>
#define inf 10000000
using namespace std;
int dis[1000]; //记录到起点的距离 
int map[1000][1000]; //图 
int vis[1000];//点是否遍历 
int path[1000];//记录路径 
int n,m,source,sink;//n个点，m条路，起点，终点 
int sum[1000]; //每个点最大权值 
int lu[1000];//有几条路通过这个点 
int ren[1000];//点的初始权值 
void dijkstar()
{
    memset(vis,0,sizeof(vis));
    memset(path,-1,sizeof(path));
    for(int i=0;i<n;i++)
    {
        dis[i]=map[source][i];
        if(i!=source)
        ren[i]=sum[i]+sum[source];
        lu[i]=1;
    }
    vis[source]=1;
    for(int i=1;i<n;i++)
    {
        int min=inf;
        int u=source;
        for(int j=0;j<n;j++)
        {
            if(dis[j]<=min&&!vis[j])
            {
                if(dis[j]==min)
                {
                    if(ren[j]>ren[u])
                    u=j;
                }
                else
                {
                    min=dis[j];
                    u=j;
                }
            }
        }
        vis[u]=1;
        for(int j=1;j<n;j++)
        {
            if(dis[j]>=dis[u]+map[u][j]&&!vis[j])
            {
                if(dis[j]!=dis[u]+map[u][j])
                {
                    lu[j]=lu[u];
                    ren[j]=ren[u]+sum[j];
                    path[j]=u;
                }
                else
                {
                    if(ren[j]<ren[u]+sum[j])
                    ren[j]=ren[u]+sum[j],path[j]=u;
                    lu[j]+=lu[u];
                }
                dis[j]=dis[u]+map[u][j];
            }
        }
    }
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&source,&sink);
    for(int i=0;i<n;i++)
    scanf("%d",&sum[i]);
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
        {
            if(i!=j)
            map[i][j]=inf;
            else
            map[i][j]=0;    
        }
    }
    for(int i=0;i<m;i++)
    {
        int u,v,w;
        scanf("%d%d%d",&u,&v,&w);
        map[u][v]=w;
        map[v][u]=w;
    }
    if(n==1)
    printf("1 %d\n0\n",sum[0]);
    else
    {
    dijkstar();
    printf("%d %d\n",lu[sink],ren[sink]);
    int ans[1000];
    int t=0;
    for(int i=path[sink];i!=-1;i=path[i])
        ans[t++]=i;
    printf("%d",source);
    for(int i=t-1;i>=0;i--)
        printf(" %d",ans[i]);
    printf(" %d\n",sink);
    }
}