一个人的旅行(floyd+dijskra+SPFA+Bellman)

转载

mob604756ee0e6c 2015-07-23 17:30:00

文章标签 i++ sed #include #define java 文章分类 代码人生

一个人的旅行

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

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

题解：floy剪枝优化；

代码：

 1 #include<stdio.h>
 2 #include<string.h>
 3 #define MIN(x,y)(x<y?x:y)
 4 const int INF=0x3f3f3f3f;
 5 const int MAXN=1010;
 6 int T,S,D,used[MAXN];
 7 int map[MAXN][MAXN],s[MAXN],e[MAXN];
 8 int main(){
 9     int a,b,c;
10      while(~scanf("%d%d%d",&T,&S,&D)){
11          memset(used,0,sizeof(used));
12          for(int i=0;i<=MAXN;i++)
13          for(int j=0;j<=MAXN;j++)
14              if(i-j)map[i][j]=INF;
15              else map[i][j]=0;
16          while(T--){
17              scanf("%d%d%d",&a,&b,&c);
18              if(c<map[a][b])map[a][b]=map[b][a]=c;
19              used[a]++;used[b]++;
20          }
21          for(int k=0;k<MAXN;k++){
22              if(!used[k])continue;
23              for(int i=0;i<MAXN;i++){
24                  if(map[i][k]==INF)continue;
25                  for(int j=0;j<MAXN;j++){
26                      map[i][j]=MIN(map[i][j],map[i][k]+map[k][j]);
27                     }
28                 }
29             }
30         for(int i=0;i<S;i++){
31              scanf("%d",&s[i]);
32          }
33          for(int i=0;i<D;i++){
34              scanf("%d",&e[i]);
35          }
36          int ans=INF;
37          for(int i=0;i<S;i++)
38              for(int j=0;j<D;j++)
39                      ans=MIN(ans,map[s[i]][e[j]]);
40                      printf("%d\n",ans);
41      }
42     return 0;
43 }

dijskra代码：

 1 #include<stdio.h>
 2 #include<string.h>
 3 #define MIN(x,y)(x<y?x:y)
 4 const int INF=0x3f3f3f3f;
 5 const int MAXN=1010;
 6 int map[MAXN][MAXN];
 7 int vis[MAXN],d[MAXN];
 8 int N,begin[MAXN],end[MAXN];
 9 void initial(){
10     //memset(map,INF,sizeof(map));
11     memset(vis,0,sizeof(vis));
12     memset(d,INF,sizeof(d));
13 }
14 void dijskra(int s){
15     initial();
16     d[s]=0;
17     int k;
18     while(true){
19             k=-1;
20         for(int i=1;i<=1000;i++){
21             if(!vis[i]&&(k==-1||d[i]<d[k]))k=i;
22         }
23         if(k==-1)break;
24         vis[k]=1;
25         for(int i=1;i<=1000;i++){
26             d[i]=MIN(d[i],d[k]+map[k][i]);
27         }
28     }
29 }
30 int main(){
31     int T,S,D;
32     int a,b,c,temp;
33     while(~scanf("%d%d%d",&T,&S,&D)){
34             memset(map,INF,sizeof(map));
35         while(T--){
36             scanf("%d%d%d",&a,&b,&c);
37         if(c<map[a][b])map[a][b]=map[b][a]=c;
38         }
39         int ans=INF;
40         for(int i=0;i<S;i++)scanf("%d",&begin[i]);
41         for(int i=0;i<D;i++)scanf("%d",&end[i]);
42         for(int i=0;i<S;i++){
43                 dijskra(begin[i]);
44             for(int j=0;j<D;j++){
45                // printf("%d\n",temp);
46                 ans=MIN(ans,d[end[j]]);
47             }
48         }
49         printf("%d\n",ans);
50     }
51 return 0;
52 }

SPFA

 1 #include<stdio.h>
 2 #include<string.h>
 3 #include<queue>
 4 #define MIN(x,y) (x<y?x:y)
 5 using namespace std;
 6 const int INF=0x3f3f3f3f;
 7 const int MAXN=1010;
 8 const int MAXM=2010;
 9 int head[MAXM];
