bzoj3282&&luogu3690【模板】Link Cut Tree (动态树)

觉得这篇blog讲得不错:传送门
模板题。lct,动态树,树的形态可以变化的树链剖分,用splay来维护每条重链上的信息。

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 300010
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*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,v[N],sum[N],fa[N],c[N][2],q[N];bool rev[N];
inline bool isroot(int x){
    return x!=c[fa[x]][0]&&x!=c[fa[x]][1];
}
inline void update(int p){
    int l=c[p][0],r=c[p][1];
    sum[p]=sum[l]^v[p]^sum[r];
}
inline void pushdown(int p){
    if(!rev[p]) return;rev[p]=0;
    int l=c[p][0],r=c[p][1];
    rev[l]^=1;swap(c[l][0],c[l][1]);
    rev[r]^=1;swap(c[r][0],c[r][1]);
}
inline void rotate(int x){
    int y=fa[x],z=fa[y],l=x==c[y][1],r=l^1;
    if(!isroot(y)) c[z][y==c[z][1]]=x;
    fa[c[x][r]]=y;fa[y]=x;fa[x]=z;
    c[y][l]=c[x][r];c[x][r]=y;update(y);update(x);
}
inline void splay(int x){
    int top=0;q[++top]=x;
    for(int xx=x;!isroot(xx);xx=fa[xx]) q[++top]=fa[xx];
    while(top) pushdown(q[top--]);
    while(!isroot(x)){
        int y=fa[x],z=fa[y];
        if(!isroot(y)){
            if(x==c[y][1]^y==c[z][1]) rotate(x);
            rotate(y);
        }rotate(x);
    }
}
inline void access(int x){
    int y=0;
    while(x){
        splay(x);c[x][1]=y;update(x);
        y=x;x=fa[x];
    }
}
inline void makeroot(int x){
    access(x);splay(x);rev[x]^=1;swap(c[x][0],c[x][1]);
}
inline int find(int x){
    access(x);splay(x);
    while(c[x][0]) x=c[x][0];return x;
}
inline void query(int x,int y){
    makeroot(x);access(y);splay(y);
    printf("%d\n",sum[y]);
}
inline void dolink(int x,int y){
    int rt1=find(x),rt2=find(y);
    if(rt1==rt2) return;
    makeroot(x);fa[x]=y;
}
inline void cut(int x,int y){
    makeroot(x);access(y);splay(y);
    if(c[y][0]==x){//???
        c[y][0]=fa[x]=0;update(y);
    }
}
inline void change(int x,int val){
    splay(x);v[x]=val;update(x);
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    for(int i=1;i<=n;++i) v[i]=read();
    while(m--){
        int op=read(),x=read(),y=read();
        if(op==0) query(x,y);
        else if(op==1) dolink(x,y);
        else if(op==2) cut(x,y);
        else change(x,y);
    }return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值