POJ1330(倍增LCA解法)

传送门
题意：给你一颗树，让你找出两点的最近公共祖先。
题解：LCA模板题，这里我使用倍增的方法解决
我使用倍增来在线求 LCA ，时间和空间复杂度分别是 O((n+q)logn) 和 O(nlogn) 。q是查询次数。
算法基本思想：
　　　　①如果 a 和 b 深度不同，先把深度调浅，使他变得和浅的那个一样
　　　　②现在已经保证了 a 和 b 的深度一样，所以我们只要把两个一起一步一步往上移动，直到他们到达同一个节点，也就是他们的最近公共祖先了。
优化：而实际上，一步一步往上移动太慢，所以我们做一个预处理，我们用fa[i][j]表示第i个节点的向上的第2^j是哪个节点。然后在查找的时候，查找两次，第一次查找的时候，从大到小，找，找到a,b最差值最小的地方，保存。然后，再一次找，两个一起上升，此时depth差距最小，去寻找最近的公共祖先。其实过程类似于用一个二进制来表示一个十进制的数字，但是查找的次数变成了，那个十进制数的对应的二进制数的位数。

code

//#include <bits/stdc++.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <cstring>
using namespace std;
const int mx = 2e5 + 7;
typedef long long ll;
typedef unsigned  long long ull;
vector<ll> e[mx];
int n,depth[mx],fa[mx][20],h[mx];

void dfs(int pre,int rt){
    depth[rt] = depth[pre]+1;
    fa[rt][0] = pre;
    for(int i=1;i<20;i++)
        fa[rt][i] = fa[fa[rt][i-1]][i-1];
    for(int i=0;i<e[rt].size();i++)
        dfs(rt,e[rt][i]);
}
int LCA(int x,int y){
    if(depth[x]<depth[y]){
        swap(x,y);
    }
    if(x==y)
        return x;
    for(int i=19;i>=0;i--) {
        if (depth[x] - (1 << i) >= depth[y])
            x = fa[x][i];
    }
    if(x==y)
        return x;
    for(int i=19;i>=0;i--){
        if(fa[x][i]!=fa[y][i])
            x = fa[x][i],y=fa[y][i];
    }
    return fa[x][0];
}
int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
#ifdef ACM_LOCAL
    freopen("input.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
#endif
    int t;
    cin>>t;
    while (t--)
    {
        memset(h,0,sizeof(h));
        memset(depth,0,sizeof(depth));
        cin>>n;
        for(int i=1;i<=n-1;i++){
            int u,v;
            cin>>u>>v;
            e[u].push_back(v);
            h[v]++;
        }
        int rt=1;
        for(int i=1;i<=n;i++){
            if (h[i]==0) {
                rt = i;
                break;
            }
        }
        dfs(0,rt);
        int u,v;
        cin>>u>>v;
        cout<<LCA(u,v)<<endl;
        for(int i=1;i<=n;i++)
            e[i].clear();
    }
    return 0;
}