原帖在这:http://blog.csdn.net/crazy_ac/article/details/7796497
LCA可以离线tarjan去求,也可以在线倍增去求。
记二维数组p[u][i]表示u的2^i个祖先是谁,i==0就是父节点。那么可以通过 p[u][i]=p[p[u][i-1]][i-1]来递推出所有的情况,然后再求lca(x,y)的时候,根据x,y的深度,先将深度大的点通过p[][放缩到与x相同的深度,做法是用二进制表示出深度差,然后从低到高扫描,若第i位为1,那么b=p[b][i],否则不变即第i位为1时,b跳转到其2^i个祖先处,看起来有点抽象,不过画个图或者单步监视一下就明白了。放缩到统一深度后,直接枚举一下答案就出来了...
这里放一个poj1330的模板,题意就是给一个数,然后给一个查询求LCA。
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;
typedef long long ll;
const int POW = 16;
const int maxn=10000;
vector<int> edge[maxn];
int p[maxn][16];
int d[maxn];
int n,m,k;
void dfs(int u,int fa){
d[u]=d[fa]+1;
p[u][0]=fa;
for(int i=1;i<POW;i++) p[u][i]=p[p[u][i-1]][i-1];
int sz=edge[u].size();
for(int i=0;i<sz;i++){
int v=edge[u][i];
if(v==fa) continue;
dfs(v,u);
}
}
int lca( int a, int b ){
if( d[a] > d[b] ) a ^= b, b ^= a, a ^= b;
if( d[a] < d[b] ){
int del = d[b] - d[a];
for( int i = 0; i < POW; i++ ) if(del&(1<<i)) b=p[b][i];
}
if( a != b ){
for( int i = POW-1; i >= 0; i-- )
if( p[a][i] != p[b][i] )
a = p[a][i] , b = p[b][i];
a = p[a][0], b = p[b][0];
}
return a;
}
int ideg[maxn];
int main()
{
// freopen("in.txt","r",stdin);
int tt;
scanf("%d",&tt);
while (tt--)
{
memset(d,0,sizeof d);
memset(p,0,sizeof p);
memset(ideg,0,sizeof ideg);
int x,y;
scanf("%d",&n);
for (int i=1; i<=n; i++)
edge[i].clear();
for (int i=1; i<n; i++)
{
scanf("%d%d",&x,&y);
ideg[y]++;
edge[x].push_back(y);
edge[y].push_back(x);
}
int root;
for (int i=1; i<=n; i++)
if (ideg[i]==0) root=i;
dfs(root,-1);
scanf("%d%d",&x,&y);
cout<<lca(x,y)<<endl;
}
return 0;
}