【块状树】bzoj3720 Gty的妹子树

块状树。
教程见：http://z55250825.blog.163.com/blog/static/1502308092014163413858/
将树按一定大小分块，分成许多子树，在每个子树的根节点记录一个数组，存储这个块中数据的有序表，查询时二分即可。

对于添加节点操作，若添加的节点的父节点所在块大小没有超过size，则添加到上面的块里，否则另起一块。

这个size值要注意，由于有二分操作，所以大小定为sqrt(n*log2(n))比较好，而不是sqrt(n)，证明略。

时间复杂度O(n*sqrt(n)*log(sqrt(n)))。

 #include<cstdio>
 #include<algorithm>
 #include<cmath>
 #include<vector>
 using namespace std;
 #define maxn 60001
 typedef vector<int>::iterator ITER;
 vector<int>List[maxn],Goto[maxn];
 int res,CH[],Num;char c;
 inline int Ge()
 {
     res=;c='*';
     while(c<''||c>'')c=getchar();
     while(c>=''&&c<=''){res=res*+(c-'');c=getchar();}
     return res;
 }
 inline void P(int x)
 {
     if(!x){putchar('');puts("");return;}Num=;
     while(x>){CH[++Num]=x%;x/=;}
     while(Num)putchar(CH[Num--]+);
     putchar('\n');
 }
 struct Graph
 {
     int v[maxn<<],first[maxn<<],next[maxn<<],en;
     void AddEdge(const int &a,const int &b)
     {v[++en]=b;next[en]=first[a];first[a]=en;}
 };
 Graph G,son;
 int top[maxn],siz[maxn],sz,w[maxn];
 bool vis[maxn];
 int n,x,y,m,op;
 int ans,val,U;
 void makeblock(int cur)
 {
     vis[cur]=true;
     for(int i=G.first[cur];i;i=G.next[i])
       if(!vis[G.v[i]])
         {
           son.AddEdge(cur,G.v[i]);
           if(siz[top[cur]]<sz)
             {
               List[top[cur]].push_back(w[G.v[i]]);
               siz[top[cur]]++;
               top[G.v[i]]=top[cur];
             }
           makeblock(G.v[i]);
         }
 }
 void makeGoto(int cur)
 {
     for(int i=son.first[cur];i;i=son.next[i])
       {
           if(top[son.v[i]]!=top[cur])
           Goto[top[cur]].push_back(son.v[i]);
         makeGoto(son.v[i]);
       }
 }
 void dfs(int cur)
 {
     ans+=( List[cur].end() - upper_bound( List[cur].begin() , List[cur].end() , val ) );
     for(ITER it=Goto[cur].begin();it!=Goto[cur].end();it++)
       dfs(*it);
 }
 void dfs_block(int cur)
 {
     if(w[cur]>val) ans++;
     for(int i=son.first[cur];i;i=son.next[i])
       if(top[son.v[i]]==top[cur]) dfs_block(son.v[i]);
       else dfs(son.v[i]);
 }
 void query()
 {
     ans=;
     if(U==top[U]) dfs(U);
     else dfs_block(U);
     P(ans);
 }
 void update()
 {
     List[top[U]].erase( lower_bound(List[top[U]].begin(),List[top[U]].end(),w[U]) );
     w[U]=val;
     List[top[U]].insert( lower_bound(List[top[U]].begin(),List[top[U]].end(),val+) , val );
 }
 void AddPoint()
 {
     n++;
     if(siz[top[U]]<sz)
       {
           top[n]=top[U];
           siz[top[n]]++;
       }
     else
       {
           top[n]=n;
           siz[n]++;
           Goto[top[U]].push_back(n);
       }
     son.AddEdge(U,n);
     w[n]=val;
     List[top[n]].insert( lower_bound(List[top[n]].begin(),List[top[n]].end(),val+) , val );
 }
 int main()
 {
     n=Ge();
     for(int i=;i<n;i++)
       {
           x=Ge();y=Ge();
           G.AddEdge(x,y);
           G.AddEdge(y,x);
       }
     sz=sqrt((double)n*log2(n));
     for(int i=;i<=n;i++)
       {
           w[i]=Ge();
           top[i]=i;
           siz[i]=;
       }
     makeblock();
     for(int i=;i<=n;i++)
       if(top[i]==i)
         {
           List[i].push_back(w[i]);
           sort(List[i].begin(),List[i].end());
         }
     makeGoto();
     m=Ge();
     for(int i=;i<=m;i++)
       {
           op=Ge();U=Ge();val=Ge();U^=ans;val^=ans;
           if(!op)query();
           else if(op==)update();
           else AddPoint();
       }
     return ;
 }