BZOJ1050 [HAOI2006]旅行comf（Kruskal算法）

SPFA是错误的：局部最优不能保证全局最优，因为要求的是比例最小，与路的长短无直接关系

可以设计m^2的算法：
先将边按长度排序，然后每次枚举最小边的长度w[i]作限定，依次加入更大的边，
当加到某条边w[j]后，s与t连通(并查集判断)，就用w[j]/w[i](double类型)更新答案 
若没有w[i]能使s与t连通，则无解 

注意：路径中分出的"杈"不用管，因为它们的长度既非最小，也非最大，去掉后s,t仍连通，不会影响答案 

#include<stdio.h>
#include<stdlib.h>
int fa[10005]={0},u[10005]={0},v[10005]={0},w[10005]={0};
int father(int x)
{
	if(fa[x]!=x) fa[x]=father(fa[x]);
	return fa[x];
}
void jh(int* a,int* b)
{
	int t=*a;
	*a=*b;
	*b=t;
}
void kp(int low,int high)
{
	int i=low,j=high,mid=w[(i+j)/2];
	while(i<j)
	{
		while(w[i]<mid) i++;
		while(w[j]>mid) j--;
		if(i<=j)
		{
			jh(&u[i],&u[j]);
			jh(&v[i],&v[j]);
			jh(&w[i],&w[j]);
			i++;
			j--;
		}
	}
	if(j>low) kp(low,j);
	if(i<high) kp(i,high);
}
int gcd(int a,int b)
{
	if(b==0) return a;
	return gcd(b,a%b);
}
int main()
{
	int n,m,i,j,s,t,min=-1,max=-1;
	scanf("%d%d",&n,&m);
	for(i=1;i<=m;i++)
		scanf("%d%d%d",&u[i],&v[i],&w[i]);
	scanf("%d%d",&s,&t);
	kp(1,m);
	for(i=1;i<=m;i++)//枚举最小边权:w[i]
	{
		if(w[i]==w[i-1]) continue;//快得不止一点点 
		for(j=1;j<=n;j++)
			fa[j]=j;
		for(j=i;j<=m;j++)
		{
			if(father(u[j])!=father(v[j])) fa[father(u[j])]=father(v[j]);
			if(father(s)==father(t)) break;
		}
		if(j<=m&&(min==-1||(double)max/(double)min>(double)w[j]/(double)w[i]))
		{
			min=w[i];
			max=w[j];
		}
	}
	if(min==-1) printf("IMPOSSIBLE");
	else
	{
		if(max%min==0) printf("%d",max/min);
		else printf("%d/%d",max/gcd(max,min),min/gcd(max,min));
	}
	return 0;
}