Vijos.lxhgww的奇思妙想(k级祖先 长链剖分)

题目链接

https://blog.bill.moe/long-chain-subdivision-notes/
http://www.cnblogs.com/zzqsblog/p/6700133.html

长链剖分模板。
6,7,8,20个点RE，没什么办法了。。（迷）

#include <cstdio>
#include <cctype>
#include <vector>
#include <algorithm>
#define gc() getchar()
#define MAXIN 300000
//#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
#define BIT 18
typedef long long LL;
const int N=3e5+5;

int Enum,H[N],nxt[N<<1],to[N<<1],fa[N][19],dep[N],son[N],mxd[N],top[N],hb[N];
std::vector<int> up[N],down[N];
char IN[MAXIN],*SS=IN,*TT=IN;

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}
inline void AE(int u,int v)
{
    to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
    to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum;
}
void DFS1(int x)
{
    int mx=dep[x];
    for(int i=1; 1<<i<=mx; ++i) fa[x][i]=fa[fa[x][i-1]][i-1];
    for(int i=H[x],v; i; i=nxt[i])
        if((v=to[i])!=fa[x][0])
            fa[v][0]=x, dep[v]=dep[x]+1, DFS1(v), mxd[v]>mx&&(mx=mxd[v],son[x]=v);
    mxd[x]=mx;
}
void DFS2(int x,int tp)
{
    top[x]=tp;
    if(son[x])
    {
        DFS2(son[x],tp);
        for(int i=H[x],v; i; i=nxt[i])
            if((v=to[i])!=fa[x][0]&&v!=son[x]) DFS2(v,v);
    }
}
inline int Query(int k,int x)
{
    if(!k) return x;
    if(k>dep[x]) return 0;
    x=fa[x][hb[k]], k^=1<<hb[k];
    return dep[x]-dep[top[x]]>=k?down[top[x]][dep[x]-dep[top[x]]-k]:up[top[x]][k-dep[x]+dep[top[x]]];
}

int main()
{
    int n=read();
    for(int i=1; i<n; ++i) AE(read(),read());
    DFS1(1), DFS2(1,1);
    for(int i=1; i<=n; ++i)
        if(top[i]==i)
        {
            int lim=mxd[i]-dep[i];
            for(int x=i,j=0; j<=lim; ++j) up[i].push_back(x),x=fa[x][0];
            for(int x=i,j=0; j<=lim; ++j) down[i].push_back(x),x=son[x];
        }
    for(int mx=1,i=2; i<=n; ++i) i>>mx&1&&(++mx), hb[i]=mx-1;
//  for(int i=2; i<=n; ++i) hb[i]=hb[i-1]+(i>>(hb[i-1]+1)&1);
    for(int ans=0,Q=read(); Q--; printf("%d\n",ans=Query(ans^read(),ans^read())));

    return 0;
}

转载于:https://www.cnblogs.com/SovietPower/p/10018890.html