[BZOJ3731]Gty的超级妹子树（树上分块）

题目：

我是超链接

题解：

因为这题太神而我太sb的缘故，调了好久还是RE？弃了弃了。
说个思路吧，下面的代码估计也就是个手残我就是懒得调了打我啊
断开一个节点与父节点的连边时，如果这个点是所在块的根节点，直接断掉就行
如果这个点不是所在块的根节点，那么就要把这个块分裂，这个点以及在块中的子树都分裂到新的块中
然后是此题卡内存 静态不用想了 开vector吧

代码：

#include <cmath>
#include <vector>
#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
const int N=200005;
const int M=400005;
int nw,tot,nxt[M],point[M],v[M],tot1,nxt1[M],point1[M],v1[M],pos[M],ans,block,cnt,a[M],f[M],visk[M];
bool ml[N];
struct Block
{
    vector<int> a;int size;
    void insert(int x){size++;a.insert(lower_bound(a.begin(),a.end(),x+1),x);}
    void del(int x){size--;a.erase(lower_bound(a.begin(),a.end(),x));}
    void change(int x,int y){del(x); insert(y);}
    int qurry(int x)
    {
        int pos=upper_bound(a.begin(),a.end(),x)-a.begin()+1;
        return a.size()-pos+1;
    }
}b[M];
inline int read(){
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while ('0'<=ch&&ch<='9') x=x*10-'0'+ch,ch=getchar();
    return x;
}
void addline(int x,int y)
{
    ++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y;
    ++tot; nxt[tot]=point[y]; point[y]=tot; v[tot]=x;
}
void addline1(int x,int y){++tot1; nxt1[tot1]=point1[x]; point1[x]=tot1; v1[tot1]=y;}
void dfs(int x,int fa)
{
    if (b[pos[fa]].size==block)
      pos[x]=++cnt,b[cnt].insert(a[x]),addline1(pos[fa],pos[x]),visk[pos[x]]=nw;
    else pos[x]=pos[fa],b[pos[x]].insert(a[x]);
    for (int i=point[x];i;i=nxt[i])
      if (v[i]!=fa) f[v[i]]=x,dfs(v[i],x);
}
void findk(int x,int y)
{
    ans+=b[x].qurry(y);
    for (int i=point1[x];i;i=nxt1[i])
      if (visk[v1[i]]==visk[x]) findk(v1[i],y);
}
void find(int x,int y)
{
    if (a[x]>y) ans++;
    for (int i=point[x];i;i=nxt[i])
      if (v[i]!=f[x] && visk[pos[v[i]]]==visk[pos[x]]) 
      {
        if (pos[v[i]]==pos[x]) find(v[i],y);
        else findk(pos[v[i]],y);//找到这个块和他管辖的块的答案 
      }
}
void push(int x)
{
    for (int i=point1[x];i;i=nxt1[i])
      if (visk[v1[i]]==visk[x]) push(v1[i]);
    visk[x]=nw;
}
void dele(int x,int fa,int cnt,int gd)
{
    b[pos[gd]].del(a[x]);   b[cnt].insert(a[x]);
    for (int i=point[x];i;i=nxt[i])
      if (v[i]!=fa && pos[v[i]]==pos[x] && visk[pos[v[i]]]==visk[pos[x]]) dele(v[i],x,cnt,fa); 
      else if (v[i]!=fa && visk[pos[v[i]]]==visk[pos[x]]) 
        push(pos[v[i]]),addline1(cnt,pos[v[i]]);
}
void dle(int x)
{
    for (int i=point1[x];i;i=nxt1[i])
      if (visk[v1[i]]==visk[x]) push(v1[i]);
    visk[x]=nw;
}
int main()
{
//  freopen("sbg.in","r",stdin);
    int n,m;n=read();
    for (int i=1;i<n;i++)
    {
        int x,y;x=read();y=read();
        addline(x,y);
    }
    for (int i=1;i<=n;i++) a[i]=read();
    block=sqrt(n);nw=1;visk[0]=1;
    dfs(1,0);m=read();
    while(m--)
    {
        int id,x,y;id=read();
        if (id==3)
        {
            x=read();x^=ans;
            if (ml[x]) continue;
            ++nw;ml[x]=1;
            if (pos[x]==pos[f[x]]) dele(x,f[x],cnt+1,f[x]),pos[x]=++cnt;
            else dle(pos[x]);
            visk[pos[x]]=nw;
            f[x]=0;
        }
        else
        {
        x=read();y=read();x^=ans; y^=ans;
        switch(id)
        {
            case 0:ans=0;find(x,y);printf("%d\n",ans);break;
            case 1:b[pos[x]].change(a[x],y);a[x]=y;break;
            case 2:
            {
                a[++n]=y;addline(x,n);f[n]=x;
                if (b[pos[x]].size==block) pos[n]=++cnt,b[cnt].insert(y),addline1(pos[x],pos[n]);
                else pos[n]=pos[x],b[pos[x]].insert(y); 
                visk[pos[n]]=visk[pos[x]];
                break;
            }
        }   
        }   
    }
}