【bzoj1050】[HAOI2006]旅行comf

1050: [HAOI2006]旅行comf

Time Limit: 10 Sec   Memory Limit: 162 MB
Submit: 3376   Solved: 1875
[ Submit][ Status][ Discuss]

Description

给你一个无向图，N(N<=500)个顶点, M(M<=5000)条边，每条边有一个权值Vi(Vi<30000)。给你两个顶点S和T，求

一条路径，使得路径上最大边和最小边的比值最小。如果S和T之间没有路径，输出”IMPOSSIBLE”，否则输出这个

比值，如果需要，表示成一个既约分数。 备注： 两个顶点之间可能有多条路径。

Input

第一行包含两个正整数，N和M。下来的M行每行包含三个正整数：x，y和v。表示景点x到景点y之间有一条双向公路

，车辆必须以速度v在该公路上行驶。最后一行包含两个正整数s，t，表示想知道从景点s到景点t最大最小速度比

最小的路径。s和t不可能相同。

1<N<=500,1<=x,y<=N，0<v<30000，0<M<=5000

Output

如果景点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

HINT

Source

[ Submit][ Status][ Discuss]



啊。。昨天晚上写了这题早上模拟赛考了一个差不多的

这题我们考虑枚举最大值和最小值

先把所有边按权值从小到大排序排好

然后对于一个在排序后数组中编号为i的边，我们枚举一个编号为j的边，如果区间[i,j]中的边能使s,t连通，那么更新答案

然后剩下的就是并查集判连通的事啦

PS.这题的区间[i,j]的转移可以通过滑动窗口优化，然而并没有优化多少。

代码：

#include<cstdio>
#include<vector>
#include<queue>
#include<ctime>
#include<algorithm>
#include<cstdlib>
#include<stack>
#include<cstring>
#include<cmath>
using namespace std;

typedef long long LL;

const int INF = 2147483647;
const int maxm = 5010;
const int maxn = 510;

struct data{
	int u,v,w;
}e[maxm];

int fa[maxn];
int n,m,y,x,s,t;

inline LL getint()
{
	LL ret = 0,f = 1;
	char c = getchar();
	while (c < '0' || c > '9')
	{
		if (c == '-') f = -1;
		c = getchar();
	}
	while (c >= '0' && c <= '9')
		ret = ret * 10 + c - '0',c = getchar();
	return ret * f;
}

inline int gcd(int a,int b)
{
	if (!b) return a;
	else return gcd(b,a % b);
}

inline bool cmp(data a,data b)
{
	return a.w < b.w;
}

inline int find(int x)
{
	return fa[x] == x ? x : fa[x] = find(fa[x]);
}

int main()
{
	n = getint(); m = getint();
	for (int i = 1; i <= m; i++)
	{
		e[i].u = getint();
		e[i].v = getint();
		e[i].w = getint();
	}
	s = getint();
	t = getint();
	sort(e + 1,e + m + 1,cmp);
	y = 1;
	x = 40000;
	for (int i = 1; i <= m; i++)
	{
		if (e[i].w == e[i - 1].w) continue;
		for (int j = 1; j <= n; j++) fa[j] = j;
		for (int j = i; j <= m; j++)
		{
			if (find(e[j].u) != find(e[j].v)) fa[fa[e[j].v]] = fa[e[j].u];
			if (e[j].w == e[j + 1].w) continue;
			if (find(s) == find(t))
			{
				if (x * e[i].w > y * e[j].w)
				{
					x = e[j].w;
					y = e[i].w;
				}
				break;
			}
		}
	}
	if (x == 40000) printf("IMPOSSIBLE");
	else
	{
		int g = gcd(x,y);
		if (y == g) printf("%d",x / g);
		else printf("%d/%d",x / g,y / g);
	}
	return 0;
}