HDU2066一个人的旅行---（多起点多终点最短路径）

http://acm.hdu.edu.cn/showproblem.php?pid=2066

一个人的旅行

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 36107    Accepted Submission(s): 12313

Problem Description

虽然草儿是个路痴（就是在杭电待了一年多，居然还会在校园里迷路的人，汗~),但是草儿仍然很喜欢旅行，因为在旅途中 会遇见很多人（白马王子，^0^），很多事，还能丰富自己的阅历，还可以看美丽的风景……草儿想去很多地方，她想要去东京铁塔看夜景，去威尼斯看电影，去阳明山上看海芋，去纽约纯粹看雪景，去巴黎喝咖啡写信，去北京探望孟姜女……眼看寒假就快到了，这么一大段时间，可不能浪费啊，一定要给自己好好的放个假，可是也不能荒废了训练啊，所以草儿决定在要在最短的时间去一个自己想去的地方！因为草儿的家在一个小镇上，没有火车经过，所以她只能去邻近的城市坐火车（好可怜啊~）。

 

Input

输入数据有多组，每组的第一行是三个整数T，S和D，表示有T条路，和草儿家相邻的城市的有S个，草儿想去的地方有D个；
接着有T行，每行有三个整数a，b，time,表示a,b城市之间的车程是time小时；(1=<(a,b)<=1000;a,b 之间可能有多条路)
接着的第T+1行有S个数，表示和草儿家相连的城市；
接着的第T+2行有D个数，表示草儿想去地方。

 

Output

输出草儿能去某个喜欢的城市的最短时间。

 

Sample Input

6 2 3

1 3 5

1 4 7

2 8 12

3 8 4

4 9 12

9 10 2

1 2

8 9 10

 

Sample Output

9

 

题意：

给你一些通路，告诉你通过这些通路所需的时间。

给你一些起点，一些终点。

求从任意一起点，到任意一终点，所需的最短时间

分析：多起点多终点最短路径问题

 

 法一：

多次使用经典dijkstra单源最短路径算法（求一起点，到其他所有点的最短路径,找所有终点中最小花费。依次操作完每个起点。。。）

#include<stdio.h>
#include<string.h>
#define inf 0x3f3f3f3f
int map1[][],dis[],visit[];///map1[][]存储所有边，
///dis[]存储第一个点到第i个点的路径长度，visit[]存储各点的访问状态
int start[],end1[];
const int n=;
void dijstra(int s1)
{
    int i,j,pos=,min1;
    memset(visit,,sizeof(visit));
    for(i=;i<=n;++i)
    {
        dis[i]=map1[s1][i];///第一个点到其他所有点的距离
    }
    visit[s1]=;
    dis[s1]=;
    for(i=;i<n;i++)///剩余n-1个点，广度优先，依次求取距离原点的 最短路径长度
    {
        min1=inf;
        for(j=;j<=n;++j)
        {
            if(!visit[j]&&min1>dis[j])
            {
                min1=dis[j];
                pos=j;
            }
        }///找到最小dis值

        visit[pos]=;
        for(j=;j<=n;++j)
        {
            if(!visit[j]&&dis[j]>dis[pos]+map1[pos][j])///1->j的距离是否可通过i->pos->j距离更短？
            dis[j]=dis[pos]+map1[pos][j];
        }
    }
}
int main()
{
    int i,j;
    int t,s,d;
    while(~scanf("%d%d%d",&t,&s,&d)&&t)
    {
        for(i=;i<=;++i)
        {
            for(j=;j<=;++j)
            {
                map1[i][j]=inf;
            }
        }
        int a,b,c;
        for(i=;i<=t;++i)///输入两点+权值
        {
            scanf("%d%d%d",&a,&b,&c);
            if(c<map1[a][b])
            map1[a][b]=map1[b][a]=c;///无向边
        }
        for(int i=;i<s;i++)
        {
            scanf("%d",&start[i]);
        }
        for(int i=;i<d;i++)
        {
            scanf("%d",&end1[i]);
        }
        int cnt=inf;
        for(int i=;i<s;i++)
        {
            dijstra(start[i]);///对每个起点使用
            for(int j=;j<d;j++)///求到多个目的地的距离中的最小
            {
                cnt=cnt>dis[end1[j]]?dis[end1[j]]:cnt;
            }
        }

        printf("%d\n",cnt);
    }
    return ;
}

网搜其他解法，抽时间理解学习，加注释。

法二：SPFA， 将起点都推入队列。求出各个终点的最短时间。然后再找到最短的；

