AAAT 笔记(P5649)

实际上去掉主函数不长于线段树 3。

对于 LCT 每个点的虚儿子。用 splay 把它们串起来(称为新 splay,虽然是共用的)。

具体来说,设 \(1\le x\le n\) 是原 LCT 的 splay 的根,\(x+n\) 就用来挂载其子树信息(虚实都包括),认 \(x+n\) 做父亲。

而每个有虚儿子的结点 \(x\) 具有 splay 上儿子 \(n<rt_x\le 2n\) 作为新 splay 的一个根来记录其虚儿子子树信息,这个 \(rt_x\) 所在的新 splay 上就有的是它下面虚儿子的 \(y+n\)

每个点维护的信息包括虚子树信息和和实 splay 信息,以及虚实的标记,当然 \(n+1\sim 2n\) 的范围没有虚儿子(也没有虚标记),虚标记传到 \(rt_x\) 的时候也要变成实标记。

虚子树信息需要一道累加上实子树信息。这样就可以写出链、子树修改查询的代码了:

node query(int x,int y){split(x,y);return sr[y];}
node query(int x){split(crt,x);return v[x]+sr[rt[x]];}
void modify(int x,int y,tag t){split(x,y);ADD(y,t,tag());}
void modify(int x,tag t){split(crt,x);ADD(rt[x],t,tag());v[x]=v[x]*t;pushup(x);}

其中 sr 是实信息和,sv 是虚信息和。

至于换根的操作,直接 makeroot 即可。而换儿子只需找到父亲 linkcut 即可。

还有问题:splay 的时候会改变当前 splay 根,因此还需支持新 splay 的添加/删除虚子树。添加很简单,挂成新的根就可以。删除只需找到前驱拼上右儿子即可。splay 前/后需要删除/重新添加这个子树的信息。

access 也涉及改变虚儿子状态,也要类似的删除/添加。

还有一个调了一个下午(这个调试真的麻烦)的问题是 link 和一般 LCT 有所不同。link(x,y) 需要先 access,splay(y) 因为插入虚儿子有更新信息环节。

