主席树+启发式合并（LT） BZOJ3123

好久没做题了，写道SBT又RE又T

查询：主席树裸题。

修改：对于两个树合并重建小的树。

注意fa[x][i]重新计算时要清空

#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
using namespace std;
inline int read()
{
    char c=getchar();int x=,sig=;
    for(;!isdigit(c);c=getchar()) if(c=='-') sig=-;
    for(;isdigit(c);c=getchar()) x=x*+c-'';
    return x*sig;
}
inline void print(int x)
{
    int buf[],l=;
    if(!x) l++;
    while(x) buf[++l]=x%,x/=;
    for(;l;l--) putchar(buf[l]+'');
    putchar('\n');
}
const int maxn=;
const int maxnode=;
int ls[maxnode],rs[maxnode],s[maxnode],ToT;
void update(int& y,int x,int l,int r,int pos)
{
     s[y=++ToT]=s[x]+;if(l==r) return;
     int mid=l+r>>;ls[y]=ls[x];rs[y]=rs[x];
     if(pos<=mid) update(ls[y],ls[x],l,mid,pos);
     else update(rs[y],rs[x],mid+,r,pos);
}
int query(int x1,int x2,int y1,int y2,int l,int r,int k)
{
    if(l==r) return l;
    int mid=l+r>>,k2=s[ls[y1]]+s[ls[y2]]-s[ls[x1]]-s[ls[x2]];
    if(k<=k2) return query(ls[x1],ls[x2],ls[y1],ls[y2],l,mid,k);
    return query(rs[x1],rs[x2],rs[y1],rs[y2],mid+,r,k-k2);
}
int n,m,q,val[maxn],tmp[maxn],root[maxn],lastans;
int first[maxn],next[maxn*],to[maxn*],e,fa[maxn][],dep[maxn];
int pa[maxn],sz[maxn];
void AddEdge(int a,int b)
{
    to[++e]=b;next[e]=first[a];first[a]=e;
    to[++e]=a;next[e]=first[b];first[b]=e;
}
void dfs(int x,int f)
{
    dep[x]=dep[fa[x][]=f]+;update(root[x],root[f],,n,val[x]);
    sz[pa[x]]++;
    for(int i=;i<;i++) fa[x][i]=fa[fa[x][i-]][i-];
    for(int i=first[x];i;i=next[i]) if(to[i]!=f) pa[to[i]]=pa[x],dfs(to[i],x);
}
int lca(int p,int q)
{
    if(dep[p]<dep[q]) swap(p,q);
    for(int i=;i>=;i--) if((<<i)<=dep[p]-dep[q]) p=fa[p][i];
    for(int i=;i>=;i--) if(fa[p][i]!=fa[q][i]) p=fa[p][i],q=fa[q][i];
    return p==q?p:fa[p][];
}
int main()
{
    read();n=read();m=read();q=read();
    for(int i=;i<=n;i++) tmp[i]=val[i]=read();
    sort(tmp+,tmp+n+);
    for(int i=;i<=n;i++) val[i]=lower_bound(tmp+,tmp+n+,val[i])-tmp;
    for(int i=;i<m;i++) AddEdge(read(),read());
    for(int i=;i<=n;i++) if(!fa[i][]) pa[i]=i,dfs(i,);
    for(int i=;i<q;i++)
    {
        char tp=getchar();
        while(!isalpha(tp)) tp=getchar();
        if(tp=='Q')
        {
            int x=read()^lastans,y=read()^lastans,k=read()^lastans;
            int c=lca(x,y);
            print(lastans=tmp[query(root[c],root[fa[c][]],root[x],root[y],,n,k)]);
        }
        else
        {
            int u=read()^lastans,v=read()^lastans;
            int x=pa[u],y=pa[v];
            if(sz[x]>sz[y]) swap(x,y),swap(u,v);
            pa[u]=y;AddEdge(u,v);dfs(u,v);
        }
    }
    return ;
}