bzoj 4811: [Ynoi2017]由乃的OJ

树链剖分，用zkw线段树维护每条链两个方向上对每一位的变换情况，由于位数较少，可以用两个unsigned long long表示

#include<cstdio>
typedef unsigned long long u64;
const int N=;
char buf[N*],*ptr=buf-;
int _(){
    int x=,c=*++ptr;
    while(c<)c=*++ptr;
    while(c>)x=x*+c-,c=*++ptr;
    return x;
}
u64 _u64(){
    u64 x=;
    int c=*++ptr;
    while(c<)c=*++ptr;
    while(c>)x=x*+c-,c=*++ptr;
    return x;
}
int es[N*],enx[N*],e0[N],ep=,tp[N];
u64 v[N];
int fa[N],sz[N],top[N],dep[N],son[N],id[N],idp=,idr[N];
int n,m,k;
u64 vs[][];
void f1(int w,int pa){
    dep[w]=dep[fa[w]=pa]+;
    sz[w]=;
    for(int i=e0[w];i;i=enx[i]){
        int u=es[i];
        if(u!=pa){
            f1(u,w);
            sz[w]+=sz[u];
            if(sz[u]>sz[son[w]])son[w]=u;
        }
    }
}
void f2(int w,int tp){
    top[w]=tp;
    id[w]=++idp;
    idr[idp]=w;
    if(son[w])f2(son[w],tp);
    for(int i=e0[w];i;i=enx[i]){
        int u=es[i];
        if(u!=fa[w]&&u!=son[w])f2(u,u);
    }
}
int s1[],p1,s2[],p2,s3[],s4[],p3,p4;
void find(int l,int r,int tp){
    p1=p2=;
    for(l+=,r+=;r-l!=;l>>=,r>>=){
        if(~l&)s1[p1++]=l^;
        if(r&)s2[p2++]=r^;
    }
    if(!tp){
        for(int i=;i<p2;++i)s4[p4++]=s2[i];
        while(p1)s4[p4++]=s1[--p1];
    }else{
        for(int i=;i<p2;++i)s3[p3++]=s2[i];
        while(p1)s3[p3++]=s1[--p1];
    }
}
#define cal(a,b,c0,c1) (a)=((~(b)&(c0))|((b)&(c1)))
void init(int w,int tp,u64 v){
    u64*a=vs[w];
    if(tp==){
        a[]=a[]=0ll&v;
        a[]=a[]=-1ll&v;
    }else if(tp==){
        a[]=a[]=0ll|v;
        a[]=a[]=-1ll|v;
    }else if(tp==){
        a[]=a[]=0ll^v;
        a[]=a[]=-1ll^v;
    }
}
void up(u64*a,u64*l,u64*r){
    cal(a[],l[],r[],r[]);
    cal(a[],l[],r[],r[]);
    cal(a[],r[],l[],l[]);
    cal(a[],r[],l[],l[]);
}
u64 _v0,_v1;
void get(int x,int y){
    p3=p4=;
    int a=top[x],b=top[y];
    while(a!=b){
        if(dep[a]>dep[b]){
            find(id[a],id[x],);
            x=fa[a],a=top[x];
        }else{
            find(id[b],id[y],);
            y=fa[b],b=top[y];
        }
    }
    if(dep[x]>dep[y]){
        find(id[y],id[x],);
    }else{
        find(id[x],id[y],);
    }
    u64 v0=0ll,v1=-1ll;
    for(int i=;i<p4;++i){
        u64*a=vs[s4[i]];
        cal(v0,v0,a[],a[]);
        cal(v1,v1,a[],a[]);
    }
    while(p3){
        u64*a=vs[s3[--p3]];
        cal(v0,v0,a[],a[]);
        cal(v1,v1,a[],a[]);
    }
    _v0=v0,_v1=v1;
}
int main(){
    fread(buf,,sizeof(buf),stdin)[buf]=;
    n=_();m=_();k=_();
    for(int i=;i<=n;++i){
        tp[i]=_();
        v[i]=_u64();
    }
    for(int i=,a,b;i<n;++i){
        a=_();b=_();
        es[ep]=b;enx[ep]=e0[a];e0[a]=ep++;
        es[ep]=a;enx[ep]=e0[b];e0[b]=ep++;
    }
    f1(,);f2(,);
    for(int i=;i<=n;++i)init(id[i]+,tp[i],v[i]);
    for(int i=;i;--i)up(vs[i],vs[i<<],vs[i<<|]);
    for(int i=;i<m;++i){
        int o=_(),x=_(),y=_();
        u64 z=_u64();
        if(o==){
            get(x,y);
            u64 r=,r0,r1;
            for(int i=;i>=;--i)if((r|1llu<<i)<=z){
                cal(r0,r,_v0,_v1);
                cal(r1,(r|1llu<<i),_v0,_v1);
                if(r0<r1)r|=1llu<<i;
            }
            cal(r0,r,_v0,_v1);
            printf("%llu\n",r0);
        }else{
            init(id[x]+,tp[x]=y,v[x]=z);
            for(int w=id[x]+>>;w!=;w>>=)up(vs[w],vs[w<<],vs[w<<^]);
        }
    }
    return ;
}