LCA+最大生成树 货车运输（洛谷 P1967）

货车运输

题目描述

A 国有 n 座城市，编号从 1 到 n，城市之间有 m 条双向道路。每一条道路对车辆都有重量限制，简称限重。

现在有 q 辆货车在运输货物， 司机们想知道每辆车在不超过车辆限重的情况下，最多能运多重的货物。

输入格式

第一行有两个用一个空格隔开的整数 n,m，表示 A 国有 n 座城市和 m 条道路。

接下来 m 行每行三个整数 x,y,z，每两个整数之间用一个空格隔开，表示从 x 号城市到 y 号城市有一条限重为 z 的道路。

注意： x$\ne$y 两座城市之间可能有多条道路 。

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

接下来 q 行，每行两个整数 x,y，之间用一个空格隔开，表示一辆货车需要从 x 城市运输货物到 y 城市，保证 x$\ne$y

输出格式

共有 q 行，每行一个整数，表示对于每一辆货车，它的最大载重是多少。
如果货车不能到达目的地，输出 -1。

---

这题可以转化为求两点之间的最小最大路径；

先转化为最大生成树，个人认为如果你不知道lca可以维护两点之间最小路径，那么是不会想到转化为最大生成树的；我就不知道；

之前学的lca只知道可以求两点之间的距离（树上就是最小距离），可以求公共祖先；没有做过这种题目；

自然不可能想到；

如果知道了这个，那么转化为最大生成树理所当然，这里讲一下怎么维护两点之间的最小路径；

首先知道lca的倍增思路也是一个区间思路，f[i][j]表示 i 结点往上走 2^j 个结点的结点；

那么是不是可以设一个mx[i][j]，表示 i 结点到 i-2^j 结点的路径最小值呢？这是不是就是一个区间思路，而且这个总路径上面区间不可能有重复；

代码：

#include<bits/stdc++.h>
#define LL long long
#define pa pair<int,LL>
#define ls k<<1
#define rs k<<1|1
#define inf 0x3f3f3f3f
using namespace std;
const int N=10100;
const int M=50100;
const int mod=1e9;
int n,m,q,head[N],cnt1,cnt2,lg[N],dep[N],f[N][40],mx[N][40],fa[N];
set<int>se;//存根 
int find(int p){
	if(p==fa[p]) return p;
	return fa[p]=find(fa[p]); 
}
struct Node{int p,q,w;}edge1[M*2];
struct Nod{int to,nex,w;}edge2[M*2];
void add1(int p,int q,int w){edge1[cnt1].w=w,edge1[cnt1].p=p,edge1[cnt1].q=q,cnt1++;}
void add2(int p,int q,int w){edge2[cnt2].w=w,edge2[cnt2].to=q,edge2[cnt2].nex=head[p],head[p]=cnt2++;} 
bool cmp(Node p,Node q){return p.w>q.w;}
void krus(){
	for(int i=1;i<=n;i++) fa[i]=i;
	sort(edge1+1,edge1+1+cnt1,cmp);
	for(int i=1;i<=cnt1;i++){
		int u=find(edge1[i].p);
		int v=find(edge1[i].q);
		if(u==v) continue;
		fa[u]=v;
		add2(edge1[i].p,edge1[i].q,edge1[i].w),add2(edge1[i].q,edge1[i].p,edge1[i].w);
	}
}
void dfs1(int p,int q){
	if(q==0) mx[p][0]=inf; 
	for(int i=head[p];~i;i=edge2[i].nex){
		if(edge2[i].to!=q){
			mx[edge2[i].to][0]=edge2[i].w;//主要是赋这个初值 
			dfs1(edge2[i].to,p);
		}
	}
} 
void dfs2(int p,int q){
	dep[p]=dep[q]+1,f[p][0]=q;
	for(int i=1;i<=lg[dep[p]];i++){
		mx[p][i]=min(mx[p][i-1],mx[f[p][i-1]][i-1]);
		f[p][i]=f[f[p][i-1]][i-1];
	}
	for(int i=head[p];~i;i=edge2[i].nex){
		if(edge2[i].to!=q) dfs2(edge2[i].to,p);
	}
}
int lca(int x,int y){
	if(find(x)!=find(y)) return -1;
	int ans=inf;
	if(dep[x]<dep[y]) swap(x,y);
	while(dep[x]>dep[y]) ans=min(ans,mx[x][lg[dep[x]-dep[y]]-1]),x=f[x][lg[dep[x]-dep[y]]-1];
	if(x==y) return ans; 
	for(int k=lg[dep[x]]-1;k>=0;k--){
		if(f[x][k]!=f[y][k]){
			ans=min(ans,min(mx[x][k],mx[y][k]));
			x=f[x][k],y=f[y][k];
		}
	}

	ans=min(ans,min(mx[x][0],mx[y][0]));
	return ans;
}
int main(){
	memset(head,-1,sizeof(head));
	memset(mx,inf,sizeof(mx));
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++){
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		add1(x,y,z),add1(y,x,z);
	}
	krus();
	for(int i=1;i<=n;i++) lg[i]=lg[i-1]+(1<<lg[i-1]==i);
	for(int i=1;i<=n;i++){
		if(se.count(find(i))) continue;
		se.insert(find(i));
		dfs1(i,0);//lca初始化 
		dfs2(i,0);//lca初始化 
	}
	scanf("%d",&q);
	while(q--){
		int x,y;
		scanf("%d%d",&x,&y);
		printf("%d\n",lca(x,y));
	}
	return 0;
}