BZOJ3723 Network——Kruskal重构树

Description

　　给你N个点的无向图 (1 <= N <= 15,000)，记为：1…N。 
　　图中有M条边 (1 <= M <= 30,000) ，第j条边的长度为： d_j ( 1 < = d_j < = 1,000,000,000).
　　现在有 K个询问 (1 < = K < = 15,000)。 
　　每个询问的格式是：A B，表示询问从A点走到B点的所有路径中，最长的边最小值是多少？

Input

　　第一行： N, M, K。 
　　第2..M+1行: 三个正整数：X, Y, and D (1 <= X <=N; 1 <= Y <= N). 表示X与Y之间有一条长度为D的边。 
　　第M+2..M+K+1行: 每行两个整数A B,表示询问从A点走到B点的所有路径中，最长的边最小值是多少？

Output

　　对每个询问，输出最长的边最小值是多少。

Sample Input

6 6 8

1 2 5

2 3 4

3 4 3

1 4 8

2 5 7

4 6 2

1 2

1 3

1 4

2 3

2 4

5 1

6 2

6 1

Sample Output

5

5

5

4

4

7

 4

5

Hint

1 <= N <= 15,000 
1 <= M <= 30,000 
1 <= d_j <= 1,000,000,000 
1 <= K <= 15,000

和noip2013d1t3类似，最小生成树+树剖可过，不过最近几场考试都考了重构树，感觉要跟上时代的潮流才行。。。

大致过程：

按照生成树的做法，每次选两个点出来之后，新建一个节点（权值为v），然后将这两个节点所在的两棵子树放在新节点的左右儿子。

整个过程用并查集维护一下。

如果还不懂，看这个dalao的blog：看图

#include<bits/stdc++.h>
using namespace std;
#define Inc(i,L,r) for(register int i=(L);i<=(r);++i)
#define Red(i,r,L) for(register int i=(r);i>=(L);--i)
const int N = 3e4+10;
struct Edge{
	int cnt,h[N],to[N<<1],next[N<<1];
	inline void add(int x,int y){
		next[++cnt]=h[x];
		to[cnt]=y;
		h[x]=cnt;
	}
}e;
struct edge{
	int x,y,w;
	bool operator <(const edge&A) const {
		return w<A.w;
	}
}a[N];
int n,m,k;
inline void init(){
	scanf("%d%d%d",&n,&m,&k);
	Inc(i,1,m){
		int x,y,z;scanf("%d%d%d",&x,&y,&z);
		a[i]=(edge){x,y,z};
	}
}
int fa[N],vl[N];
inline int getfa(int x){
	return (x^fa[x])?fa[x]=getfa(fa[x]):x;
}
inline void rebuild(){
	Inc(i,1,n<<1)fa[i]=i;
	sort(a+1,a+m+1);
	Inc(i,1,m){
		int fx=getfa(a[i].x),fy=getfa(a[i].y);
		if(fx==fy)continue;
		vl[++n]=a[i].w;
		fa[fx]=n,fa[fy]=n;
		e.add(n,fx),e.add(fx,n);
		e.add(n,fy),e.add(fy,n);
	}
}
int dep[N],p[N][17];
inline void dfs(int x,int fa,int depth){
	for(int P=e.h[x];P;P=e.next[P]){
		int to=e.to[P];
		if(to==fa)continue;
		dfs(to,p[to][0]=x,dep[to]=depth+1);
	}
}
inline void ST(){
	memset(p,-1,sizeof(p));
	dfs(n,0,0);
	int maxdep=log2(n);
	Inc(j,1,maxdep)
		Inc(i,1,n)if(~p[i][j-1])p[i][j]=p[p[i][j-1]][j-1];
}
inline int lca(int x,int y){
	if(dep[x]<dep[y])swap(x,y);
	Red(i,log2(dep[x]),0)if(dep[x]-(1<<i)>=dep[y])x=p[x][i];
	if(x==y)return y;
	Red(i,log2(dep[x]),0)if(p[x][i]^p[y][i])x=p[x][i],y=p[y][i];
	return p[x][0];
}
inline void solv(){
	Inc(i,1,k){
		int x,y;scanf("%d%d",&x,&y);
		cout<<vl[lca(x,y)]<<"\n";
	}
}
int main(){
	init();
	rebuild();
	ST();
	solv();
	return 0;
}