BZOJ1050 BSOJ2016 【HAOI2006】旅行 CODEVS 1001 舒适的路线

2016 -- 【HAOI2006】旅行 CODEVS 1001 舒适的路线

Description

　　Z小镇是一个景色宜人的地方，吸引来自各地的观光客来此旅游观光。Z小镇附近共有N个景点（编号为1,2,3,…,N），这些景点被M条道路连接着，所有道路都是双向的，两个景点之间可能有多条道路。也许是为了保护该地的旅游资源，Z小镇有个奇怪的规定，就是对于一条给定的公路Ri，任何在该公路上行驶的车辆速度必须为Vi。速度变化太快使得游客们很不舒服，因此从一个景点前往另一个景点的时候，大家都希望选择行使过程中最大速度和最小速度的比尽可能小的路线，也就是所谓最舒适的路线。

Input

　　第一行包含两个正整数，N和M。

　　接下来的M行每行包含三个正整数：x，y和v。表示景点x到景点y之间有一条双向公路，车辆必须以速度v在该公路上行驶。

　　最后一行包含两个正整数s，t，表示想知道从景点s到景点t最大最小速度比最小的路径。s和t不可能相同。

Output

　　如果景点s到景点t没有路径，输出“IMPOSSIBLE”。否则输出一个数，表示最小的速度比。如果需要，输出一个既约分数。

Sample Input

4 2

1 2 1

3 4 2

1 4

Sample Output

IMPOSSIBLE

Sample Input

　　3 3
　　1 2 10
　　1 2 5
　　2 3 8
　　1 3

Sample Output

　　5/4

Sample Output

　　3 2
　　1 2 2
　　2 3 4

　　1 3

Sample Output

　　2

　　数据范围
　　　　1 < N ≤ 500
　　　　1 ≤ x,y ≤ N，0 < v < 30000，x≠y

　　　　0 < M ≤ 5000

这道题一开始没什么想法，后来别人提醒才知道是生成树，我们根据生成树的性质，枚举起点即可，跟刘汝佳书上的一道题很像，没想起来实在不应该。

而且在编程的时候也犯了一些小错误，比如说循环变量乱用，还有在判断 s 和 t 是否联通时，我们不能直接用fa[s]==fa[t],而应使用getfa(s)==getfa(t)，因为我们在加入一条边之后，可能fa[s]和fa[t]还没有更新，指向的fa还是之前的，还是对并查集理解不深刻啊。

代码如下，还是比较简单

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdio>
#define eps 0.000001
using namespace std;
struct nodee{int fr,to,val;
}e[10005];
int maxn=0x3fffffff,minn=1,maxt=0,mint=0x3fffffff,m,n,fa[100005],s,t;
bool cmp(nodee a,nodee b)
{
	return a.val<b.val;
}
int gcd(int a,int b)
{
	if(b==0)return a;
	return gcd(b,a%b);
}
void check()
{
	double b1=(maxn*1.0)/(minn*1.0);
	double b2=(maxt*1.0)/(mint*1.0);
	if(b2-b1<eps)maxn=maxt,minn=mint;
	return ;
}
int getfa(int x)
{
	if(x==fa[x])return x;
	return fa[x]=getfa(fa[x]);
}
void K()
{
	sort(e+1,e+m+1,cmp);
	int bj=0;
	for(int i=1;i<=m;i++)
	{
		int c1=0;
		for(int j=1;j<=n;j++)
	      fa[j]=j;
	    maxt=e[i].val,mint=e[i].val;
		for(int j=i;j<=m;j++)
		{
			int f1=getfa(e[j].fr),f2=getfa(e[j].to);
			if(f1!=f2)
			{
				fa[f1]=f2;
				getfa(e[j].to);getfa(e[j].fr);
				mint=min(mint,e[j].val);
				maxt=max(maxt,e[j].val);
				c1++;
			}
			if(getfa(s)==getfa(t))break;
		}
		if(fa[s]==fa[t])
		{
			check();
			bj=1;
		}
	}
	if(bj)
	{
		int G=gcd(minn,maxn);
		if(minn/G!=1)
		  cout<<maxn/G<<"/"<<minn/G;
		else cout<<maxn/G;
	}
	else cout<<"IMPOSSIBLE\n";
}
int main(){
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	  cin>>e[i].fr>>e[i].to>>e[i].val;
	cin>>s>>t;
	K();
	return 0;
}