bzoj3720 Gty的妹子树

Description

我曾在弦歌之中听过你，
檀板声碎，半出折子戏。
舞榭歌台被风吹去，
岁月深处尚有余音一缕……

Gty神(xian)犇(chong)从来不缺妹子……
他来到了一棵妹子树下，发现每个妹子有一个美丽度……
由于Gty很哲♂学，他只对美丽度大于某个值的妹子感兴趣。
他想知道某个子树中美丽度大于k的妹子个数。
某个妹子的美丽度可能发生变化……
树上可能会出现一只新的妹子……

维护一棵初始有n个节点的有根树(根节点为1)，树上节点编号为1-n，每个点有一个权值wi。
支持以下操作：
0 u x          询问以u为根的子树中，严格大于x的值的个数。(u^=lastans,x^=lastans)
1 u x          把u节点的权值改成x。(u^=lastans,x^=lastans)
2 u x          添加一个编号为"当前树中节点数+1"的节点，其父节点为u，其权值为x。(u^=lastans,x^=lastans)
最开始时lastans=0。

Input

输入第一行包括一个正整数n(1<=n<=30000)，代表树上的初始节点数。
接下来n-1行，每行2个整数u,v，为树上的一条无向边。
任何时刻，树上的任何权值大于等于0，且两两不同。
接下来1行，包括n个整数wi，表示初始时每个节点的权值。
接下来1行，包括1个整数m(1<=m<=30000)，表示操作总数。
接下来m行，每行包括三个整数 op,u,v：
op,u,v的含义见题目描述。
保证题目涉及的所有数在int内。

Output

对每个op=0，输出一行，包括一个整数，意义见题目描述。

Sample Input

2
1 2
10 20
1
0 1 5

Sample Output

2

正解：块状树。

比较$simple$的块状树吧，不过挺容易写错的。。

为了防止被卡，我直接设块的大小为$245$，然后直接按照$size$分块就行了。

加入的话就看它父亲所在的块的大小是否大于$size$，如果大于就分一个新块。

然后我们对于每个块，弄一个动态数组，把块内的权值排好序插进去。

修改的话直接把整个块的动态数组清空，再重新排序加入就行了。

查询的话我们还要把相邻的块连边，按照当前子树是否包括一个整块分类讨论一下就行了。。

//It is made by wfj_2048~
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
#define inf (1<<30)
#define N (500010)
#define il inline
#define RG register
#define ll long long
#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)

using namespace std;

struct edge{ int nt,to; }G[N],g[N];

int Head[N],head[N],sz[N],st[N],top[N],fa[N],bl[N],a[N],w[N],T,n,m,num,Num,ans,cnt,block;

vector <int> val[N];

il int gi(){
    RG int x=0,q=1; RG char ch=getchar();
    while ((ch<'0' || ch>'9') && ch!='-') ch=getchar();
    if (ch=='-') q=-1,ch=getchar();
    while (ch>='0' && ch<='9') x=x*10+ch-48,ch=getchar();
    return q*x;
}

il void insert(RG int from,RG int to){
    g[++num]=(edge){head[from],to},head[from]=num; return;
}

il void link(RG int from,RG int to){
    G[++Num]=(edge){Head[from],to},Head[from]=Num; return;
}

il void dfs1(RG int x,RG int p){
    RG int tp=T; fa[x]=p,st[++T]=x;
    for (RG int i=head[x];i;i=g[i].nt) if (g[i].to!=p) dfs1(g[i].to,x);
    if (T-tp>=block){
    ++cnt; while (T>tp) a[++sz[cnt]]=w[st[T]],bl[st[T--]]=cnt;
    sort(a+1,a+sz[cnt]+1); for (RG int j=1;j<=sz[cnt];++j) val[cnt].push_back(a[j]);
    }
    return;
}

il void dfs2(RG int x,RG int p){
    for (RG int i=head[x],v;i;i=g[i].nt){
    v=g[i].to; if (v==p) continue; dfs2(v,x);
    if (bl[x]!=bl[v]) top[bl[v]]=v,link(bl[x],bl[v]);
    }
    return;
}

il int find(RG int x,RG int k){
    RG int ans=0,l=0,r=sz[x]-1,mid;
    while (l<=r){
    mid=(l+r)>>1;
    if (val[x][mid]<=k) ans=mid+1,l=mid+1;
    else r=mid-1;
    }
    return sz[x]-ans;
}

il int calc1(RG int x,RG int k){
    RG int res=find(x,k);
    for (RG int i=Head[x];i;i=G[i].nt) res+=calc1(G[i].to,k);
    return res;
}

il int calc2(RG int x,RG int p,RG int k){
    RG int res=w[x]>k,v;
    for (RG int i=head[x];i;i=g[i].nt){
    v=g[i].to; if (v==p) continue;
    if (bl[x]==bl[v]) res+=calc2(v,x,k);
    else res+=calc1(bl[v],k);
    }
    return res;
}

il int query(RG int x,RG int k){
    RG int ans=0;
    if (top[bl[x]]==x) ans=calc1(bl[x],k); else ans=calc2(x,fa[x],k);
    return ans;
}

il void work(){
    n=gi(); for (RG int i=1,u,v;i<n;++i) u=gi(),v=gi(),insert(u,v),insert(v,u);
    for (RG int i=1;i<=n;++i) w[i]=gi(); block=245,dfs1(1,0);
    if (T){
    top[++cnt]=1; while (T) a[++sz[cnt]]=w[st[T]],bl[st[T--]]=cnt;
    sort(a+1,a+sz[cnt]+1); for (RG int i=1;i<=sz[cnt];++i) val[cnt].push_back(a[i]);
    }
    dfs2(1,0),m=gi();
    while (m--){
    RG int opt=gi(),x=gi()^ans,k=gi()^ans;
    if (!opt) ans=query(x,k),printf("%d
",ans);
    if (opt==1){
        RG int B=bl[x],ccnt=0;
        for (RG int i=0;i<sz[B];++i)
        if (val[B][i]!=w[x]) a[++ccnt]=val[B][i];
        a[++ccnt]=w[x]=k,val[B].clear(),sort(a+1,a+ccnt+1);
        for (RG int i=1;i<=ccnt;++i) val[B].push_back(a[i]);
    }
    if (opt==2){
        insert(x,++n),fa[n]=x,w[n]=k;
        if (sz[bl[x]]>=block){
        top[++cnt]=n,sz[cnt]=1,link(bl[x],cnt);
        bl[n]=cnt,val[cnt].push_back(k);
        } else{
        RG int B=bl[x],ccnt=0; bl[n]=B,a[++ccnt]=w[n];
        for (RG int i=0;i<sz[B];++i) a[++ccnt]=val[B][i];
        val[B].clear(),sort(a+1,a+(++sz[B])+1);
        for (RG int i=1;i<=sz[B];++i) val[B].push_back(a[i]);
        }
    }
    }
    return;
}

int main(){
    File("gty");
    work();
    return 0;
}