2017 ACM-ICPC 亚洲区（西安赛区）网络赛 G. Xor(LCA+K级父亲暴力卡过)

题目链接：2017 ACM-ICPC 亚洲区（西安赛区）网络赛   G. Xor

题意：

给你一棵有n个节点的树，然后有Q个询问，每次询问给出两个点x,y,k。将x->y路径上经过的点放进一个数组a里，将询问[0],a[k],a[2*k],a[m*k]的值异或起来的值为多少。

题解：

预处理lca和k=1的答案，k>1的用倍增暴力往上跳，然后就完了，不过有点慢，需要加个读入挂才能过。

#include<bits/stdc++.h>
#define F(i,a,b) for(int i=(a);i<=(b);++i)
using namespace std;

const int N=5e5+5,DEG=17;
int ed,g[N],nxt[2*N],v[2*N],fa[N][DEG],dep[N],val[N];
int a[N],n,q,vis[N],TIM;

inline void adg(int x,int y){v[++ed]=y,nxt[ed]=g[x],g[x]=ed;}

void dfs(int u,int pre){
    dep[u]=dep[pre]+1,fa[u][0]=pre,val[u]=val[pre]^a[u];
    F(i,1,DEG-1)fa[u][i]=fa[fa[u][i-1]][i-1];
    for(int i=g[u];i;i=nxt[i])if(v[i]!=pre)dfs(v[i],u);
}

int LCA(int a,int b){
    if(dep[a]>dep[b])a^=b,b^=a,a^=b;
    if(dep[a]<dep[b])F(i,0,DEG-1)if((dep[b]-dep[a])&(1<<i))b=fa[b][i];
    if(a!=b)for(int i=DEG-1;i<0?a=fa[a][0]:0,i>=0;i--)
    if(fa[a][i]!=fa[b][i])a=fa[a][i],b=fa[b][i];
    return a;
}


int kth_father(int x,int k)
{
    while(k)x=fa[x][__builtin_ctz(k&-k)],k^=k&-k;
    return x;
}

void solve(int x,int y,int k)
{
    TIM++;
    if(x==y){printf("%d
",a[x]);return;}
    int lca=LCA(x,y);
    if(k==1)
    {
        if(lca==x)printf("%d
",val[y]^val[lca]^a[lca]);
        else if(lca==y)printf("%d
",val[x]^val[lca]^a[lca]);
        else printf("%d
",val[y]^val[x]^a[lca]);
        return;
    }
    int now=kth_father(x,k),ans=a[x],pre=now;
    vis[x]=TIM;
    while(dep[now]>=dep[lca])
    {
        if(vis[now]!=TIM)ans^=a[now],vis[now]=TIM;
        pre=now,now=kth_father(now,k);
    }
    int Long=dep[x]+dep[y]-2*dep[lca]+1;
    if((Long-1)<k){printf("%d
",ans);return;}
    int last=(Long-1)/k*k+1;
    int dis=Long-last;
    now=kth_father(y,dis);
    if(vis[now]!=TIM&&dep[now]>=dep[lca])ans^=a[now],vis[now]=TIM;
    now=kth_father(now,k);
    while(dep[now]>=dep[lca])
    {
        if(vis[now]!=TIM)ans^=a[now],vis[now]=TIM;
        now=kth_father(now,k);
    }
    printf("%d
",ans);
}

int main(){
    while(read(n),read(q))
    {
        int x,y,k;
        F(i,1,n)g[i]=0;ed=0;
        F(i,2,n)
        {
            read(x),read(y);
            adg(x,y),adg(y,x);
        }
        F(i,1,n)read(a[i]);
        dfs(1,0);
        F(i,1,q)
        {
            read(x),read(y),read(k);
            solve(x,y,k);
        }
    }
    return 0;
}