树上距离HDU2586（在线RMQ算法）

题干请参考
树上距离HDU2586（LCA）
分析：
两个节点的LCA一定是两个节点之间欧拉序列中深度最小的节点，寻找深度最小值是可以使用RMQ(区间最值）算法。

欧拉序列就是指在深度遍历过程中依次经过的节点记录下来，把回溯时经过的节点也记录下来，一个节点可能被记录多次，相当于从树根开始，一笔画出一个经过所有节点的回路。
如下图：

该树的欧拉序列为1 2 4 6 8 6 9 6 4 2 5 7 5 2 1 3 1，搜索时得到6和5首次出现的下标i、j，然后查询该区间深度最小的节点，为6和5号节点的最近公共祖先。

算法实现：
1、深度遍历，得到三个数组，每个节点首次出现的下标pos[ ],深度遍历得到的欧拉序列seq[ ],每个节点的深度dep[ ]。
2、根据欧拉序列的深度，创建ST。
3、查询[l,r]区间最小深度，最小深度节点就是LCA。
4、求两个节点之间的距离。
代码如下：

#include<iostream>
#include<algorithm>
#include<cmath>//使用log函数
using namespace std;
const int maxn=40005;
int n,m,tot,F[maxn<<1][20];//F(i,j)表示区间[i，i+2^j-1]的最值，区间长度为2^j
int head[maxn],dist[maxn],pos[maxn],seq[maxn<<1],dep[maxn<<1],cnt;//头结点
bool vis[maxn];
struct Edge{
	int to,c,next;
}e[maxn<<1];
 
void add(int u,int v,int w){
	e[++cnt].to=v;
	e[cnt].c=w;
	e[cnt].next=head[u];
	head[u]=cnt;
}

void dfs(int u,int d){//dfs序
	vis[u]=true;
	//表示u节点已访问，访问v时v-u边就不做处理，这样可以保证一条边只处理一次。
	//（以u、v节点的边被加了两次，当成u-v和v-u两条边） 
	pos[u]=++tot;//u首次出现的下标 
	seq[tot]=u;//dfs遍历序列
	dep[tot]=d;//深度
	for(int i=head[u];i;i=e[i].next){
		int v=e[i].to,w=e[i].c;
		if(vis[v])
			continue;
		dist[v]=dist[u]+w;
		dfs(v,d+1);
		seq[++tot]=u;//递归返回到u节点，dfs遍历序列
		dep[tot]=d;//返回到u节点的深度和经过u节点的深度是一样的 
	}
}/*(1)请初学的同学注意，递归机制每一次调用都会给这次调用保存所有变量的值
(2)在dfs(v,d+1)之前的语句，是在调用之前执行的语句。在 dfs(v,d+1)之后的语句，
是递归调用返回来执行的语句。
（3）在dfs(v,d+1)语句前后，相同变量的值都是相同的。 
*/ 

void ST_create(){
	for(int i=1;i<=tot;i++)//初始化 
		F[i][0]=i;//记录下标，不是最小深度 
	int k=log2(tot);
	for(int j=1;j<=k;j++)
		for(int i=1;i<=tot-(1<<j)+1;i++)//tot-2^j+1
			if(dep[F[i][j-1]]<dep[F[i+(1<<(j-1))][j-1]])
				F[i][j]=F[i][j-1];
			else
				F[i][j]=F[i+(1<<(j-1))][j-1];
}

int RMQ_query(int l,int r){//求区间[l..r]的最值
	int k=log2(r-l+1);
	if(dep[F[l][k]]<dep[F[r-(1<<k)+1][k]])
		return 	F[l][k];
	else
		return 	F[r-(1<<k)+1][k];//返回深度最小的结点下标 
}

int LCA(int x,int y){//求x，y的最近公共祖先
	int l=pos[x],r=pos[y];//读取第一次出现的下标
	if(l>r)
		swap(l,r);
	return seq[RMQ_query(l,r)];//返回结点 
} 

void print(){
	for(int i=1;i<=tot;i++)
		cout<<seq[i]<<" ";
	cout<<endl;
}

int main(){
	int x,y,T,lca;
	cin>>T;
	while(T--){
		cin>>n>>m;
		for(int i=1;i<=n;i++)//初始化 
			head[i]=vis[i]=dist[i]=0;
		cnt=0,tot=0;
		for(int i=1;i<n;i++){//输入一棵树的n-1边
			int x,y,z;
			cin>>x>>y>>z;
			add(x,y,z);
			add(y,x,z);
		}
		dfs(1,1);
		print();
		ST_create();//创建ST表
		for(int i=1;i<=m;i++){
			cin>>x>>y;
			lca=LCA(x,y);
			cout<<dist[x]+dist[y]-2*dist[lca]<<endl;//输出x y的距离 
		}	
	}
	return 0;
}
/*
1
9 8
1 3 2
1 2 3
2 5 1
2 4 6
4 6 2
5 7 5
6 9 7
6 8 8
输出 1 2 4 6 8 6 9 6 4 2 5 7 5 2 1 3 1
6 5
输出 9
*/