洛谷3613睡觉困难综合征（LCT维护链信息（前后缀）+贪心）

这个题目还是很好啊QWQ很有纪念意义

首先，如果在序列上且是单次询问的话，就是一个非常裸的贪心了QWQ这也是NOI当时原题的问题和数据范围

我们考虑上树的话，应该怎么做？

我的想法是，对于每一位建一个LCT来做，然后对一个点x维护，当前是0到链顶的结果，当前是1到顶的结果，当前是0到底的结果，当前是1到底的结果（之所以维护后两个是因为(reverse)操作的存在）

这样的话，对于每次询问，我就可以直接(split)，然后贪心的询问。

不过很可惜，这个算法的复杂度是(O(qklogn))的
并不能跑过（虽然没试过硬卡常行不行）

那么应该怎么做呢QWQ 这时候，我选择了看题解。

首先，经过仔细思考，我们会发现，QWQ对于题目来说，其实并不需要按位来做，可以直接维护(f_0,f_1,g_0,g_1)四个值，分别表示所有位0到链顶的结果，所有二进制位全是1的数到顶的结果，所有位全是0到底的结果，所有位全是1到底的结果

那么这个东西应该怎么转移呢？

这里合并的大概思想是 全0走到中间节点之后，是0的那几位，相当于后面的全0走的，是1的那几位，相当于是后面全1走的 ，全1的也是同理

而之所以这么写，是因为运用了与的美妙的性质，只有两边都是1，最终结果才是1。

所以，假设对于当前是全零0来说，考虑前一半，是计算0的那几位，要是采用了&与这个运算，如果最后走完，(f_0)的某一位是0，那么最一开始无论是什么，最终&完都是0，也就不会存在一开始是1的那些位影响答案。

要是这一位是1，我们要满足与完 是1的话，同时去除一开始是1的位的贡献，我们就需要先取反，然后再&

其他的其实也是同理

Node merge(Node a,Node b) //我们默认a在前，b在后 
{
    Node c;
    c.f0=(~a.f0 & b.f0) + (a.f0 & b.f1); //这里合并的大概思想是 全0走到中间节点之后，是0的那几位，相当于后面的全0走的，是1的那几位，相当于是后面全1走的 
    c.f1=(~a.f1 & b.f0) + (a.f1 & b.f1); //而之所以这么写，是因为运用了与的美妙的性质，只有两边都是1，最终结果才是1。所以，假设对于当前是0来说，只有后面全零弄出来是1，我们当前的的答案才是1，那么为了出来1，我们就需要先取反，然后再& 
    return c;
}

void update(unsigned long long x)
{
    zheng[x]=fan[x]=val[x];
    if (ch[x][0])
    {
        zheng[x]=merge(zheng[x],zheng[ch[x][0]]);
        fan[x]=merge(fan[ch[x][0]],fan[x]); 
    }
    if (ch[x][1])
    {
        zheng[x]=merge(zheng[ch[x][1]],zheng[x]);
        fan[x]=merge(fan[x],fan[ch[x][1]]); 
    }
}

其他的LCT操作的话，也就和别的没什么区别了。

这里再讲一下底下那个贪心的过程

unsigned long long ans=ling;
     	split(x,y);
     	//cout<<1<<endl;
     	for (long long i=k-1;i>=0;i--)
     	{
     		if (fan[y].f0&power[i]) 
            {
              ans=ans+power[i];
              continue;
            }
            if (fan[y].f1 & power[i])
            {
                if (power[i]>z) continue;
                ans=ans+power[i];
                z-=power[i];
            }
            //cout<<i<<endl;
        }
        cout<<ans<<"
";

从高位向低位考虑，如果当前位填0，最终结果是1的话，那么就填0。
如果当前位填1，最终结果才能是1。我们就需要比较一下剩余的值是否比这个位填1的数大，大的话才能填

以此类推

QWQ

记得开(unsigned long long)

上代码

// luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<set>
#define mk makr_pair
#define ll long long
#define lc ch[x][0]
#define rc ch[x][1]

using namespace std;

