【洛谷】1967 [noip2013]货车运输 最小生成树+LCA

题目传送门

这题和上一题的解题思路是一模一样的，也是先对给出的边做一遍最小生成树，哦不对，这题是最大生成树。

然后还是把无根树转化成有根树，求每个点的倍增祖先。

然后，这题的判断有点特殊，其实也挺简单的，可能是我为了中考失了智吧……

就是判断给出的两个节点的祖先是否相同，若不同，则输出-1；若相同，则输出这两个节点到它们LCA的路径上的最小边权。

其他只要想进去就挺简单的，总之细心一点吧。

附上AC代码：

#include <cstdio>
#include <cctype>
#include <algorithm>
#include <cmath>
#define N 10010
using namespace std;

struct note{
	int x,y,w;
	bool operator < (const note lyf) const {return w>lyf.w;}
}a[50010];
struct side{
	int to,w,nt;
}s[N<<1];
int n,m,q,fa[N],fx,fy,h[N],num,size[N],mx[N],rt,d[N],f[N][15],k,f_min[N][15],x,y;

inline char nc(){
	static char ch[100010],*p1=ch,*p2=ch;
	return p1==p2&&(p2=(p1=ch)+fread(ch,1,100010,stdin),p1==p2)?EOF:*p1++;
}

inline void read(int &a){
	static char c=nc();int f=1;
	for (;!isdigit(c);c=nc()) if (c=='-') f=-1;
	for (a=0;isdigit(c);a=a*10+c-'0',c=nc());
	a*=f;return;
}

inline int gf(int x){
	return !fa[x]||x==fa[x]?x:fa[x]=gf(fa[x]);
}

inline void add(int x,int y,int w){
	s[++num]=(side){y,w,h[x]},h[x]=num;
	s[++num]=(side){x,w,h[y]},h[y]=num;
}

inline void so(int x,int fa){
	size[x]=1,mx[x]=0;
	for (int i=h[x]; i; i=s[i].nt)
		if (s[i].to!=fa)
			so(s[i].to,x),size[x]+=size[s[i].to],mx[x]=min(mx[x],size[s[i].to]);
	mx[x]=min(mx[x],n-size[x]);
	if (mx[x]<mx[rt]) rt=x;
	return;
}

inline void build(int x){
	d[x]=d[f[x][0]]+1;
	for (int i=h[x]; i; i=s[i].nt)
		if (s[i].to!=f[x][0])
			f[s[i].to][0]=x,f_min[s[i].to][0]=s[i].w,build(s[i].to);
	return;
}

inline int query(int x,int y){
	if (d[x]<d[y]) swap(x,y);
	int ans=2e9;
	for (int i=k; i>=0; --i)
		if (d[f[x][i]]>=d[y]) ans=min(ans,f_min[x][i]),x=f[x][i];
	if (x==y) return ans;
	for (int i=k; i>=0; --i)
		if (f[x][i]!=f[y][i]){
			ans=min(ans,min(f_min[x][i],f_min[y][i]));
			x=f[x][i],y=f[y][i];
		}
	ans=min(ans,min(f_min[x][0],f_min[y][0]));
	return ans;
}

int main(void){
	read(n),read(m);
	for (int i=1; i<=m; ++i) read(a[i].x),read(a[i].y),read(a[i].w);
	sort(a+1,a+1+m);
	for (int i=1; i<=m; ++i){
		fx=gf(a[i].x),fy=gf(a[i].y);
		if (fx!=fy) fa[fx]=fy,add(a[i].x,a[i].y,a[i].w);
	}
	mx[0]=n,so(1,0),build(rt);
	for (int j=1; j<=(k=log(n)/log(2)+1); ++j)
		for (int i=1; i<=n; ++i)
			f[i][j]=f[f[i][j-1]][j-1],f_min[i][j]=min(f_min[i][j-1],f_min[f[i][j-1]][j-1]);
	read(q);
	while (q--){
		read(x),read(y);
		printf("%d\n",gf(x)!=gf(y)?-1:query(x,y)); 
	}
	return 0;
}