codevs 1001 舒适的路线
2006年
时间限制: 2 s
空间限制: 128000 KB
题目等级 : 钻石 Diamond
题目描述 Description
Z小镇是一个景色宜人的地方，吸引来自各地的观光客来此旅游观光。
Z小镇附近共有
N(1<N≤500)个景点（编号为1,2,3,…,N），这些景点被M（0<M≤5000）条道路连接着，所有道路都是双向的，两个景点之间可能有多条道路。也许是为了保护该地的旅游资源，Z小镇有个奇怪的规定，就是对于一条给定的公路Ri，任何在该公路上行驶的车辆速度必须为Vi。频繁的改变速度使得游客们很不舒服，因此大家从一个景点前往另一个景点的时候，都希望选择行使过程中最大速度和最小速度的比尽可能小的路线，也就是所谓最舒适的路线。
输入描述 Input Description
第一行包含两个正整数，N和M。
接下来的M行每行包含三个正整数：x，y和v（1≤x,y≤N，0 最后一行包含两个正整数s，t，表示想知道从景点s到景点t最大最小速度比最小的路径。s和t不可能相同。
输出描述 Output Description
如果景点s到景点t没有路径，输出“IMPOSSIBLE”。否则输出一个数，表示最小的速度比。如果需要，输出一个既约分数。
样例输入 Sample Input
样例1
4 2
1 2 1
3 4 2
1 4
样例2
3 3
1 2 10
1 2 5
2 3 8
1 3
样例3
3 2
1 2 2
2 3 4
1 3
样例输出 Sample Output
样例1
IMPOSSIBLE
样例2
5/4
样例3
2
数据范围及提示 Data Size & Hint
N(1<N≤500)
M（0<M≤5000）
Vi在int范围内
最最基本的思想：枚举最小边，把和它组合最好的最大边找出来，这样找到所有的最小最大值组合，从其中找到比例最小的。
网上的题解+注释：
1 /*
2 其实这个题目看似是类似于最短路的问题，但实际上与最短路无关。
3 因为路是双向的，所以我们为了最优的比例，当然可以去走一条无关的路。
4 可以设计m^2的算法：
5 先将边按长度排序，然后每次枚举最小边的长度w[i]作限定，依次加入更大的边，
6 当加到某条边w[j]后，s与t连通(并查集判断)，就用w[j]/w[i](double类型)更新答案
7 若没有w[i]能使s与t连通，则无解
8
9 注意：路径中分出的"杈"不用管，因为即使无关，我们也可以走（仅仅是为了那个比例），去掉后s,t仍连通，不会影响答案
10 */
11 #include<stdio.h>
12 #include<stdlib.h>
13 int fa[10005]={0},u[10005]={0},v[10005]={0},w[10005]={0};
14 int father(int x)
15 {
16 if(fa[x]!=x) fa[x]=father(fa[x]);
17 return fa[x];
18 }
19 void jh(int* a,int* b)
20 {
21 int t=*a;
22 *a=*b;
23 *b=t;
24 }
25 void kp(int low,int high)
26 {
27 int i=low,j=high,mid=w[(i+j)/2];
28 while(i<j)
29 {
30 while(w[i]<mid) i++;
31 while(w[j]>mid) j--;
32 if(i<=j)
33 {
34 jh(&u[i],&u[j]);
35 jh(&v[i],&v[j]);
36 jh(&w[i],&w[j]);
37 i++;
38 j--;
39 }
40 }
41 if(j>low) kp(low,j);
42 if(i<high) kp(i,high);
43 }
44 int gcd(int a,int b)
45 {
46 if(b==0) return a;
47 return gcd(b,a%b);
48 }
49 int main()
50 {
51 int n,m,i,j,s,t,min=-1,max=-1;
52 scanf("%d%d",&n,&m);
53 for(i=1;i<=m;i++)
54 scanf("%d%d%d",&u[i],&v[i],&w[i]);
55 scanf("%d%d",&s,&t);
56 kp(1,m);
57 for(i=1;i<=m;i++)//枚举最小边权:w[i]
58 {
59 if(w[i]==w[i-1]) continue;/*快得不止一点点，这个剪枝，之所以正确，就是因为可以走无关的路，所以最小边的大小不变的话，再建一棵树，结果还是相同的。*/
60 for(j=1;j<=n;j++)
61 fa[j]=j; /*每次都要建树，所以每次都把并查集初始化*/
62 for(j=i;j<=m;j++)
63 {
64 if(father(u[j])!=father(v[j])) fa[father(u[j])]=father(v[j]);
65 if(father(s)==father(t)) break; /*用并查集检验s-t是否联通来找的最长边，然后最短边是枚举的，那么就可以很轻易地求出比值了*/
66 }
67 if(j<=m&&(min==-1||(double)max/(double)min>(double)w[j]/(double)w[i])) /*比值用double，更准*/
68 {
69 min=w[i];
70 max=w[j];
71 }
72 }
73 if(min==-1) printf("IMPOSSIBLE");
74 else
75 {
76 if(max%min==0) printf("%d",max/min);
77 else printf("%d/%d",max/gcd(max,min),min/gcd(max,min));
78 }
79 return 0;
80 }
my代码：
1 #include<iostream>
2 using namespace std;
3 #include<cstdio>
4 #include<algorithm>
5 #define N 505
6 #define M 5005
7 struct Edge{
8 int u,v,w;
9 bool operator<(Edge P)
10 const{return w<P.w;}
11 }edge[M];
12 int n,m,s,t,maxx=-1,minn=-1,father[N];
13 int gcd(int a,int b)
14 {
15 if(!b) return a;
16 return gcd(b,a%b);
17 }
18 void input()
19 {
20 scanf("%d%d",&n,&m);
21 for(int i=1;i<=m;++i)
22 {
23 scanf("%d%d%d",&edge[i].u,&edge[i].v,&edge[i].w);
24 }
25 scanf("%d%d",&s,&t);
26 }
27 int find(int x)
28 {
29 if(father[x]==x) return x;
30 return father[x]=find(father[x]);
31 }
32 bool kruskal()
33 {
34 sort(edge+1,edge+m+1);
35 for(int i=1;i<=m;++i)
36 {
37 if(edge[i].w==edge[i-1].w) continue;
38 int j;
39 for(j=1;j<=n;++j)
40 father[j]=j;
41 for(j=i;j<=m;++j)/*j循环要从i开始，因为此时i这条边还没有加入到生成树中*/
42 {
43 int x1=find(edge[j].u),x2=find(edge[j].v);
44 if(x1!=x2)
45 {
46 father[x2]=x1;
47 }
48 if(find(s)==find(t)) break;
49 }
50 if(j<=m&&(minn==-1||((double)maxx/(double)minn>(double)edge[j].w/(double)edge[i].w)))
51 {
52 minn=edge[i].w;
53 maxx=edge[j].w;
54 }
55 }
56 if(minn==-1) return false;
57 return true;
58 }
59 int main()
60 {
61 input();
62 if(!kruskal())
63 printf("IMPOSSIBLE\n");
64 else {
65 if(maxx%minn==0) printf("%d",maxx/minn);
66 else printf("%d/%d",maxx/gcd(maxx,minn),minn/gcd(maxx,minn));
67 }
68 return 0;
69 }
时间： 2024-12-23 15:30:32