bzoj3510 首都（lct维护虚边信息+树的重心性质）

首先关于树的重心的一些性质：
1.树的重心就是一个点，其所有的子树中最大的子树节点数最少。
2.根据树的重心，每棵子树的大小不超过n的一半。
3.树中所有点到某个点的距离和中，到重心的距离和是最小的，如果有两个距离和，他们的距离和一样。
4.把两棵树通过某两点相连得到一颗新的树，新的树的重心必然在连接原来两棵树重心的路径上。
5.一棵树添加或者删除一个节点，树的重心最多只移动一条边的位置。
性质3可以通过性质2来证明。性质5也可以由性质2得到。一棵树最多有两个重心

因此这题就是要动态的维护每棵树的重心。我们合并两棵树时，可以启发式合并，把sz较小的那棵树拆掉，一个个的扔进大树里。这样相当于每次只添加了一个叶子节点，我们可以直接维护重心是否改变。（记新加的叶子节点所在子树为x，则当子树x的大小*2超过了总结点数时，重心变为x，否则不变。或者相当，取编号较小者。）因此我们就要维护子树大小。可以lct维护虚边信息来做。具体见bzoj4530
复杂度$O(nlog^2n)$

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>
using namespace std;
#define inf 0x3f3f3f3f
#define ll long long
#define N 100010
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,ans,pa[N],G[N],fa[N],c[N][2],sz[N],isz[N],szp[N],h[N],num=0;
bool rev[N];int q[N];
struct edge{
    int to,next;
}data[N<<1];
inline int find(int x){return x==pa[x]?x:pa[x]=find(pa[x]);}
inline void add(int x,int y){
    data[++num].to=y;data[num].next=h[x];h[x]=num;
    data[++num].to=x;data[num].next=h[y];h[y]=num;
}
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];
    sz[p]=sz[l]+sz[r]+isz[p]+1;
}
inline void pushdown(int p){
    if(!rev[p]) return;rev[p]=0;swap(c[p][0],c[p][1]);
    int l=c[p][0],r=c[p][1];rev[l]^=1;rev[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);
            else rotate(y);
        }rotate(x);
    }
}
inline void access(int x){
    int y=0;while(x){splay(x);isz[x]+=sz[c[x][1]]-sz[y];c[x][1]=y;update(x);y=x;x=fa[x];}
}
inline void makeroot(int x){
    access(x);splay(x);rev[x]^=1;
}
inline void link(int x,int y){
    makeroot(y);fa[x]=y;isz[y]+=sz[x];
    int rt=G[find(y)];makeroot(rt);access(x);splay(rt);
    int xx=c[rt][1],size=sz[rt];pushdown(xx);
    while(c[xx][0]) xx=c[xx][0],pushdown(xx);
    access(xx);
    if((isz[xx]+1<<1)>size||(isz[xx]+1<<1)==size&&xx<G[find(y)])
        ans^=G[find(y)],G[find(y)]=xx,ans^=xx;
}
inline void dfs(int x,int Fa){
    isz[x]=fa[x]=c[x][1]=c[x][0]=0;sz[x]=1;link(x,Fa);
    for(int i=h[x];i;i=data[i].next){
        int y=data[i].to;if(y==Fa) continue;dfs(y,x);
    }
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    for(int i=1;i<=n;++i) pa[i]=i,szp[i]=1,G[i]=i,ans^=G[i];
    while(m--){
        char op[5];scanf("%s",op+1);
        if(op[1]=='X'){printf("%d\n",ans);continue;}
        if(op[1]=='Q'){printf("%d\n",G[find(read())]);continue;}
        int x=read(),y=read(),xx=find(x),yy=find(y);
        if(szp[xx]>szp[yy]) swap(xx,yy),swap(x,y);
        pa[xx]=yy;szp[yy]+=szp[xx];ans^=G[xx];dfs(x,y);add(x,y);
    }return 0;
}