洛谷P1967 货车运输

题目描述

A 国有 n 座城市，编号从 1 到 n，城市之间有 m 条双向道路。每一条道路对车辆都有重量限制，简称限重。现在有 q 辆货车在运输货物， 司机们想知道每辆车在不超过车辆限重的情况下，最多能运多重的货物。

输入输出格式

输入格式：

输入文件名为 truck.in。

输入文件第一行有两个用一个空格隔开的整数 n，m，表示 A 国有 n 座城市和 m 条道

路。 接下来 m 行每行 3 个整数 x、 y、 z，每两个整数之间用一个空格隔开，表示从 x 号城市到 y 号城市有一条限重为 z 的道路。注意： x 不等于 y，两座城市之间可能有多条道路 。

接下来一行有一个整数 q，表示有 q 辆货车需要运货。

接下来 q 行，每行两个整数 x、y，之间用一个空格隔开，表示一辆货车需要从 x 城市运输货物到 y 城市，注意： x 不等于 y 。

输出格式：

输出文件名为 truck.out。

输出共有 q 行，每行一个整数，表示对于每一辆货车，它的最大载重是多少。如果货

车不能到达目的地，输出-1。

输入输出样例

输入样例#1： 

4 3
1 2 4
2 3 3
3 1 1
3
1 3
1 4
1 3

输出样例#1： 

3
-1
3

说明

对于 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。

最大生成树+LCA+暴搜
因为要树上暴搜，所以LCA用倍增。。。
附代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<queue>
#define MAXN 10010
#define MAXM 50010
#define MAX 999999999
using namespace std;
int n,m,q,c=1;
int head[MAXN],fa[MAXN],deep[MAXN],f[MAXN][20];
struct node1{
	int next,to,w;
}a[MAXM<<1];
struct node2{
	int u,v,w;
}b[MAXM];
inline int read(){
	int date=0,w=1;char c=0;
	while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}
	while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();}
	return date*w;
}
bool cmp(const node2 &x,const node2 &y){
	return x.w>y.w;
}
int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
void uniun(int x,int y){x=find(x);y=find(y);if(x!=y)fa[y]=x;}
void add(int u,int v,int w){
	a[c].to=v;a[c].w=w;
	a[c].next=head[u];
	head[u]=c++;
	a[c].to=u;a[c].w=w;
	a[c].next=head[v];
	head[v]=c++;
}
void kruskal(){
	for(int i=1;i<=m;i++)
	if(find(b[i].u)!=find(b[i].v)){
		uniun(b[i].u,b[i].v);
		add(b[i].u,b[i].v,b[i].w);
	}
}
void buildtree(int rt){
	int will;
	for(int i=head[rt];i;i=a[i].next){
		will=a[i].to;
		if(!deep[will]){
			deep[will]=deep[rt]+1;
			f[will][0]=rt;
			buildtree(will);
		}
	}
}
void step(){
	for(int i=1;i<=19;i++)
	for(int j=1;j<=n;j++)
	f[j][i]=f[f[j][i-1]][i-1];
}
int LCA(int x,int y){
	if(deep[x]<deep[y])swap(x,y);
	for(int i=19;i>=0;i--)
	if(deep[f[x][i]]>=deep[y])
	x=f[x][i];
	if(x==y)return x;
	for(int i=19;i>=0;i--)
	if(f[x][i]!=f[y][i]){
		x=f[x][i];
		y=f[y][i];
	}
	return f[x][0];
}
void solve(int x,int y){
	int fa=LCA(x,y),s1=MAX,s2=MAX,ans;
	if(fa==0){
		printf("-1\n");
		return;
	}
	for(int i=x;i!=fa;i=f[i][0]){
		int path;
		for(int j=head[i];j;j=a[j].next)if(a[j].to==f[i][0])path=a[j].w;
		s1=min(s1,path);
	}
	for(int i=y;i!=fa;i=f[i][0]){
		int path;
		for(int j=head[i];j;j=a[j].next)if(a[j].to==f[i][0])path=a[j].w;
		s2=min(s2,path);
	}
	ans=min(s1,s2);
	printf("%d\n",ans);
}
int main(){
	int u,v,w;
	n=read();m=read();
	for(int i=1;i<=n;i++)fa[i]=i;
	for(int i=1;i<=m;i++){b[i].u=read();b[i].v=read();b[i].w=read();}
	sort(b+1,b+m+1,cmp);
	kruskal();
	deep[1]=1;
	buildtree(1);
	step();
	q=read();
	while(q--){
		u=read();v=read();
		solve(u,v);
	}
	return 0;
}