#include<bits/stdc++.h>
using namespace std;
#define int long long 
#define endl "\n"
const int maxn=4e5+5,INF=1e12;
struct tag{
	int k,b;
	tag(int x=1,int y=0){k=x,b=y;}
}rtag[maxn],vtag[maxn];
struct node{
	int mn,mx,sze,sum;
	node(int x=INF,int y=-INF,int z=0,int w=0){mn=x,mx=y,sze=z,sum=w;}
}v[maxn],sr[maxn],sv[maxn];
int tg[maxn],ch[maxn][2],fa[maxn],a[maxn],st[maxn],tp=0,sze[maxn],rt[maxn],n,q;
tag operator ^(tag A,tag B){return {B.k*A.k,B.k*A.b+B.b};}
node operator *(node A,tag B){
	if(!A.sze)return {INF,-INF,0,0};
	return {B.k*A.mn+B.b,B.k*A.mx+B.b,A.sze,B.k*A.sum+B.b*A.sze};
}
node operator +(node A,node B){return {min(A.mn,B.mn),max(A.mx,B.mx),A.sze+B.sze,A.sum+B.sum};}
#define ls(x) (ch[x][0])
#define rs(x) (ch[x][1])
int nroot(int x){return ls(fa[x])==x||rs(fa[x])==x;}
int whs(int x){return rs(fa[x])==x;}
void pushup(int x){
	if(x>n){
		v[x]=sr[x-n]+sv[x-n];
		sr[x]=v[x]+sr[ls(x)]+sr[rs(x)];
	}else{
		sr[x]=v[x]+sr[ls(x)]+sr[rs(x)];
		sv[x]=sr[rt[x]]+sv[ls(x)]+sv[rs(x)];
		pushup(x+n);
	}
}
void rot(int x){
	int y=fa[x],z=fa[y],k=whs(x),w=ch[x][k^1];
	if(nroot(y))ch[z][whs(y)]=x;fa[y]=x,fa[x]=z;
	ch[y][k]=w;if(w)fa[w]=y;ch[x][k^1]=y;
	pushup(y),pushup(x);
}
void ADD(int x,tag tr,tag tv){
	if(!x)return ;
	v[x]=v[x]*tr,sr[x]=sr[x]*tr,rtag[x]=rtag[x]^tr;
	sv[x]=sv[x]*tv,vtag[x]=vtag[x]^tv;
}
void FLIP(int x){swap(ls(x),rs(x)),tg[x]^=1;}
void pushdown(int x){
	if(tg[x])FLIP(ls(x)),FLIP(rs(x)),tg[x]=0;
	if(x>n){
		ADD(ls(x),rtag[x],tag());
		ADD(rs(x),rtag[x],tag());
		ADD(x-n,rtag[x],rtag[x]);
		rtag[x]=tag();
	}else{
		ADD(ls(x),rtag[x],vtag[x]);
		ADD(rs(x),rtag[x],vtag[x]);
		ADD(rt[x],vtag[x],tag());
		rtag[x]=tag(),vtag[x]=tag();
	}
}
void del(int x,int y);
void ins(int x,int y);
void splay(int x);
void ins(int x,int y){
	if(y<=n)return ;
	if(!rt[x])pushup(y),rt[x]=y,pushup(x);
	else{
		int p=rt[x];
		ls(y)=p,fa[p]=y,rt[x]=y;
		pushup(y);pushup(x);
	}
}
void del(int x,int y){
	if(y<=n)return ;
	splay(y),fa[ls(y)]=fa[rs(y)]=0;
	if(ls(y)){
		int p=ls(y);while(rs(p))pushdown(p),p=rs(p);
		splay(p);rs(p)=rs(y),fa[rs(y)]=p;
		pushup(p),rt[x]=p;
	}else rt[x]=rs(y);
	ls(y)=rs(y)=0;
	pushup(y),pushup(x);
}
void splay(int x){
	int y=x;while(nroot(y))y=fa[y];
	if(y<=n&&fa[y])del(fa[y],y+n);
	y=x;st[tp=1]=y;while(nroot(y))y=fa[y],st[++tp]=y;
	while(tp)pushdown(st[tp--]);
	while(nroot(x)){
		int y=fa[x];
		if(nroot(y))rot(whs(y)==whs(x)?y:x);
		rot(x);
	}	
	if(x<=n&&fa[x])ins(fa[x],x+n);
}
int crt,mkc=0;
void access(int x){for(int y=0;x;y=x,x=fa[x])splay(x),del(x,y+n),ins(x,rs(x)+n),rs(x)=y,pushup(x);}
void makeroot(int x){access(x),splay(x),FLIP(x);}
void split(int x,int y){
	makeroot(x);
	access(y);splay(y);
}
void link(int x,int y){
	if(!mkc)makeroot(x);
	access(y),splay(y);
	fa[x]=y;
	ins(y,x+n);
}
void cut(int x,int y){split(x,y);fa[x]=ch[y][0]=0;pushup(y);pushup(x);}
int findroot(int x){access(x),splay(x);while(ls(x))pushdown(x),x=ls(x);splay(x);return x;}
void cgf(int x,int y){
	split(crt,x);
	int f=ls(x);
	if(!f)return ;
	while(rs(f))pushdown(f),f=rs(f);
	cut(x,f);
	split(x,y);
	if(findroot(y)==x)link(x,f);
	else link(x,y);
}
int U[maxn],V[maxn];
node query(int x,int y){split(x,y);return sr[y];}
node query(int x){
	split(crt,x);
	return v[x]+sr[rt[x]];
}
void modify(int x,int y,tag t){split(x,y);ADD(y,t,tag());}
void modify(int x,tag t){split(crt,x);ADD(rt[x],t,tag());v[x]=v[x]*t;pushup(x);}
int currt(int x){x=findroot(x);if(!fa[x])return x;return currt(fa[x]);}
signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	cin>>n>>q;for(int i=1;i<n;i++)cin>>U[i]>>V[i];
	for(int i=1,x;i<=n;i++)cin>>x,v[i]={x,x,1,x},pushup(i);cin>>crt;
	for(int i=1;i<n;i++)link(U[i],V[i]);
	for(int i=1;i<=q;i++){
		int tp,x,y,z;
		cin>>tp>>x;
		if(tp==0)cin>>y,modify(x,{0,y});
		if(tp==1)crt=x;
		if(tp==2)cin>>y>>z,modify(x,y,{0,z});
		if(tp==3)cout<<query(x).mn<<endl;
		if(tp==4)cout<<query(x).mx<<endl;
		if(tp==5)cin>>y,modify(x,{1,y});
		if(tp==6)cin>>y>>z,modify(x,y,{1,z});
		if(tp==7)cin>>y,cout<<query(x,y).mn<<endl;
		if(tp==8)cin>>y,cout<<query(x,y).mx<<endl;
		if(tp==9)cin>>y,cgf(x,y);
		if(tp==10)cin>>y,cout<<query(x,y).sum<<endl;
		if(tp==11)cout<<query(x).sum<<endl;
	}
	return 0;
}
posted @ 2025-01-05 21:35  British_Union  阅读(40)  评论(1)    收藏  举报