hdu5390 tree

　　先求出dfs序，然后建立线段树，线段树每个节点套一个set。

　　修改操作只需要改被子树区间完全覆盖的线段树节点，将其节点中set的原来的值删除，加入新值。

　　询问操作查询单点到根的所有节点上的set中与查询值异或起来最大的那个。

　　查询set中的数与x异或的最大值，可以从高位到低位枚举二进制位，根据x的二进制位，查询一些set中数值的存在情况，例如加入x的二进制第y位为1，那么如果set中存在第y位为0的数字，明显可以使得答案更大。下面代码中这段的代码类似于二分。

　　时间复杂度O(nlogn^3)

　　

　　代码

#include<cstdio>
#include<set>
#include<algorithm>
#include<cstring>
#define N 500010
using namespace std;
int n,m,i,a,b,typ;
int dp,p[N],pre[N],tt[N];
int tot,L[N],R[N],l[N],r[N],v[N],stack[N],deep,flag;
multiset<int> Set[N];
void link(int x,int y)
{
    dp++;pre[dp]=p[x];p[x]=dp;tt[dp]=y;
}
void dfs()
{
    int i;
    while (deep)
    {
        if (L[stack[deep]]==0)
        {
            tot++;
            L[stack[deep]]=tot;
        }
        if (p[stack[deep]])
        {
            i=p[stack[deep]];
            deep++;stack[deep]=tt[i];
            p[stack[deep-1]]=pre[i];
        }
        else
        {
            R[stack[deep]]=tot;
            deep--;
        }
    }
}
void build(int x,int a,int b)
{
    int m;
    l[x]=a;r[x]=b;
    Set[x].clear();
    Set[x].insert(-1);
    Set[x].insert(1<<30);
    if (b-a>1)
    {
        m=(a+b)>>1;
        build(2*x,a,m);
        build(2*x+1,m,b);
    }
}
void change(int x,int a,int b,int c,int d)
{
    if ((a<=l[x])&&(r[x]<=b))
    {
        Set[x].insert(d);
        if (flag)
        Set[x].erase(c);
        return;
    }
    int m=(l[x]+r[x])>>1;
    if (a<m) change(2*x,a,b,c,d);
    if (m<b) change(2*x+1,a,b,c,d);
}
int Q(int x,int y)
{
    int l,r,m,p,q;
    if (Set[x].size()==2) return 0;
    l=0;r=(1<<30)-1;
    while (l!=r)
    {
        m=(l+r)>>1;
        multiset<int>::iterator it=Set[x].upper_bound(m);
        p=*it;
        q=*(--it);
        if (p>r)
            r=m;
        else
        if (q<l)
            l=m+1;
        else
        {
            if ((l^y)>((m+1)^y))
            r=m;
            else
            l=m+1;
        }
    }
    return l^y;
}
int query(int x,int a,int b,int c)
{
    int ans=Q(x,c);
    if ((a<=l[x])&&(r[x]<=b))
    return ans;
    int m=(l[x]+r[x])>>1;
    if (a<m) ans=max(ans,query(2*x,a,b,c));
    if (m<b) ans=max(ans,query(2*x+1,a,b,c));
    return ans;
}
int main()
{
    int test;
    scanf("%d",&test);
    while (test--)
    {
    scanf("%d%d",&n,&m);
    dp=0;
    for (i=1;i<=n;i++)
    {
        L[i]=0;
        p[i]=0;
    }
    for (i=2;i<=n;i++)
    {
        scanf("%d",&a);
        link(a,i);
    }
    tot=0;
    deep=1;stack[1]=1;
    dfs();

    build(1,0,n);
    flag=0;
    for (i=1;i<=n;i++)
    {
        scanf("%d",&v[i]);
        change(1,L[i]-1,R[i],0,v[i]);
    }
    flag=1;
    for (i=1;i<=m;i++)
    {
        scanf("%d",&typ);
        if (typ==0)
        {
            scanf("%d%d",&a,&b);
            change(1,L[a]-1,R[a],v[a],b);
            v[a]=b;
        }
        else
        {
            scanf("%d",&a);
            printf("%d
",query(1,L[a]-1,L[a],v[a]));
        }
    }
    }
}