bzoj1050 旅行comf

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不可能相同。

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

【數據範圍】

1<  N < = 500

1 < = x, y < = N，0 < v < 30000，x ≠ y

0 < M < =5000

Source

枚舉最大邊，然後求最小生成樹，記錄最優值。

然後，就沒了。

代碼：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
inline int read()
{
	int x = 0, f = 1;
	char ch = getchar();
	while(ch < '0' || ch > '9')
	{
		if(ch == '-') f = -1;
		ch = getchar();
	}
	while(ch >= '0' && ch <= '9')
	{
		x = x * 10 + ch - '0';
		ch = getchar();
	}
	return x * f;
}
struct Edge
{
	int u, v, w;
}a[5010];
int n, m, f[510];
int gcd(int a, int b)
{
	if(a > b) swap(a, b);
	while(a)
	{
		int c = a;
		a = b % a;
		b = c;
	}
	return b;
}

int find(int x)
{
	if(f[x] == x) return f[x];
	else return f[x] = find(f[x]);
}

int cmp(Edge a, Edge b)
{
	return a.w < b.w;
}

int main()
{
	n = read();
	m = read();
	for(int i = 1; i <= m; i ++)
		a[i].u = read(), a[i].v = read(), a[i].w = read();
	sort(a + 1, a + m + 1, cmp);
	int s, t, start = 1, ansmx = 1, ansmn = 0;
	s = read();
	t = read();
	while(start <= m)
	{
		int mn, mx, x;
		mn = mx = -1;
		for(int i = 1; i <= n; i ++)
			f[i] = i;
		for(x = start; x <= m; x ++)
		{
			int u = find(a[x].u);
			int v = find(a[x].v);
			f[u] = v;
			if(find(s) == find(t))
			{
				mx = a[x].w;
				break;
			}
		}
		if(mx == -1)
		{
			if(!ansmn) {printf("IMPOSSIBLE\n"); return 0;}
			else break;
		}		
		for(int i = 1; i <= n; i ++)
			f[i] = i;
		for(; x > 0; x --)
		{
			int u = find(a[x].u);
			int v = find(a[x].v);
			f[u] = v;
			if(find(s) == find(t)) 
			{
				mn = a[x].w;
				break;
			}
		}
		start = x + 1;
		if(mn == -1)
		{
			if(!ansmn) {printf("IMPOSSIBLE\n"); return 0;}
			else break;
		}		
		int r = gcd(mx, mn);
		mx /= r;
		mn /= r;
		if(ansmx * mn > ansmn * mx)
		{
			ansmn = mn;
			ansmx = mx;
		}
	}
	if(ansmn == 1) printf("%d\n", ansmx);
	else printf("%d/%d\n", ansmx, ansmn);
	
	return 0;
}