HDU 4757 Tree 可持久化字典树 trie

http://acm.hdu.edu.cn/showproblem.php?pid=4757

给出一棵树，每个节点有权值，每次查询节点 (u,v) 以及 val，问 u 到 v 路径上的某个节点与 val 异或最大的值是多少。

和可持久化线段树差不多，看代码吧。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
#include<vector>
using namespace std;
#define LL long long
const int maxn=100010;
const int maxm=3000010;
int n,m;
int a[maxn]={},fa[maxn][20]={},dep[maxn]={},rt[maxn]={};
int sig[maxm][2]={},sz[maxm]={},cnt=0;
struct nod{int y,next;}e[maxn*2];int head[maxn],tot=0;
inline void init(int x,int y){e[++tot].y=y;e[tot].next=head[x];head[x]=tot;}
inline int getnew(){
    sig[++cnt][0]=0;sig[cnt][1]=0;
    sz[cnt]=0;return cnt;
}
inline void bui(int x,int pa,int val){
    rt[x]=getnew();pa=rt[pa];
    x=rt[x];
    for(int i=15;i>=0;--i){
        int ch=(val>>i)&1;
        if(!sig[x][ch]){
            int id=getnew();
            sig[x][ch]=id; sig[x][!ch]=sig[pa][!ch];
            sz[sig[x][ch]]=sz[sig[pa][ch]];
        }
        x=sig[x][ch];pa=sig[pa][ch];
        ++sz[x];
    }
}
void dfs(int x){
    bui(x,fa[x][0],a[x]);
    for(int i=1;fa[x][i-1];++i)fa[x][i]=fa[fa[x][i-1]][i-1];
    for(int i=head[x];i;i=e[i].next){
        if(e[i].y==fa[x][0])continue;
        fa[e[i].y][0]=x;dep[e[i].y]=dep[x]+1;dfs(e[i].y);
    }
}
inline int getlca(int x,int y){
    if(dep[x]<dep[y])swap(x,y);
    for(int i=19;i>=0;--i)if(dep[fa[x][i]]>=dep[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 Query(int x,int y,int val){
    int lc=getlca(x,y); int res=a[lc]^val;
    x=rt[x];y=rt[y];lc=rt[lc];
    int ret=0;
    for(int i=15;i>=0;--i){
        int ch=(val>>i)&1;
        if(sz[sig[x][!ch]]+sz[sig[y][!ch]]-2*sz[sig[lc][!ch]]>0){
            ret+=1<<i;
            ch=!ch;
        }
        x=sig[x][ch];y=sig[y][ch];lc=sig[lc][ch];
    }
    return max(res,ret);
}
int main(){
    while(~scanf("%d%d",&n,&m)){
        memset(head,0,sizeof(head));memset(rt,0,sizeof(rt));
        memset(sig,0,sizeof(sig));memset(fa,0,sizeof(fa));
        int x,y,z;tot=0;cnt=0;
        for(int i=1;i<=n;i++)scanf("%d",&a[i]);
        for(int i=1;i<n;i++){scanf("%d%d",&x,&y);init(x,y);init(y,x);}
        dep[1]=1;dfs(1);
        for(int i=1;i<=m;i++){
            scanf("%d%d%d",&x,&y,&z);
            printf("%d
",Query(x,y,z));
        }
    }
    return 0;
}