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

visit[pos]=1;
for(j=1;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=1;i<=1002;++i)
{
for(j=1;j<=1002;++j)
{
map1[i][j]=inf;
}
}
int a,b,c;
for(i=1;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=0;i<s;i++)
{
scanf("%d",&start[i]);
}
for(int i=0;i<d;i++)
{
scanf("%d",&end1[i]);
}
int cnt=inf;
for(int i=0;i<s;i++)
{
dijstra(start[i]);///对每个起点使用
for(int j=0;j<d;j++)///求到多个目的地的距离中的最小
{
cnt=cnt>dis[end1[j]]?dis[end1[j]]:cnt;
}
}

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

 

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

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

///SPFA shortest path faster algorithm
#include <cstdio>
#include <vector>
#include <algorithm>
#include <queue>
#define INF 0xfffffff
const int maxn = 1010;
using namespace std;
struct node{
int v,len;
node(int v=0,int len=0):v(v),len(len){}
};
int st[1010];
int goal[1010];
vector<node>G[maxn];
int minDist[maxn];
int inqueue[maxn];
int n,m;
int T,S,D;
void init(){
for(int i=0;i<maxn;i++){
minDist[i]=INF;
inqueue[i]=0;
G[i].clear();
}
}
int Dijkstra(){
queue<int >Q;
for(int i=0;i<S;i++){ ///将所有的起点 都推入队列
inqueue[st[i]]=true;
minDist[st[i]]=0;
Q.push(st[i]);
}
while(!Q.empty()){
int vex = Q.front();
Q.pop();
inqueue[vex]=0;
for(int i=0;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]=1;
Q.push(v);
}
}

}
}
int Min = INF;
for(int i=0;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=0;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=0;i<S;i++){
scanf("%d",&st[i]);
}
for(int i=0;i<D;i++){
scanf("%d",&goal[i]);
}
Dijkstra();
}
return 0;
}

 

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

 

法三：

【Floyd-Warshall算法】

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

 

#include <cstdio>
#include <vector>
#include <queue>
#define INF 0xfffff
const int maxn = 1010;
using namespace std;
int map[maxn][maxn];
int st[maxn];
int goal[maxn];
void init(){
for(int i=1;i<maxn;i++){
for(int j=1;j<maxn;j++){
if(i==j) map[i][j]=0;
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=-1;
for(int i=0;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=0;i<S;i++)
scanf("%d",&st[i]);
for(int i=0;i<D;i++)
scanf("%d",&goal[i]);
int ans=INF;
for(int i=0;i<S;i++){
for(int j=0;j<D;j++){
if(map[st[i]][goal[j]]<ans)
ans=map[st[i]][goal[j]];
}
}
printf("%d/n",ans);
}
return 0;
}