洛谷P4332 [SHOI2014]三叉神经树（LCT）

传送门

FlashHu大佬太强啦%%%

首先，我们可以根据每一个点的权值为$1$的儿子的个数把每个点记为$0~3$，表示这一个点的点权

先考虑一下暴力的过程，假设从$0$变为$1$，先更改一个叶子结点，然后不断地往上更新，如果更改之后父亲的儿子中权值为$1$的儿子个数大于权值为$0$的儿子个数，那么就继续更新父亲，直到不能更新为止

不难发现，我们每一个更改的都是树上的一条链。而且，只有当点权为$1$时，才能因为修改使点权变为$2$，使输出改变，从而更新祖先。那么，我们需要的就是对于每一个节点，从它向上的连续的点权为$1$的一部分

如果是从$1$变为$0$同理，我们要维护从每一个节点向上的连续的点权为$2$的一部分

那么我们对于每一个节点，先access，然后记录splay中最深的点权不为$1$和$2$的点，然后把它给splay到根，再对右子树进行区间修改，对他自己单点修改

然后注意特判，如果整棵splay里都没有点权不为$1$或$2$的点，那么直接区间修改

然后怎么记录点权嘞？我们先把所有叶子结点的值*2，那么只有0和2两种值，然后每一个节点只要加上它所有叶子结点的值/2就是这个点的点权了

 //minamoto
 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 using namespace std;
 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
 char buf[<<],*p1=buf,*p2=buf;
 inline int read(){
     #define num ch-'0'
     char ch;bool flag=;int res;
     while(!isdigit(ch=getc()))
     (ch=='-')&&(flag=true);
     for(res=num;isdigit(ch=getc());res=res*+num);
     (flag)&&(res=-res);
     #undef num
     return res;
 }
 char sr[<<];int C=-;
 inline void Ot(){fwrite(sr,,C+,stdout),C=-;}
 inline void print(int x){
     sr[++C]=x|,sr[++C]='\n';
 }
 const int N=,M=;
 int fa[M],ch[N][],n1[N],n2[N],tag[N],v[M];
 int d[N],q[M],h=,t,ans;
 #define ls ch[x][0]
 #define rs ch[x][1]
 inline bool isrt(int x){return ch[fa[x]][]!=x&&ch[fa[x]][]!=x;}
 inline void pup(int x){
     //先右儿子再自己最后左儿子
     //因为&是短路运算符，有一个为假就会自动跳出
     if(!(n1[x]=n1[rs])&&!(n1[x]=x*(v[x]!=))) n1[x]=n1[ls];
     if(!(n2[x]=n2[rs])&&!(n2[x]=x*(v[x]!=))) n2[x]=n2[ls];
 }
 inline void pdown(int x,int t){
     //被区间修改的要么都是1要么都是2，直接反转信息
     v[x]^=,swap(n1[x],n2[x]),tag[x]+=t;
 }
 void all(int x){
     if(!isrt(x)) all(fa[x]);
     if(tag[x]) pdown(ls,tag[x]),pdown(rs,tag[x]),tag[x]=;
 }
 void rotate(int x){
     int y=fa[x],z=fa[y],d=ch[y][]==x;
     if(!isrt(y)) ch[z][ch[z][]==y]=x;
     fa[x]=z,fa[y]=x,fa[ch[x][d^]]=y,ch[y][d]=ch[x][d^],ch[x][d^]=y,pup(y);
 }
 void splay(int x){
     all(x);
     for(int y=fa[x],z=fa[y];!isrt(x);y=fa[x],z=fa[y]){
         if(!isrt(y))
         ((ch[y][]==x)^(ch[z][]==y))?rotate(x):rotate(y);
         rotate(x);
     }
     pup(x);
 }
 inline void access(int x){
     for(int y=;x;x=fa[y=x])
     splay(x),rs=y,pup(x);
 }
 int main(){
 //    freopen("testdata.in","r",stdin);
     int n=read();
     for(int i=;i<=n;++i)
     d[fa[read()]=fa[read()]=fa[read()]=i]=;
     for(int i=n+;i<=*n+;++i) v[q[++t]=i]=read()<<;
     while(h<=t){
         int x=q[h++];if(x<=n) pup(x);
         v[fa[x]]+=v[x]>>;
         if(!(--d[fa[x]])) q[++t]=fa[x];
     }
     ans=v[]>>;
     int q=read();
     while(q--){
         int x=read(),k=(v[x]^=)-;
         access(x=fa[x]),splay(x);
         if(~k?n1[x]:n2[x]){
             splay(x=(~k?n1[x]:n2[x]));
             pdown(rs,k),pup(rs);
             v[x]+=k,pup(x);
         }
         else pdown(x,k),pup(x),ans^=;
         print(ans);
     }
     Ot();
     return ;
 }