[bzoj3720][非树分块算法]Gty的妹子树

前言

这道题我博客有了，最近用了非树分块做法去弄，感谢dwjshift大神教会我正确的姿势。

分块

没有修改和加点操作我们是可以用可持久化线段树合并做的。
对于修改和加点操作进行分块，每跨过一块就暴力重构可持久化线段树。
对于一个询问操作，如果询问是块内的新加点，则暴力扫其整个子树获得答案。
否则，先在线段树内找到没有该块内操作影响时的答案，然后枚举块内的操作，并统计影响。这个可以用时间戳来判断是否造成影响。
然后就是这样了。

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
using namespace std;
const int maxn=30000+10,maxtot=2000000,maxd=2147483647,c=3000;
int root[maxn*2],sum[maxtot],left[maxtot],right[maxtot];
int f[maxn*2],g[maxn*2],fa[maxn*2],fat[maxn*2],h[maxn*2],go[maxn*2],next[maxn*2],a[maxn*2],b[maxn*2];
int sta[c+10][4];
//0 type(1 change 0 add) 1 point 2 new data 3 old data(type=1)
bool bz[maxn*2];
int i,j,k,l,t,n,m,cnt,tot,top,dian,euler,ans,x;
bool czy;
void add(int x,int y){
    fa[y]=x;
    go[++cnt]=y;
    next[cnt]=h[x];
    h[x]=cnt;
}
int newnode(int x){
    sum[++tot]=sum[x];
    left[tot]=left[x];
    right[tot]=right[x];
    return tot;
}
int merge(int a,int b,int l,int r){
    if (!a||!b) return a+b;
    if (l==r){
        int c=newnode(a);
        sum[c]+=sum[b];
        return c;
    }
    int mid=(l+r)/2;
    int c=newnode(a);
    left[c]=merge(left[a],left[b],l,mid);
    right[c]=merge(right[a],right[b],mid+1,r);
    sum[c]=sum[left[c]]+sum[right[c]];
    return c;
}
void insert(int &x,int l,int r,int a){
    x=newnode(x);
    sum[x]++;
    if (l==r) return;
    int mid=(l+r)/2;
    if (a<=mid) insert(left[x],l,mid,a);else insert(right[x],mid+1,r,a);
}
int query(int x,int l,int r,int a,int b){
    if (a>b) return 0;
    if (l==a&&r==b) return sum[x];
    int mid=(l+r)/2;
    if (b<=mid) return query(left[x],l,mid,a,b);
    else if (a>mid) return query(right[x],mid+1,r,a,b);
    else return query(left[x],l,mid,a,mid)+query(right[x],mid+1,r,mid+1,b);
}
void dfs(int x){
    root[x]=0;
    bz[x]=1;
    f[x]=++euler;
    int t=h[x];
    while (t){
        dfs(go[t]);
        root[x]=merge(root[x],root[go[t]],0,maxd);
        t=next[t];
    }
    insert(root[x],0,maxd,a[x]);
    g[x]=++euler;
}
void dg(int x,int y){
    if (b[x]>y) ans++;
    int t=h[x];
    while (t){
        dg(go[t],y);
        t=next[t];
    }
}
void rebuild(){
    int i;
    fo(i,1,dian) a[i]=b[i];
    tot=top=euler=0;
    dfs(1);
}
bool isfa(int x,int y){
    return f[x]<=f[y]&&g[x]>=g[y];
}
int main(){
    freopen("data7.in","r",stdin);freopen("gty.out","w",stdout);
    czy=1;
    scanf("%d",&n);
    dian=n;
    fo(i,1,n-1){
        scanf("%d%d",&j,&k);
        if (j>k) swap(j,k);
        add(j,k);
    }
    fo(i,1,n) scanf("%d",&a[i]),b[i]=a[i];
    rebuild();
    scanf("%d",&m);
    while (m--){ 
        scanf("%d%d%d",&t,&j,&x);
        if (czy) j^=ans,x^=ans;
        if (t==0){
            if (bz[j]){
                ans=query(root[j],0,maxd,x+1,maxd);
                fo(i,1,top){
                    k=sta[i][1];
                    if (isfa(j,k)){
                        if (sta[i][2]>x) ans++;
                        if (sta[i][3]>x) ans--;
                    }
                }
            }
            else{
                ans=0;
                dg(j,x);
            }
            printf("%d\n",ans);
        }
        else{
            if (top==c) rebuild();
            sta[++top][2]=x;
            if (t==1){
                if (bz[j]) sta[top][1]=j;else sta[top][1]=fat[j];
                sta[top][0]=1;
                sta[top][3]=b[j];
                b[j]=x;
            }
            else{
                ++dian;
                if (bz[j]) fat[dian]=j;else fat[dian]=fat[j];
                sta[top][0]=0;
                sta[top][1]=fat[dian];
                sta[top][3]=0;
                b[dian]=a[dian]=x;
                add(j,dian);
            }
        }
    }
    fclose(stdin);fclose(stdout);
}