bzoj Gty的超级妹子树 块状树

Gty的超级妹子树

Time Limit: 7 Sec  Memory Limit: 32 MB
Submit: 500  Solved: 122
[Submit][Status][Discuss]

Description

我曾在青山之中遇过你，

新竹做杖，鬓插紫茱萸。

跣足踏过无边丝雨，

又拾起燕川雪片片落如席……

Gty神(xian)犇(chong)从来不缺妹子……

他又来到了一棵妹子树下，发现每个妹子有一个美丽度……

由于Gty很 哲♂学 也很 机♂智，他只对美丽度大于某个值的妹子感兴趣。

他想知道某个子树中美丽度大于x的妹子个数。

某个妹子的美丽度可能发生变化……

树上可能会出现一只新的妹子……

但是……树枝可能会断裂，于是，Gty惊讶地发现，他的面前变成了一片妹子树组成的森林……

维护一棵初始有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)

3 u            删除节点u与其父节点之间的路径。此时u的父节点变成叶子节点，u变成分裂出的树的根。(u^=lastans)

最开始时lastans=0。

Input

输入第一行包括一个正整数n(1<=n<=100000)，代表树上的初始节点数。

接下来n-1行，每行2个整数u,v，为树上的一条无向边。

任何时刻，树上的任何权值大于等于0，且两两不同。

接下来1行，包括n个整数wi，表示初始时每个节点的权值。

接下来1行，包括1个整数m(1<=m<=100000)，表示操作总数。

接下来m行，每行最开始包括一个整数op,

若op=3，该行还会有一个整数u；

若op不等于3，该行还会有两个整数u,x；

op,u,x的范围见题目描述。

保证题目涉及的所有数在int内。

Output

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

Sample Input

2
1 2
10 20
1
0 1 5

Sample Output

2

HINT

数据范围见描述。

Source

题解：分块的方式真的不行，被菊花图卡的飞起，这里加了一个删除操作，就是将一个块分成两个，如果刚好一个块的话就直接分出，其实

　　　貌似改的地方不多。

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <cmath>
 #include <cstdlib>
 #include <algorithm>
 #include <vector>
 #define N 200006
 #define M 400006

 using namespace std;
 inline int read(){
     int ret=;char ch=getchar();
     while (ch<''||ch>'') ch=getchar();
     while (''<=ch&&ch<=''){
         ret=ret*-+ch;
         ch=getchar();
     }
     return ret;
 }

 const int KK=;
 vector<int> seq[N];

 struct edge{
     int adj,next;
     edge(){}
     edge(int _adj,int _next):adj(_adj),next(_next){}
 } e[M];
 int n,g[N],m;
 void AddEdge(int u,int v){
     e[++m]=edge(v,g[u]);g[u]=m;
     e[++m]=edge(u,g[v]);g[v]=m;
 }
 int w[N];

 vector<int> nxt[N];
 int fa[N],bl[N],sz[N],cnt,top[N];
 void refresh(int u,int p){
     for (;p>&&seq[u][p]<seq[u][p-];--p)swap(seq[u][p],seq[u][p-]);
     for (;p<sz[u]-&&seq[u][p]>seq[u][p+];++p)swap(seq[u][p],seq[u][p+]);
 }
 void newnode(int u){
     if (sz[bl[fa[u]]]==KK){
         sz[bl[u]=++cnt]=;
         top[cnt]=u;
     }
     else bl[u]=bl[fa[u]];
     ++sz[bl[u]];
 }
 void dfs(int u){
     newnode(u);
     for (int i=g[u];i;i=e[i].next){
         int v=e[i].adj;
         if (v==fa[u]) continue;
         fa[v]=u;
         dfs(v);
     }
 }
 void dfs_pre(int u){
     seq[bl[u]].push_back(w[u]);++sz[bl[u]];
     for (int i=g[u];i;i=e[i].next)if (fa[e[i].adj]==u){
         int v=e[i].adj;
         if (bl[v]==bl[u]) dfs_pre(v);
         else nxt[bl[u]].push_back(bl[v]);
     }
 }
 void prepare(int id){
     seq[id].clear();vector<int>(seq[id]).swap(seq[id]);
     nxt[id].clear();vector<int>(nxt[id]).swap(nxt[id]);
     sz[id]=;
     dfs_pre(top[id]);
     sort(seq[id].begin(),seq[id].end());
 }

 int solve(int u,int lmt){
     int l=-,r=sz[u],mid;
     while (l+<r){
         mid=l+r>>;
         if (seq[u][mid]<=lmt) l=mid;
         else r=mid;
     }
     int ret=sz[u]-l-;
     for (int j=nxt[u].size()-;j>=;--j) ret+=solve(nxt[u][j],lmt);
     return ret;
 }
 int query(int u,int lmt){
     if (bl[u]!=bl[fa[u]]) return solve(bl[u],lmt);
     int ret=w[u]>lmt;
     for (int i=g[u];i;i=e[i].next)if (fa[e[i].adj]==u)ret+=query(e[i].adj,lmt);
     return ret;
 }

 void modify(int u,int x){
     int p;
     for (int i=;i<sz[bl[u]];++i)
         if (seq[bl[u]][i]==w[u]){p=i;break;}
     w[u]=x;seq[bl[u]][p]=x;
     refresh(bl[u],p);
 }

 void create(int u,int x){
     AddEdge(u,++n);fa[n]=u;w[n]=x;
     int tmp=cnt;
     newnode(n);
     seq[bl[n]].push_back(w[n]);
     refresh(bl[n],sz[bl[n]]-);
     if (tmp<cnt) nxt[bl[u]].push_back(cnt);
 }

 void dfs_update(int u,int last,int now){
     bl[u]=now;
     for (int i=g[u];i;i=e[i].next)if (fa[e[i].adj]==u)
         if (bl[e[i].adj]==last) dfs_update(e[i].adj,last,now);
 }

 void cut(int u){
     int lst=bl[u];
     if (bl[fa[u]]!=bl[u]){
         int f=bl[fa[u]];
         for (vector<int>::iterator it=nxt[f].begin();it!=nxt[f].end();++it)
             if ((*it)==lst){nxt[f].erase(it);break;}
         fa[u]=;return;
     }
     top[++cnt]=u;
     fa[u]=;
     dfs_update(u,lst,cnt);
     prepare(lst);
     prepare(cnt);
 }

 int main(){
     n=read();
     memset(g,,sizeof(g));m=;
     for (int i=;i<n;++i) AddEdge(read(),read());
     for (int i=;i<=n;++i) w[i]=read();
     fa[]=;bl[]=;sz[]=KK;cnt=;
     dfs();
     int lastans=;
     for (int i=;i<=cnt;++i) prepare(i);
     for (int Q=read();Q;Q--){
         int op=read(),u=read()^lastans,x;
         if (op<) x=read()^lastans;
         else cut(u);
         if (op==) printf("%d\n",lastans=query(u,x));
         else if (op==) modify(u,x);
         else if (op==) create(u,x);
     }
 }