树链剖分,用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 ;
}

bzoj 4811: [Ynoi2017]由乃的OJ的更多相关文章

  1. BZOJ 4811 [Ynoi2017]由乃的OJ ——Link-Cut Tree

    直接维护按照顺序经过每一段,初始的1可以变成什么,初始为0可以变成什么. 然后答案就可以和起床困难综合征一样贪心处理了. 写起来并不好写. 发现交换左右子树之后答案会改变,GG 调了一天,最后还是T掉 ...

  2. bzoj 4811 由乃的OJ

    bzoj 4811 由乃的OJ 考虑树链剖分. 树剖后用一颗线段树维护一段连续区间,类似于一个函数,各位上进入 \(0/1\) ,输出的数字分别是什么.注意到最多只有 \(64\) 位,可以用一个 \ ...

  3. 【BZOJ4811】[Ynoi2017]由乃的OJ 树链剖分+线段树

    [BZOJ4811][Ynoi2017]由乃的OJ Description 由乃正在做她的OJ.现在她在处理OJ上的用户排名问题.OJ上注册了n个用户,编号为1-",一开始他们按照编号排名. ...

  4. BZOJ 3595: [Scoi2014]方伯伯的Oj SBT+可持久化Treap

    3595: [Scoi2014]方伯伯的Oj Time Limit: 6 Sec  Memory Limit: 256 MBSubmit: 102  Solved: 54[Submit][Status ...

  5. BZOJ4811 [Ynoi2017]由乃的OJ

    本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000 作者博客:http://www.cnblogs.com/ljh2000-jump/ ...

  6. 【bzoj4811】[Ynoi2017]由乃的OJ 树链剖分/LCT+贪心

    Description 给你一个有n个点的树,每个点的包括一个位运算opt和一个权值x,位运算有&,l,^三种,分别用1,2,3表示. 每次询问包含三个数x,y,z,初始选定一个数v.然后v依 ...

  7. BZOJ 4810 [Ynoi2017]由乃的玉米田 (莫队 + bitset)

    题目链接  BZOJ 4810 首先对询问离线, 莫队算法处理. 首先我们可以用bitset维护处当前区间中是否存在某个数. 对于询问1, 我们可以用 ((f >> q[i].x) &am ...

  8. bzoj4811 [Ynoi2017]由乃的OJ 树链剖分+位运算

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4811 因为位运算的结果有可合并性,所以可以树链剖分,线段树维护: 细节很多,特别要注意从左往 ...

  9. BZOJ 3595: [Scoi2014]方伯伯的Oj Splay + 动态裂点 + 卡常

    Description 方伯伯正在做他的Oj.现在他在处理Oj上的用户排名问题. Oj上注册了n个用户,编号为1-”,一开始他们按照编号排名.方伯伯会按照心情对这些用户做以下四种操作,修改用户的排名和 ...

随机推荐

  1. Java中的关键字

    1)48个关键字:abstract.assert.boolean.break.byte.case.catch.char.class.continue.default.do.double.else.en ...

  2. T-SQL 基础

    什么是T-SQL语言? T-SQL语言是SQL语言的扩展和增强,是与SQL server交流沟通的语言之一. T-SQL的组成部分,数据定义语言-DDL(对数据库对象的创建和管理),数据控制语言-DC ...

  3. Linux安装模式AppImage,Flatpak,Snap整理

    本文只谈Linux世界用户较多的前2大主要分支, RedHat Red Hat Enterprise Linux 简称RHEL rpm (RedHat, CentOS, Fedora, Oracle. ...

  4. Full Schema Stitching with Apollo Server

    转自: https://tomasalabes.me/blog/nodejs/graphql/apollo/2018/09/18/schema-stitiching-apollo.html Full ...

  5. Explicit

    Prefixing the explicit keyword to the constructor prevents the compiler from using that constructor ...

  6. Tanks!Tutorial 学习

    using UnityEngine; namespace Complete { public class CameraControl : MonoBehaviour { /// <summary ...

  7. laya的那些坑

    游戏运行在chrome里面 听不见声音 游戏运行在chrome里面  听不见声音:其它浏览器可以听见声音开发者模式提示如下: The AudioContext was not allowed to s ...

  8. optimization.splitChunks 中,chunks 的3个值:all、async、initial 的含义

    chunks 参数值 含义 all 把动态和非动态模块同时进行优化打包:所有模块都扔到 vendors.bundle.js 里面. initial 把非动态模块打包进 vendor,动态模块优化打包 ...

  9. React 中的 Component、PureComponent、无状态组件 之间的比较

    React 中的 Component.PureComponent.无状态组件之间的比较 table th:first-of-type { width: 150px; } 组件类型 说明 React.c ...

  10. jmeter ---断言持续时间(Duration to Assert )和断言文件大小

    jmeter ---断言持续时间(Duration to Assert )和断言文件大小 1.断言持续时间(Duration to Assert ) Duration to Assert -- 允许的 ...