倍增法求LCA

一般用于求解树上两结点距离。LCA是两节点的最近公共祖先。

倍增法求解LCA实际上是将上升的步数用二进制来表示，经过向上跳两结点深度的比较来确定是否可以上升，最终找到LCA。

先建树维护几个数组

father[][];   记录祖先结点编号 
deep[];       记录结点深度
dis[][]记录结点和祖先结点的距离 

mxdeep 代表最大层数   mxdeep=log(n)/log(2);

代码：

void dfs(int root)
{
	for(int i=1;i<=mxdeep;i++)
	{
		father[root][i]=father[father[root][i-1]][i-1];
		dis[root][i]=dis[father[root][i-1]][i-1]+dis[root][i-1];
		if(!father[root][i])   //优化，当为0的时候说明已经更新到了根结点位置 
		break;
	}
	for(int i=0;i<V[root].size();i++)
	{
		tree &e=G[V[root][i]];
		if(e.to!=father[root][0]) //说明e.to是root的儿子
		{
			deep[e.to]=deep[root]+1;
			father[e.to][0]=root;  //建立父子关系
			dis[e.to][0]=e.w;      //更新到父亲的距离
			dfs(e.to); 
		} 
	} 
}

也可以把维护father过程单独拿出来写法类似RMQ

void beizeng()  
{
	for(int i=1;i<=mxdeep;i++)
	{
		for(int j=1;j<=n;j++)
		{
			father[j][i]=father[father[j][i-1]][i-1];
			dis[j][i]=dis[father[j][i-1]][i-1]+dis[j][i-1];
		}
	}
}

这里说一下为什么 father[root][i]=father[father[root][i-1]][i-1]。

因为father[root][i]代表 结点root的第 2^i 倍祖先的结点编号，所以father[root][0] 代表着root的父亲。这里的倍并不代表倍数而是距离，所以 2^(i-1)+2^(i-1)=2^i。

接着让深度大的先上移，但不能超过深度小的，之后如果没有移动到一起，再两个一起上移。

最后可能出现   root

                      u       v

那么答案就是u或者v的父亲。

LCA代码：

int lca(int u,int v)
{
	int ans=0;
	if(deep[u]>deep[v])  //让u结点不在v结点下方，方便操作 
	swap(u,v);
	for(int i=mxdeep;i>=0;i--)  //将 v结点上移，但不能高于u结点 
	{
		if(deep[u]<deep[v]&&deep[u]<=deep[father[v][i]])
		{
		//	printf("! %d %d %d %d\n",u,v,i,father[v][i]);
			ans+=dis[v][i];
			v=father[v][i];
		}
	}
	if(v==u)
	return ans;
	for(int i=mxdeep;i>=0;i--)
	{
		if(father[v][i]!=father[u][i])  //一起移动到共同父亲处 
		{
			ans+=dis[v][i]+dis[u][i];
			v=father[v][i];
			u=father[u][i];
		}
	}
	if(v!=u) ans+=dis[v][0]+dis[u][0];  // 如果没有移动到同一处，可能是   1->u   1->v ，那么向上移动一个位置就可以了 
	return ans;
}