bzoj 4811: [Ynoi2017]由乃的OJ

#include<cstdio>
typedef unsigned long long u64;
const int N=100007;
char buf[N*100],*ptr=buf-1;
int _(){
int x=0,c=*++ptr;
while(c<48)c=*++ptr;
while(c>47)x=x*10+c-48,c=*++ptr;
return x;
}
u64 _u64(){
u64 x=0;
int c=*++ptr;
while(c<48)c=*++ptr;
while(c>47)x=x*10+c-48,c=*++ptr;
return x;
}
int es[N*2],enx[N*2],e0
,ep=2,tp
;
u64 v
;
int fa
,sz
,top
,dep
,son
,id
,idp=0,idr
;
int n,m,k;
u64 vs[262777][4];
void f1(int w,int pa){
dep[w]=dep[fa[w]=pa]+1;
sz[w]=1;
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[107],p1,s2[107],p2,s3[10007],s4[10007],p3,p4;
void find(int l,int r,int tp){
p1=p2=0;
for(l+=131071,r+=131073;r-l!=1;l>>=1,r>>=1){
if(~l&1)s1[p1++]=l^1;
if(r&1)s2[p2++]=r^1;
}
if(!tp){
for(int i=0;i<p2;++i)s4[p4++]=s2[i];
while(p1)s4[p4++]=s1[--p1];
}else{
for(int i=0;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==1){
a[0]=a[2]=0ll&v;
a[1]=a[3]=-1ll&v;
}else if(tp==2){
a[0]=a[2]=0ll|v;
a[1]=a[3]=-1ll|v;
}else if(tp==3){
a[0]=a[2]=0ll^v;
a[1]=a[3]=-1ll^v;
}
}
void up(u64*a,u64*l,u64*r){
cal(a[0],l[0],r[0],r[1]);
cal(a[1],l[1],r[0],r[1]);
cal(a[2],r[2],l[2],l[3]);
cal(a[3],r[3],l[2],l[3]);
}
u64 _v0,_v1;
void get(int x,int y){
p3=p4=0;
int a=top[x],b=top[y];
while(a!=b){
if(dep[a]>dep[b]){
find(id[a],id[x],0);
x=fa[a],a=top[x];
}else{
find(id[b],id[y],1);
y=fa[b],b=top[y];
}
}
if(dep[x]>dep[y]){
find(id[y],id[x],0);
}else{
find(id[x],id[y],1);
}
u64 v0=0ll,v1=-1ll;
for(int i=0;i<p4;++i){
u64*a=vs[s4[i]];
cal(v0,v0,a[2],a[3]);
cal(v1,v1,a[2],a[3]);
}
while(p3){
u64*a=vs[s3[--p3]];
cal(v0,v0,a[0],a[1]);
cal(v1,v1,a[0],a[1]);
}
_v0=v0,_v1=v1;
}
int main(){
fread(buf,1,sizeof(buf),stdin)[buf]=0;
n=_();m=_();k=_();
for(int i=1;i<=n;++i){
tp[i]=_();
v[i]=_u64();
}
for(int i=1,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(1,0);f2(1,1);
for(int i=1;i<=n;++i)init(id[i]+131072,tp[i],v[i]);
for(int i=131071;i;--i)up(vs[i],vs[i<<1],vs[i<<1|1]);
for(int i=0;i<m;++i){
int o=_(),x=_(),y=_();
u64 z=_u64();
if(o==1){
get(x,y);
u64 r=0,r0,r1;
for(int i=63;i>=0;--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]+131072,tp[x]=y,v[x]=z);
for(int w=id[x]+131072>>1;w!=1;w>>=1)up(vs[w],vs[w<<1],vs[w<<1^1]);
}
}
return 0;
}