PTA L2-001 紧急救援

题目

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

输入格式:

输入第一行给出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计算最短路径，再用dfs搜索，计算最短路径的条数并保存路径和经过的最多救援队数（可以进行剪枝优化）

AC代码

#include<stdio.h>
#include<string.h>
#define inf 1<<29
#include<algorithm>
using namespace std;
int n,m,s,d;
int C[505][505],minl;//存放两个城市之间的关系
int V[505],v;//存放每个城市的救援队数量
int dis[505];//到每一个点的最短距离
int vis[505];//标记
int line[505],lo,tem[505];//记录路径
int num;
void dijkstra()
{
    memset(vis,0,sizeof(vis));
    vis[s]=1;
    for(int i=0; i<n; ++i)
        dis[i]=C[s][i];
    for(int i=1; i<n; ++i) //进行n-1轮（除起点外有n-1个点）
    {
        int mi=inf,mid;
        for(int j=0; j<n; ++j)
            if(!vis[j]&&dis[j]<mi)
                mi=dis[mid=j];//找到每一轮未走过的最近点（不可能有另一条路径绕到该点）
        vis[mid]=1;
        for(int i=0; i<n; ++i)
            if(dis[i]>dis[mid]+C[mid][i])//可以绕路使距离更短
                dis[i]=dis[mid]+C[mid][i];
    }
    minl=dis[d];
    return;
}
void dfs(int now,int ll,int res,int step)
{
    tem[step]=now;
    if(ll>minl)
        return;
    if(ll==minl&&now==d)//找到最短路径
    {
        ++num;
        if(v<res)
        {
            lo=step;
            v=res;
            memcpy(line,tem,sizeof(tem));
        }
        return;
    }
    for(int i=0; i<n; ++i)
    {
        if(!vis[i])
        {
            vis[i]=1;
            dfs(i,ll+C[now][i],res+V[i],step+1);
            vis[i]=0;
        }
    }
}
int main()
{
    int i,j;
    scanf("%d %d %d %d",&n,&m,&s,&d);
    for(int i=0; i<n; ++i)
        scanf("%d",&V[i]);
    for(int i=0; i<n; ++i)
        for(int j=0; j<n; ++j)
            if(i==j)
                C[i][j]=0;
            else
                C[i][j]=inf;
    for(int i=0; i<m; ++i)
    {
        int n1,n2,l;
        scanf("%d %d %d",&n1,&n2,&l);
        if(l<C[n1][n2])
        C[n1][n2]=C[n2][n1]=l;
    }
    dijkstra();
    memset(vis,0,sizeof(vis));
    v=0;
    num=0;
    vis[s]=1;
    dfs(s,0,V[s],0);
    printf("%d %d\n",num,v);
    for(int i=0; i<=lo; ++i)
    {
        if(i)
            putchar(' ');
        printf("%d",line[i]);
    }
    putchar('\n');
    return 0;
}