【LCA倍增】POJ1330-Nearest Common Ancestors

【知识点：离线算法&在线算法】

一个离线算法，在开始时就需要知道问题的所有输入数据，而且在解决一个问题后就要立即输出结果。

一个在线算法是指它可以以序列化的方式一个个的处理输入，也就是说在开始时并不需要已经知道所有的输入。

【倍增-在线算法】

水题……这道题范围卡得很紧…数组少了一位就会WA。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
const int MAXN=10000+10;
const int DEG=20;
vector<int> E[MAXN];
int n;
int anc[MAXN][DEG],dep[MAXN];

void addedge(int u,int v)
{
    E[u].push_back(v);
}

void dfs(int v,int p,int d)
{
    anc[v][0]=p;
    dep[v]=d;
    for (int i=0;i<E[v].size();i++)
        dfs(E[v][i],v,d+1);
}

void getanc()
{
    for (int i=1;i<DEG;i++)
        for (int j=1;j<=n;j++)
        {
            anc[j][i]=anc[anc[j][i-1]][i-1];
        }
}

int swim(int x,int H)
{
    for (int i=0;H>0;i++)
    {
        if (H&1) x=anc[x][i];
        H/=2;
    }
    return x;
}

int LCA(int u,int v)
{
    if (dep[u]<dep[v]) swap(u,v);
    u=swim(u,dep[u]-dep[v]);
    if (u==v) return u;
    for (int i=DEG-1;i>=0;i--)
    {
        if (anc[u][i]!=anc[v][i])
        {
            u=anc[u][i];
            v=anc[v][i];
        }
    }
    return anc[u][0];
}

void init()
{
    for (int i=0;i<MAXN;i++) vector<int>().swap(E[i]);
    scanf("%d",&n);
    int notrt[MAXN];
    memset(notrt,0,sizeof(notrt));
    for (int i=0;i<n-1;i++)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        addedge(u,v);
        notrt[v]=1;
    }
    int rt;
    for (int i=1;i<=n;i++) if (!notrt[i])    //这里注意下标是从1开始的
    {
        rt=i;
        break;
    }
    dfs(rt,0,0);
    getanc();
}

void query()
{
    int u,v;
    scanf("%d%d",&u,&v);
    cout<<LCA(u,v)<<endl;
}

int main()
{
    int T;
    scanf("%d",&T);
    for (int kase=0;kase<T;kase++)
    {
        init();
        query();
    }
    return 0;
}

 【待更新】