【YNOI2017】【BZOJ4811】【洛谷P3631】由乃的OJ / 睡觉困难综合征（树链剖分）（线段树）（贪心）

BZOJ传送门

洛谷传送门

---

解析：

稍微利用一点位运算技巧我们可以令一个状态记录为二元组$(x_0,x_1)$，其中$x_0$每一位表示这一位初始为$0$的数在最后这一位是$0$还是$1$，$x_1$同理。

定义状态加法为$(a_0,a_1)+(b_0,b_1)=(c_0,c_1)$，表示两个状态为$(a_0,a_1)$和$(b_0,b_1)$的区间拼接在一起后形成的状态，则稍加推导可以得到：

$$c_0=(a_0\&b_0)|(a_0\&b_1),c_1=(a_1\&b_1)|(a_1\&b_0)$$

于是就可以上树剖或者LCT来维护了，注意我们需要同时维护两个方向的答案，这个运算不满足交换律的。

提取一条链的$x_0,x_1$后就可以直接按位贪心了，和起床困难综合征这道题一个套路。

---

代码：

#include<bits/stdc++.h>
#define re register
#define gc get_char
#define cs const

typedef unsigned long long ull;

namespace IO{
	inline char get_char(){
		static cs int Rlen=1<<20|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++; 
	}
	
	template<typename T>
	inline T get(){
		re char c;
		while(!isdigit(c=gc()));re T num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return num;
	} 
}
using namespace IO;

using std::cout;
using std::cerr;

cs int N=1e5+5;

int n,m;

struct node{
	ull x0,x1;
	node():x0(0),x1(~0ull){}
	node(cs ull &x_0,cs ull &x_1):x0(x_0),x1(x_1){}
	void init(int typ,ull val){
		switch(typ){
			case 1:{
				x0=0;
				x1=val;
				break;
			}
			case 2:{
				x0=val;
				x1=~0ull;
				break;
			}
			case 3:{
				x0=val;
				x1=(~0ull)^val;
				break;
			}
		}
	}
	friend node operator+(cs node &a,cs node &b){
		return node((a.x0&b.x1)|((~a.x0)&b.x0),(a.x1&b.x1)|((~a.x1)&b.x0));
	}
};

int last[N],nxt[N<<1],to[N<<1],ecnt;
inline void addedge(int u,int v){
	nxt[++ecnt]=last[u],last[u]=ecnt,to[ecnt]=v;
	nxt[++ecnt]=last[v],last[v]=ecnt,to[ecnt]=u;
}

int fa[N],top[N],dep[N],siz[N],son[N];

void dfs1(int u){
	siz[u]=1;
	for(int re e=last[u],v=to[e];e;v=to[e=nxt[e]])if(v!=fa[u]){
		fa[v]=u;
		dep[v]=dep[u]+1;
		dfs1(v);
		siz[u]+=siz[v];
		if(siz[son[u]]<siz[v])son[u]=v;
	}
}

int pos[N],in[N],dfs_clock;

void dfs2(int u){
	pos[in[u]=++dfs_clock]=u;
	if(son[u]){
		top[son[u]]=top[u];
		dfs2(son[u]);
	}
	else return ;
	for(int re e=last[u],v=to[e];e;v=to[e=nxt[e]])if(v!=fa[u]&&v!=son[u]){
		top[v]=v;
		dfs2(v);
	}
}

int opt[N];
ull val[N];

node ln[N<<2],rn[N<<2];

inline void pushup(int k){
	ln[k]=ln[k<<1]+ln[k<<1|1];
	rn[k]=rn[k<<1|1]+rn[k<<1];
}

inline void build(int k,int l,int r){
	if(l==r){
		ln[k].init(opt[pos[l]],val[pos[l]]);
		rn[k]=ln[k];
		return ;
	}
	int mid=(l+r)>>1;
	build(k<<1,l,mid);
	build(k<<1|1,mid+1,r);
	pushup(k);
}

inline void modify(int k,int l,int r,cs int &p,cs int &typ,cs ull &vl){
	if(l==r){
		opt[pos[l]]=typ,val[pos[l]]=vl;
		ln[k].init(typ,vl);
		rn[k]=ln[k];
		return ;
	}
	int mid=(l+r)>>1;
	if(p<=mid)modify(k<<1,l,mid,p,typ,vl);
	else modify(k<<1|1,mid+1,r,p,typ,vl);
	pushup(k);
}

inline node queryl(int k,int l,int r,cs int &ql,cs int &qr){
	if(ql<=l&&r<=qr)return ln[k];
	int mid=(l+r)>>1;
	if(qr<=mid)return queryl(k<<1,l,mid,ql,qr);
	if(mid<ql)return queryl(k<<1|1,mid+1,r,ql,qr);
	return queryl(k<<1,l,mid,ql,qr)+queryl(k<<1|1,mid+1,r,ql,qr);
}

inline node queryr(int k,int l,int r,cs int &ql,cs int &qr){
	if(ql<=l&&r<=qr)return rn[k];
	int mid=(l+r)>>1;
	if(qr<=mid)return queryr(k<<1,l,mid,ql,qr);
	if(mid<ql)return queryr(k<<1|1,mid+1,r,ql,qr);
	return queryr(k<<1|1,mid+1,r,ql,qr)+queryr(k<<1,l,mid,ql,qr);
}

inline node path(int u,int v){
	node l,r;
	while(top[u]!=top[v]){
		dep[top[u]]>dep[top[v]]?
		(l=l+queryr(1,1,n,in[top[u]],in[u]),u=fa[top[u]]):
		(r=queryl(1,1,n,in[top[v]],in[v])+r,v=fa[top[v]]); 
	}
	(dep[u]>dep[v])?(l=l+queryr(1,1,n,in[v],in[u])):(l=l+queryl(1,1,n,in[u],in[v]));
	return l+r;
}

signed main(){
	n=get<int>(),m=get<int>();get<int>();
	for(int re i=1;i<=n;++i){
		opt[i]=get<int>();
		val[i]=get<ull>();
	}
	for(int re i=1;i<n;++i)addedge(get<int>(),get<int>());
	dfs1(1),top[1]=1,dfs2(1);
	build(1,1,n);
	while(m--)switch(get<int>()){
		case 1:{
			int u=get<int>(),v=get<int>();
			ull lim=get<ull>();
			node t=path(u,v);
			ull res=0;
			for(int re i=63;~i;--i){
				if(t.x0&(1ull<<i))res|=1ull<<i;
				else if(lim>=(1ull<<i)&&(t.x1&(1ull<<i)))res|=1ull<<i,lim-=1ull<<i;
			}
			cout<<res<<"\n";
			break;
		}
		case 2:{
			int u=get<int>(),typ=get<int>();
			ull x=get<ull>();
			modify(1,1,n,in[u],typ,x);
			break;
		} 
	} 
	return 0;
}