///SPFA shortest path faster algorithm
#include <cstdio>
#include <vector>
#include <algorithm>
#include <queue>
#define INF 0xfffffff
const int maxn = ;
using namespace std;
struct node{
    int v,len;
    node(int v=,int len=):v(v),len(len){}
};
int st[];
int goal[];
vector<node>G[maxn];
int minDist[maxn];
int inqueue[maxn];
int n,m;
int T,S,D;
void init(){
    for(int i=;i<maxn;i++){
        minDist[i]=INF;
        inqueue[i]=;
        G[i].clear();
    }
}
int Dijkstra(){
    queue<int >Q;
    for(int i=;i<S;i++){ ///将所有的起点 都推入队列
        inqueue[st[i]]=true;
        minDist[st[i]]=;
        Q.push(st[i]);
    }
    while(!Q.empty()){
        int vex = Q.front();
        Q.pop();
        inqueue[vex]=;
        for(int i=;i<G[vex].size();i++){
            int v = G[vex][i].v;
            if(G[vex][i].len+minDist[vex]<minDist[v]){ ///松弛操作，找到最短边
                minDist[v]=G[vex][i].len+minDist[vex];
                if(!inqueue[v])
                {
                    inqueue[v]=;
                    Q.push(v);
                }
            }

        }
    }
    int Min = INF;
    for(int i=;i<D;i++){ ///在终点找到用时最小的终点
        Min = min(minDist[goal[i]],Min);
    }
    printf("%d\n",Min);
}
int main(){

    while(scanf("%d%d%d",&T,&S,&D)!=EOF){
        init();
        int a,b,len;
        for(int i=;i<T;i++){
            scanf("%d%d%d",&a,&b,&len);
            G[a].push_back(node(b,len));
            G[b].push_back(node(a,len));
        }
        for(int i=;i<S;i++){
            scanf("%d",&st[i]);
        }
        for(int i=;i<D;i++){
            scanf("%d",&goal[i]);
        }
        Dijkstra();
    }
    return ;
}

#include <cstdio>
#include <cstring>
#include <queue>
#include <vector>
#include <algorithm>
using namespace std;
int t, s, d, k;
int dis[];
bool vis[];
int goal[];
struct edge
{
    int v, w, next;
}e[];
int head[];
queue <int> q;
void adde(int u, int v, int w)
{
    e[k].v = v;
    e[k].w = w;
    e[k].next = head[u];
    head[u] = k++;
}
void input()
{
    int u, v, w;
    for (int i = ; i <= t; i++)
    {
        scanf("%d%d%d", &u, &v, &w);
        adde(u, v, w);
        adde(v, u, w);
    }
    for (int i = ; i < s; i++)
    {
        scanf("%d", &v);
        adde(, v, );
        adde(v, , );
    }
    for (int i = ; i < d; i++)
        scanf("%d", &goal[i]);
}
void init()
{
    k = ;
    memset(head, -, sizeof(head));
    memset(dis, 0x3f, sizeof(dis));
    memset(vis, , sizeof(vis));
    dis[] = ;
    while (!q.empty())
        q.pop();
}
void spfa()
{
    q.push();
    vis[] = ;
    while (!q.empty())
    {
        int u = q.front();
        q.pop();
        vis[u] = ;
        for (int i = head[u]; i != -; i = e[i].next)
        {
            int v = e[i].v;
            if (dis[v] > dis[u] + e[i].w)
            {
                dis[v] = dis[u] + e[i].w;
                if (!vis[v])
                {
                    vis[v] = ;
                    q.push(v);
                }
            }
        }
    }
}
void solve()
{
    int ans = 0x3f3f3f3f;
    spfa();
    for (int i = ; i < d; i++)
        ans = min(ans, dis[goal[i]]);
    printf("%d\n", ans);
}
int main()
{
    while (scanf("%d%d%d", &t, &s, &d) == )
    {
        init();
        input();
        solve();
    }
    return ;
}

法三：

【Floyd-Warshall算法】

效率是O（n^3） 用到了DP的思想

#include <cstdio>
#include <vector>
#include <queue>
#define INF 0xfffff
const int maxn = ;
using namespace std;
int map[maxn][maxn];
int st[maxn];
int goal[maxn];
void init(){
    for(int i=;i<maxn;i++){
        for(int j=;j<maxn;j++){
            if(i==j) map[i][j]=;
            else
                map[i][j]=INF;
        }
    }
}
int main(){
    int T,S,D;
    while(scanf("%d%d%d",&T,&S,&D)!=EOF){
        init();
        int a,b,len;
        int Min=INF,Max=-;
        for(int i=;i<T;i++){
            scanf("%d%d%d",&a,&b,&len);
            if(a>Max) Max = a;
            if(b>Max) Max = b;
            if(b<Min) Min = b;
            if(a<Min) Min =a;
            if(len<map[a][b])
            {
                map[a][b]=map[b][a]=len;
            }
        }
        for(int k=Min;k<=Max;k++){
            for(int i=Min;i<=Max;i++){
                if(map[i][k]!=INF) ///没有这句优化就超时了
                for(int j=Min;j<=Max;j++){
                    if(map[i][j]>map[i][k]+map[k][j])
                        map[i][j]=map[i][k]+map[k][j];
                }
            }
        }
        for(int i=;i<S;i++)
            scanf("%d",&st[i]);
        for(int i=;i<D;i++)
            scanf("%d",&goal[i]);
        int ans=INF;
        for(int i=;i<S;i++){
            for(int j=;j<D;j++){
                if(map[st[i]][goal[j]]<ans)
                    ans=map[st[i]][goal[j]];
            }
        }
        printf("%d/n",ans);
    }
    return ;
}