10 struct Edge{
11     int from,to,value,next;
12 };
13 Edge edg[MAXM];
14 int T;
15 int S;
16 int D;
17 int enu;
18 int s[MAXN],e[MAXN];
19 int vis[MAXN],dis[MAXN],used[MAXN];
20 queue<int>dl;
21 void initial(){
22     memset(dis,INF,sizeof(dis));
23     memset(vis,0,sizeof(vis));
24     while(!dl.empty()){
25         dl.pop();
26     }
27 }
28 void SPFA(int sx){
29     initial();
30     vis[sx]=1;dis[sx]=0;dl.push(sx);
31     while(!dl.empty()){
32         int k=dl.front();
33         dl.pop();
34         vis[k]=0;
35         for(int i=head[k];i!=-1;i=edg[i].next){
36             int v=edg[i].to;
37             if(dis[k]+edg[i].value<dis[v]){
38                 dis[v]=dis[k]+edg[i].value;
39                 if(!vis[v])
40                     dl.push(v),vis[v]=1;
41             }
42         }
43     }
44 }
45 void add(int u,int v,int value){
46     Edge E={u,v,value,head[u]};
47     edg[enu]=E;
48     head[u]=enu++;
49 }
50 void get(){
51     int a,b,c;
52     while(T--){
53         scanf("%d%d%d",&a,&b,&c);
54         add(a,b,c);
55         add(b,a,c);
56     }
57     for(int i=0;i<S;i++)scanf("%d",&s[i]);
58     for(int i=0;i<D;i++)scanf("%d",&e[i]);
59 }
60 void print(){
61     int ans=INF;
62     for(int i=0;i<S;i++){
63         SPFA(s[i]);
64         for(int j=0;j<D;j++){
65             ans=MIN(ans,dis[e[j]]);
66         }
67     }
68     printf("%d\n",ans);
69 }
70 int main(){
71     while(~scanf("%d%d%d",&T,&S,&D)){
72         memset(head,-1,sizeof(head));
73         enu=0;
74         get();
75         print();
76     }
77     return 0;
78 }

Bellman+剪枝：14ms水过；

 1 #include<stdio.h>
 2 #include<string.h>
 3 #define MIN(x,y)(x<y?x:y)
 4 const int INF=0x3f3f3f3f;
 5 const int MAXN=1010;
 6 const int MAXM=100010;
 7 int dis[MAXN],s[MAXN],e[MAXN],used[MAXN];
 8 struct Node{
 9     int s,e,dis;
10 };
11 Node dt[MAXM];
12 int T,S,D,top,ans;
13 void Bellman(){
14     int u,v,w;
15     memset(dis,INF,sizeof(dis));
16     for(int i=0;i<S;i++)dis[s[i]]=0;
17         for(int j=1;j<MAXN;j++){
18             if(!used[j])continue;
19     for(int i=0;i<top;i++){
20             u=dt[i].s;v=dt[i].e;w=dt[i].dis;
21             dis[v]=MIN(dis[v],dis[u]+w);
22         }
23     }ans=INF;
24     for(int i=0;i<D;i++)ans=MIN(dis[e[i]],ans);
25     printf("%d\n",ans);
26 }
27 int main(){
28     int a,b,c;
29     while(~scanf("%d%d%d",&T,&S,&D)){
30         top=0;
31         memset(used,0,sizeof(used));
32         while(T--){
33         scanf("%d%d%d",&a,&b,&c);
34         dt[top].s=a;dt[top].e=b;dt[top++].dis=c;
35         dt[top].s=b;dt[top].e=a;dt[top++].dis=c;
36         used[a]++;used[b]++;
37         }
38         for(int i=0;i<S;i++)scanf("%d",&s[i]);
39         for(int i=0;i<D;i++)scanf("%d",&e[i]);
40         Bellman();
41     }
42     return 0;
43 }