一个人的旅行(floyd+dijskra+SPFA+Bellman)
一个人的旅行
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个数,表示草儿想去地方。
接着有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
题解: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 }