前言
这道题我博客有了,最近用了非树分块做法去弄,感谢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);
}