传送门

Description

「搞 OI 不如种田。」

小 D 在家种了一棵二叉树，第 ii 个结点的权值为 \(a_i\)。

小 D 为自己种的树买了肥料，每天给树施肥。

可是几天后，小 D 却发现树上有几个结点枯死了，他这才发现，自己买的肥料是二叉搜索树专用版。

二叉搜索树是一种二叉树，满足每个结点的权值大于左子树内所有点的权值，小于右子树内所有点的权值。

二叉搜索树专用版肥料是这么工作的：首先，假设所有节点权值互不相同（小 D 的二叉树可能不满足），每种权值对应一种肥料，所有肥料会从根进入树中，如果一种肥料对应的权值等于当前结点权值，这种肥料会被当前结点完全吸收，否则若肥料对应的权值小于当前结点权值，肥料会流向左子树，否则流向右子树，如果流向的子树为空，肥料只好流失蒸发了。显然，如果树是二叉搜索树，所有节点都能吸收到肥料。

小 D 觉得自己还能抢救一下，他会进行若干次操作，每次操作修改一个点的权值或者翻转一个子树（子树内所有节点左右儿子互换）。在操作过程中，他有时会想知道一个点当前是否能吸收到肥料，以决定之后如何操作，请你帮帮可怜的小 D 吧。

Solution

我们把每条边当成是一个元素，这个元素是一个可行的区间

对于一个点，它能够可行当且仅当它属于它到根节点路径上所有元素区间的并

可以求出dfs序后用线段树维护

对于取反操作，我们同时维护一下线段树上每个元素取反后的区间，修改时直接打标记就行

Code

#include<bits/stdc++.h>
#define ll long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
inline int read()
{int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}return x*f;
}
const int MN=1e6+5,inf=1<<30;
int n,m,a[MN],c[MN][2],dfn[MN],fdfn[MN],top[MN],siz[MN],mx[MN],fa[MN],R[MN];
bool ch[MN];
inline void dfs1(int x,int f)
{siz[x]=1;fa[x]=f;ch[x]=c[f][1]==x;if(c[x][0]) dfs1(c[x][0],x);if(c[x][1]) dfs1(c[x][1],x);siz[x]+=siz[c[x][0]]+siz[c[x][1]];mx[x]=siz[c[x][0]]>siz[c[x][1]]?0:1;
}
inline void dfs2(int x,int f,int tp)
{static int ind=0;dfn[x]=++ind;top[x]=tp;fdfn[ind]=x;if(c[x][mx[x]]) dfs2(c[x][mx[x]],x,tp);if(c[x][mx[x]^1]) dfs2(c[x][mx[x]^1],x,c[x][mx[x]^1]);R[x]=ind;
}
class Test
{bool rev[MN<<2],son[MN<<1];struct F{int l,r,_l,_r;}T[MN<<2];F merge(F o,F oo){return (F){max(o.l,oo.l),min(o.r,oo.r),max(o._l,oo._l),min(o._r,oo._r)};}bool P(int x,F y){return x>=y.l&&x<=y.r;}void Rev(int x){rev[x]^=1;son[x]^=1;std::swap(T[x].l,T[x]._l);std::swap(T[x].r,T[x]._r);}void pushdown(int x){if(!rev[x])return;Rev(x<<1);Rev(x<<1|1);rev[x]=0;}#define L fdfn[l]void update(int x,int l,int r,int u){if(l==r){if(son[x]) T[x].l=a[fa[L]]+1,T[x].r=inf,T[x]._l=1,T[x]._r=a[fa[L]]-1;else T[x]._l=a[fa[L]]+1,T[x]._r=inf,T[x].l=1,T[x].r=a[fa[L]]-1;return;}int mid=(l+r)>>1;pushdown(x);if(u<=mid) update(x<<1,l,mid,u);else update(x<<1|1,mid+1,r,u);T[x]=merge(T[x<<1],T[x<<1|1]);}void re(int x,int l,int r,int a,int b){if(a>b) return;if(l==a&&r==b){Rev(x);return;}int mid=(l+r)>>1;pushdown(x);if(b<=mid) re(x<<1,l,mid,a,b);else if(a>mid) re(x<<1|1,mid+1,r,a,b);else re(x<<1,l,mid,a,mid),re(x<<1|1,mid+1,r,mid+1,b);T[x]=merge(T[x<<1],T[x<<1|1]);}F Query(int x,int l,int r,int a,int b){if(a==l&&b==r) return T[x];int mid=(l+r)>>1;pushdown(x);if(b<=mid) return Query(x<<1,l,mid,a,b);else if(a>mid) return Query(x<<1|1,mid+1,r,a,b);return merge(Query(x<<1,l,mid,a,mid),Query(x<<1|1,mid+1,r,mid+1,b));}public:void Build(int x,int l,int r){if(l==r){son[x]=ch[L];if(L==1) T[x].l=T[x]._l=1,T[x].r=T[x]._r=inf;else if(ch[L]) T[x].l=a[fa[L]]+1,T[x].r=inf,T[x]._l=1,T[x]._r=a[fa[L]]-1;else T[x]._l=a[fa[L]]+1,T[x]._r=inf,T[x].l=1,T[x].r=a[fa[L]]-1;return;}int mid=(l+r)>>1;Build(x<<1,l,mid);Build(x<<1|1,mid+1,r);T[x]=merge(T[x<<1],T[x<<1|1]);}void upd(int x,int y){a[x]=y;if(c[x][0]) update(1,1,n,dfn[c[x][0]]);if(c[x][1]) update(1,1,n,dfn[c[x][1]]);}bool query(int x){int o=a[x];F ans=(F){1,inf,1,inf};while(x&&ans.l<=ans.r) ans=merge(ans,Query(1,1,n,dfn[top[x]],dfn[x])),x=fa[top[x]];return P(o,ans);}void reve(int x){re(1,1,n,dfn[x]+1,R[x]);}#undef L
}T;
int main()
{n=read();m=read();register int i,opt,x;for(i=1;i<=n;++i) a[i]=read(),c[i][0]=read(),c[i][1]=read();dfs1(1,0);dfs2(1,0,1);T.Build(1,1,n);while(m--){opt=read();x=read();if(opt==1) T.upd(x,read());if(opt==2) T.reve(x);if(opt==3) puts(T.query(x)?"YES":"NO");}return 0;
}

Blog来自PaperCloud，未经允许，请勿转载，TKS！