inline unsigned long long read()
{
  unsigned long long x=0,f=1;char ch=getchar();
  while (!isdigit(ch)) {if (ch=='-') f=-1;ch=getchar();}
  while (isdigit(ch)) {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
  return x*f;
}

const int maxn = 1e6+1e2;

struct Node{
	unsigned long long f0,f1;
};
unsigned long long ch[maxn][3];
unsigned long long rev[maxn],fa[maxn];
Node zheng[maxn],fan[maxn];
unsigned long long n,m;
Node val[maxn];
unsigned long long ling =0;

unsigned long long son(unsigned long long x)
{
	if (ch[fa[x]][0]==x) return 0;
	else return 1; 
}

bool notroot(unsigned long long x)
{
	return ch[fa[x]][0]==x || ch[fa[x]][1]==x;
}

void reverse(unsigned long long x)
{
	swap(ch[x][0],ch[x][1]);
	swap(zheng[x].f0,fan[x].f0);
	swap(zheng[x].f1,fan[x].f1);
	rev[x]^=1;
}

Node merge(Node a,Node b) //我们默认a在前，b在后 
{
	Node c;
	c.f0=(~a.f0 & b.f0) + (a.f0 & b.f1); //这里合并的大概思想是 全0走到中间节点之后，是0的那几位，相当于后面的全0走的，是1的那几位，相当于是后面全1走的 
	c.f1=(~a.f1 & b.f0) + (a.f1 & b.f1); //而之所以这么写，是因为运用了与的美妙的性质，只有两边都是1，最终结果才是1。所以，假设对于当前是0来说，只有后面全零弄出来是1，我们当前的的答案才是1，那么为了出来1，我们就需要先取反，然后再& 
	return c;
}

void update(unsigned long long x)
{
	zheng[x]=fan[x]=val[x];
	if (ch[x][0])
	{
		zheng[x]=merge(zheng[x],zheng[ch[x][0]]);
		fan[x]=merge(fan[ch[x][0]],fan[x]); 
	}
	if (ch[x][1])
	{
		zheng[x]=merge(zheng[ch[x][1]],zheng[x]);
		fan[x]=merge(fan[x],fan[ch[x][1]]); 
	}
}

void pushdown(unsigned long long x)
{
	if (rev[x])
	{
		if (ch[x][1]) reverse(ch[x][1]);
		if (ch[x][0]) reverse(ch[x][0]);
		rev[x]=0;
	}
}

void rotate(unsigned long long x)
{
	unsigned long long y=fa[x],z=fa[y];
	unsigned long long b=son(x),c=son(y);
	if (notroot(y)) ch[z][c]=x;
	fa[x]=z;
	ch[y][b]=ch[x][!b];
	fa[ch[x][!b]]=y;
	ch[x][!b]=y;
	fa[y]=x;
	update(y);
	update(x);
}

unsigned long long st[maxn];

void splay(unsigned long long x)
{
	unsigned long long y=x,cnt=0;
	st[++cnt]=y;
	while (notroot(y)) y=fa[y],st[++cnt]=y;
	while (cnt) pushdown(st[cnt--]);
	while (notroot(x))
	{
		unsigned long long y=fa[x],z=fa[y];
		unsigned long long b=son(x),c=son(y);
		if (notroot(y))
		{
			if (b==c) rotate(y);
			else rotate(x); 
		}
		rotate(x);
	 } 
	 update(x);
}

void access(unsigned long long x)
{
	for (unsigned long long y=0;x;y=x,x=fa[x])
	{
		splay(x);
		ch[x][1]=y;
		update(x);
	}
}

void makeroot(unsigned long long x)
{
	access(x);
	splay(x);
	reverse(x);
}

unsigned long long findroot(unsigned long long x)
{
	access(x);
	splay(x);
    while (ch[x][0])
	{
		pushdown(x);
		x=ch[x][0];
	} 
	return x; 
}

void split(unsigned long long x,unsigned long long y)
{
	makeroot(x);
	access(y);
	splay(y);
}

void link(unsigned long long x,unsigned long long y)
{
	makeroot(x);
	if(findroot(y)!=x) fa[x]=y;
}

unsigned long long power[maxn];

unsigned long long ymh = 2;
unsigned long long k;

signed main()
{
  n=read(),m=read(),k=read();
  power[0]=1;
  for (unsigned long long i=1;i<64;i++) power[i]=power[i-1]*ymh;
  
  for (unsigned long long i=1;i<=n;i++)
  {
  	unsigned long long y=read(),x=read();
  	if (y==1)
  	{
  		val[i]=(Node){ling,x};
	}
	if (y==2)
	{
		val[i]=(Node){x,~ling};
	}
	if (y==3)
	{
		val[i]=(Node){x,~x};
	} 
  }
  for (unsigned long long i=1;i<n;i++)
  {
  	unsigned long long x=read(),y=read();
  	link(x,y);
  } 
  
  //cout<<1<<endl;
  for (unsigned long long i=1;i<=m;i++)
  {
  	 unsigned long long opt=read(),x=read(),y=read(),z=read();
  	 if (opt==2)
  	 {
	   splay(x);
  	   if (y==1) val[x]=(Node){ling,z};
	   if (y==2) val[x]=(Node){z,~ling};
	   if (y==3) val[x]=(Node){z,~z}; 
     }
     if (opt==1)
     {
     	unsigned long long ans=ling;
     	split(x,y);
     	//cout<<1<<endl;
     	for (long long i=k-1;i>=0;i--)
     	{
     		if (fan[y].f0&power[i]) 
			{
			  ans=ans+power[i];
			  continue;
			}
			if (fan[y].f1 & power[i])
			{
				if (power[i]>z) continue;
				ans=ans+power[i];
				z-=power[i];
			}
			//cout<<i<<endl;
		}
		cout<<ans<<"
";
	 }
  }
  return 0;
}