题目描述 Description
A 国有 n 座城市,编号从 1 到 n,城市之间有 m 条双向道路。每一条道路对车辆都有重量限制,简称限重。现在有 q 辆货车在运输货物,司机们
想知道每辆车在不超过车辆限重的情况下,最多能运多重的货物。
输入描述 Input Description
第一行有两个用一个空格隔开的整数 n,m,表示 A 国有 n 座城市和 m 条道路。
接下来 m 行每行 3 个整数 x、y、z,每两个整数之间用一个空格隔开,表示从 x 号城市到 y 号城市有一条限重为 z 的道路。注意:x 不等于 y,两座城市之间可能有多条道路。
接下来一行有一个整数 q,表示有 q 辆货车需要运货。
接下来 q 行,每行两个整数 x、y,之间用一个空格隔开,表示一辆货车需要从 x 城市运输货物到 y 城市,注意:x 不等于 y。
输出描述 Output Description
输出共有 q 行,每行一个整数,表示对于每一辆货车,它的最大载重是多少。如果货车不能到达目的地,输出-1。
样例输入 Sample Input
4 3
1 2 4
2 3 3
3 1 1
3
1 3
1 4
1 3
样例输出 Sample Output
3
-1
3
数据范围及提示 Data Size & Hint
对于 30%的数据,0 < n < 1,000,0 < m < 10,000,0 < q < 1,000;
对于 60%的数据,0 < n < 1,000,0 < m < 50,000,0 < q < 1,000;
对于 100%的数据,0 < n < 10,000,0 < m < 50,000,0 < q < 30,000,0 ≤ z ≤ 100,000。
算法的话应该很明显了:生成树加树上倍增,也没必要去纠结最大最小,生成树瓶颈树。
反正就是按边从大到小排序,之后搜生成树,生成完毕就break。记得每个点Dfs一次,
之后就直接倍增跑LCA。把每次sum(a,b)换成max(a,b)就行。
#include <cstdio>
#include <algorithm>
using namespace std;
const int maxx = 50000 + 100;
const int Inf = 100000 + 500;
template<class T> T chkmin(const T &a,const T &b){return a<b? a:b;}
int head[maxx],depth[maxx],father[maxx];
int grand[maxx][20+2],dis[maxx][20+2];
bool done[maxx];
int n,m,x,y,z,num,root,cnt,t;
inline int Read(){
int x=0,f=1;char c=getchar();
while(c>'9'||c<'0') {if(c=='-') f=-1;c=getchar();}
while(c>='0'&&c<='9') {x=x*10+c-'0'; c=getchar();}
return x*f;
}
struct Edge{
int next;
int to;
int value;
}Edges[maxx<<1];
struct Road{
int x;
int y;
int v;
}Roads[maxx<<1];
int cmp(Road a,Road b){
return a.v>b.v;
}
int find(int x){
if(father[x]!=x)
father[x]=find(father[x]);
return father[x];
}
void unionn(int x,int y){
int r1=find(x),r2=find(y);
if(r1!=r2)
father[r1]=r2;
}
void Add(int x,int y,int z){
Edges[++num].to=y;
Edges[num].next=head[x];
Edges[num].value=z;
head[x]=num;
}
void Dfs(int x){
done[x]=true;
for(int i=1;i<=20;i++){
if((1<<i) > depth[x]) break;
grand[x][i]=grand[grand[x][i-1]][i-1];
dis[x][i]=chkmin(dis[x][i-1],dis[grand[x][i-1]][i-1]);
}
for(int i=head[x];i;i=Edges[i].next){
int now=Edges[i].to;
if(done[now]) continue;
depth[now]=depth[x]+1;
grand[now][0]=x;
dis[now][0]=Edges[i].value;
Dfs(now);
}
}
inline int Lca(int x,int y){
int Ans=Inf;
if(depth[x]>depth[y]) swap(x,y);
int d=depth[y]-depth[x];
for(int i=0;i<=20;i++){
if((1<<i) & d) {
Ans = chkmin(Ans,dis[y][i]);
y=grand[y][i];
}
}
for(int i=20;i>=0;i--){
if(grand[x][i]!=grand[y][i]){
Ans = chkmin(Ans,dis[x][i]);
Ans = chkmin(Ans,dis[y][i]);
x=grand[x][i];
y=grand[y][i];
}
}
if(x==y) return Ans;
else{
Ans = chkmin(Ans,dis[x][0]);
Ans = chkmin(Ans,dis[y][0]);
return Ans;
}
}
int main(){
n=Read();m=Read();
for(int i=1;i<=n;i++)
father[i]=i;
for(int i=1;i<=m;i++)
Roads[i].x=Read(),Roads[i].y=Read(),Roads[i].v=Read();
sort(Roads+1,Roads+m+1,cmp);
for(int i=1;i<=m;i++){
if(find(Roads[i].x)!=find(Roads[i].y)){
Add(Roads[i].x,Roads[i].y,Roads[i].v);
Add(Roads[i].y,Roads[i].x,Roads[i].v);
unionn(Roads[i].x,Roads[i].y);
cnt++;
}
if(cnt==n-1)
break;
}
for(int i=1;i<=n;i++)
Dfs(i);
t = Read();
while( t-- ){
x=Read();y=Read();
if(find(x)!=find(y))
printf("-1\n");
else
printf("%d\n",Lca(x,y));
}
return 0;
}