【bzoj3786】星系探索

ETT模版题。

真正的Eular-Tour-Tree维护的是树的欧拉序。

由于各种原因，没人知道怎么维护欧拉序，所以我写的是个假的，维护dfs序的。

本质还是用Splay维护序列。

然后因为我常数太差，压着线跑过去的……

#include<bits/stdc++.h>
#define N 100010
using namespace std;
typedef long long ll;
const int M=;
int n,m,a[M];
struct Edge{int u,v,next;}G[M];
int head[M],tot=,tpos[M],rt;
inline void addedge(int u,int v){
    G[++tot].u=u;G[tot].v=v;G[tot].next=head[u];head[u]=tot;
    G[++tot].u=v;G[tot].v=u;G[tot].next=head[v];head[v]=tot;
}
int size[M],c[M][],val[M],addv[M],fa[M],w[M],s[M],top;
ll sumv[M];
inline void pushup(int x){
    size[x]=size[c[x][]]+size[c[x][]]+w[x];
    sumv[x]=sumv[c[x][]]+sumv[c[x][]]+w[x]*val[x];
}
inline void puttag(int x,ll v){val[x]+=v;addv[x]+=v;sumv[x]+=v*size[x];}
inline void pushdown(int x){
    if(!addv[x])return ;
    if(c[x][])puttag(c[x][],addv[x]);
    if(c[x][])puttag(c[x][],addv[x]);
    addv[x]=;
}
inline void rotate(int x,int &k){
    int y=fa[x],z=fa[y],l,r;
    if(c[y][]==x)l=;else l=;r=l^;
    if(y==k)k=x;else{if(c[z][]==y)c[z][]=x;else c[z][]=x;}
    fa[x]=z;fa[y]=x;fa[c[x][r]]=y;
    c[y][l]=c[x][r];c[x][r]=y;
    pushup(y);pushup(x);
}
inline void splay(int x,int &k){
    s[top=]=x;for(int i=x;i!=rt;i=fa[i])s[++top]=fa[i];
       for(int i=top;i;--i)pushdown(s[i]);
       while(x!=k){
           int y=fa[x],z=fa[y];
           if(y!=k)
               if(c[z][]==y^c[y][]==x)rotate(x,k);
               else rotate(y,k);
           rotate(x,k);
       }
}
inline void ins(int x,int v,int ww){
    if(!rt)rt=x;else fa[x]=rt,c[rt][]=x;
    val[x]=v;w[x]=ww;pushup(x);splay(x,rt);
}
inline int tmax(int x){while(c[x][])x=c[x][];return x;}
inline int tmin(int x){while(c[x][])x=c[x][];return x;}
inline ll query(int x){
    splay(x,rt);return sumv[c[x][]]+1LL*w[x]*val[x];
}
inline void change(int x,int y){
    splay(x,rt);int t1=tmax(c[rt][]);
    splay(x+N,rt);int t2=tmin(c[rt][]);
    splay(t1,rt);splay(t2,c[t1][]);
    int z=c[t2][];fa[z]=c[t2][]=;
    splay(y,rt);int t3=tmin(c[y][]);
    splay(t3,c[y][]);fa[z]=t3;c[t3][]=z;
}
inline void add(int x,int v){
    splay(x,rt);int t1=tmax(c[rt][]);
       splay(x+N,rt);int t2=tmin(c[rt][]);
       splay(t1,rt);splay(t2,c[t1][]);
       puttag(c[t2][],v);
}
void build(int u,int f){
    ins(u+,a[u],);
    for(int i=head[u];i;i=G[i].next){
        int v=G[i].v;if(v==f)continue;
        build(v,u);
    }
    ins(u+N+,a[u],-);
}
int main(){
    n=io.read();
    for(int i=;i<=n;i++){int u=io.read();addedge(u,i);}
    for(int i=;i<=n;i++)a[i]=io.read();
    ins(,,);build(,);ins(N*+,,);
    m=io.read();char s[];
    while(m--){
        scanf("%s",s);
        if(s[]=='Q'){int x=read();printf("%lld\n",query(x+));}
        if(s[]=='C'){int x=read(),y=read();change(x+,y+);}
        if(s[]=='F'){int x=read(),y=read();add(x+,y);}
    }
    return ;
}