LCA最近公共祖先

倍增做法(其他的我不会

疯狂wa re 的原因：
1、数组开小了
2、题意从0出发 而初始化的j从1开始
3、ios解除输入和scanf混用(zz)
4、doubling 放在了dfs之前，应该在dfs之后

鬼知道我经历了什么

一些自己的理解

dp[i][j] 表示节点i 向上找 第 2^i 个点。这也就是倍增思维。比方说 节点a 和节点b相差了 13。13 表示成2进制就是 1 0 1 1 ， 那a到b 只需要先走最开头的1 (也就是先加上8) ， 则剩下的距离就是 a+8 和 b 的关系。完成了状态的转移
dp[i][j] = dp[ dp[i][j-1] ] [j-1] ， dp[i] [j-1] 放在上面的例子 就是 a+8 ，也就是a向上走8个单位的节点，
a到b 的距离为13 ，那么转移成 a -> a+8 -> b 。

例题

1、曼谷区域赛

#include<bits/stdc++.h>

#define rep(i,a,b) for(int i=a;i<b;i++)
typedef long long ll;
const int maxn = 400050;
const int K = 20;

using namespace std;
int T, n, m, kase = 1, p;
int pa[maxn][30], deep[maxn], sz[maxn];
vector<int> G[maxn];

int dfs(int u,int f)
{
    sz[u] = 1;
    for(int i=0,j=G[u].size(); i<j; i++)
    {
        int v = G[u][i];
        if(v == f)
            continue;
        deep[v] = deep[u] + 1;
        pa[v][0] = u;
        sz[u] += dfs(v,u);
    }
    return sz[u];
}
void doubling()
{
    rep(i,1,K) rep(j,1,n+1)
    {
        if(pa[j][i-1] == 0 || pa[pa[j][i-1]][i-1] == 0)
            pa[j][i] = 0;
        else
            pa[j][i] = pa[pa[j][i-1]][i-1];
    }
}
int lca(int x,int y)
{
    if(deep[x] < deep[y])
        swap(x,y);
    for(int i = K-1; ~i; i--)
        if(deep[pa[x][i]] > deep[y])
            x = pa[x][i];
    if(deep[x] != deep[y])
        x = pa[x][0];
    for(int i = K-1; ~i; i--)
    {
        if(pa[x][i] != pa[y][i])
            x = pa[x][i], y = pa[y][i];
    }
    if(x != y)
        x = pa[x][0], y = pa[y][0];
    return x;
}

int query(int x, int p)
{
    // p 找到 x 下面最接近的地方
    for(int i=K-1; ~i; i--)
    {
        if(deep[pa[p][i]] > deep[x])
            p = pa[p][i];
    }
    // cout<<"find d == "<<p<<endl;
    return p;
}

void solve(int x,int p)
// 关键点  首都
{
    int LCA = lca(x,p);
    //cout<<"LCA == "<<LCA<<endl;

    if(LCA == x)
    {
        int d = query(x,p);
        printf("%d\n",n-sz[d]);
    }
    else
    {
        printf("%d\n", sz[x]);
    }
}

int main()
{

    cin >> T;
    int cas = 1;
    while(T--)
    {

        deep[1] = 1;
        cin >> n >> m >> p;
        rep(i,0,n+1)
        G[i].clear();
        rep(i,1,n)
        {
            int a,b;
            scanf("%d%d",&a,&b);
            G[a].push_back(b);
            G[b].push_back(a);
        }
        dfs(1,0);
        doubling();
        printf("Case #%d:\n",cas++);
//        rep(i,1,n+1)
//        {
//            cout<<"sz dp[i][0] deep[i] "<<sz[i]<<' '<<pa[i][0]<<' '<<deep[i]<<endl;
//        }
        while(m--)
        {
            int opt,a;
            scanf("%d%d",&opt,&a);
            if(opt == 0)
            {
                p = a;

            }
            else //if(opt == 2)
            {
                if(a == p)
                {
                    printf("%d\n",n);
                    continue;
                }
                solve(a,p);
            }

//            else
//            {
//                cout<<lca(opt,a)<<endl;
//            }
        }

    }
    return 0;
}


/*

10
8 100 1
1 2
2 3
1 4
4 5
4 6
6 7
6 8
*/