bzoj 3720 Gty的妹子树 树分块？瞎搞

Gty的妹子树

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 2149  Solved: 781
[Submit][Status][Discuss]

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

HINT

2017.9.28新加数据一组by GXZlegend,未重测

Source

传说中的树分块，但是会被暴力卡，我貌似没找出什么方法，

和根一起一起到达根号n就就分块，这个没什么用，会被菊花图卡掉的。

 #include<cstring>
 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #include<cmath>

 #define N 60007
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(!isdigit(ch)){if(ch=='-')f=-;ch=getchar();}
     while(isdigit(ch)){x=(x<<)+(x<<)+ch-'';ch=getchar();}
     return x*f;
 }

 int n,m,block,ans,tot;
 int a[N],fa[N],bl[N];
 int cnt,hed[N],nxt[N<<],rea[N<<];
 int Cnt,Hed[N],Nxt[N<<],Rea[N<<];

 struct Node
 {
     int a[],siz;
     void insert(int x){a[++siz]=x;}
     void change(int x,int y)
     {
         int t=upper_bound(a+,a+siz+,x)-a;t--;
         a[t]=y;
         sort(a+,a+siz+);
     }
     int query(int x)
     {
         int t=upper_bound(a+,a+siz+,x)-a;
         return siz-t+;
     }
 }blo[N];
 inline void add(int u,int v)
 {
     nxt[++cnt]=hed[u];
     hed[u]=cnt;
     rea[cnt]=v;
 }
 inline void add_two_way(int u,int v)
 {
     add(u,v);
     add(v,u);
 }
 inline void ins(int u,int v)
 {
     Nxt[++Cnt]=Hed[u];
     Hed[u]=Cnt;
     Rea[Cnt]=v;
 }
 void dfs(int u)
 {
     if (blo[bl[fa[u]]].siz==block)
         blo[bl[u]=++tot].insert(a[u]),ins(bl[fa[u]],tot);
     else blo[bl[u]=bl[fa[u]]].insert(a[u]);
     for (int i=hed[u];i!=-;i=nxt[i])
     {
         int v=rea[i];
         if (v==fa[u]) continue;
         fa[v]=u,dfs(v);
     }
 }
 void block_dfs(int u,int z)
 {
     ans+=blo[u].query(z);
     for (int i=Hed[u];i!=-;i=Nxt[i])
     {
         int v=Rea[i];
         block_dfs(v,z);
     }
 }
 void solve(int u,int z)
 {
     if (a[u]>z) ans++;
     for (int i=hed[u];i!=-;i=nxt[i])
     {
         int v=rea[i];
         if (v==fa[u]) continue;
         if (bl[v]==bl[u]) solve(v,z);
         else block_dfs(bl[v],z);
     }
 }
 int main()
 {
     memset(hed,-,sizeof(hed));
     memset(Hed,-,sizeof(Hed));
     n=read(),block=sqrt(n);
     for (int i=;i<n;i++)
         add_two_way(read(),read());
     for (int i=;i<=n;i++) a[i]=read();
     dfs();
     for (int i=;i<=tot;i++)
         sort(blo[i].a+,blo[i].a+blo[i].siz+);

     m=read();
     for (int i=;i<=m;i++)
     {
         int opt=read(),x=read(),y=read();
         x^=ans,y^=ans;
         if (opt==)
         {
             ans=;
             solve(x,y);
             printf("%d\n",ans);
         }
         else if (opt==)
         {
             blo[bl[x]].change(a[x],y);
             a[x]=y;
         }
         else
         {
             a[++n]=y,add_two_way(x,n),fa[n]=x;
             if (blo[bl[x]].siz==block)
                 blo[bl[n]=++tot].insert(y),ins(bl[x],tot);
             else
             {
                 blo[bl[n]=bl[x]].insert(y);
                 sort(blo[bl[n]].a+,blo[bl[n]].a+blo[bl[n]].siz+);
             }
         }
     }
 }