【块状树】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[20],Num;char c;
inline int Ge()
{
    res=0;c='*';
    while(c<'0'||c>'9')c=getchar();
    while(c>='0'&&c<='9'){res=res*10+(c-'0');c=getchar();}
    return res;
}
inline void P(int x)
{
    if(!x){putchar('0');puts("");return;}Num=0;
    while(x>0){CH[++Num]=x%10;x/=10;}
    while(Num)putchar(CH[Num--]+48);
    putchar('\n');
}
struct Graph
{
    int v[maxn<<1],first[maxn<<1],next[maxn<<1],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=0;
    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+1) , 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+1) , val );
}
int main()
{
    n=Ge();
    for(int i=1;i<n;i++)
      {
          x=Ge();y=Ge();
          G.AddEdge(x,y);
          G.AddEdge(y,x);
      }
    sz=sqrt((double)n*log2(n));
    for(int i=1;i<=n;i++)
      {
          w[i]=Ge();
          top[i]=i;
          siz[i]=1;
      }
    makeblock(1);
    for(int i=1;i<=n;i++)
      if(top[i]==i)
        {
          List[i].push_back(w[i]);
          sort(List[i].begin(),List[i].end());
        }
    makeGoto(1);
    m=Ge();
    for(int i=1;i<=m;i++)
      {
          op=Ge();U=Ge();val=Ge();U^=ans;val^=ans;
          if(!op)query();
          else if(op==1)update();
          else AddPoint();
      }
    return 0;
}