bzoj 3123（主席树+启发式合并）

传送门
题解：对点到根维护权值线段树，每次连边就进行启发式合并（用size启发），询问就按区间第k大的方式询问，此处不再赘述。

吐槽三点：
1.说好的多组数据其实只有一组，醉了。。。
2.1e9的数据按理说是需要离散化的，但是离散化了反而要RE（可能是个人原因）
3.这个代码片的int为啥是畸形的。。。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXN=1e5+2;
int n,m,q,sb,maxx=0,minx=1e9+7;
int a[MAXN];//int b[MAXN];
int fa[MAXN],siz[MAXN];
int root[MAXN],sum[MAXN*260],lc[MAXN*260],rc[MAXN*260],tim=0;
int head[MAXN],edge=0,dep[MAXN],anc[MAXN][25];
struct EDGE {
    int v,nxt;
}e[MAXN<<1];
inline int read() {
    int x=0,f=1;char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
    while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*f;
}
inline void adde(int u,int v) {
    e[edge].nxt=head[u],e[edge].v=v,head[u]=edge++;
    e[edge].nxt=head[v],e[edge].v=u,head[v]=edge++;
}
int find(int x) {return fa[x]==x?x:fa[x]=find(fa[x]);}
void modify(int &rt,int pre,int l,int r,int x) {
    rt=++tim;
    sum[rt]=sum[pre]+1;
    if (l==r) return ;
    lc[rt]=lc[pre],rc[rt]=rc[pre];
    int mid=(l+r)>>1;
    if (x<=mid) modify(lc[rt],lc[pre],l,mid,x);
    else modify(rc[rt],rc[pre],mid+1,r,x);
}
int query(int a,int b,int c,int d,int l,int r,int k) {
    if (l==r) return l;
    int temp=sum[lc[a]]+sum[lc[b]]-sum[lc[c]]-sum[lc[d]];
    int mid=(l+r)>>1;
    if (k<=temp) return query(lc[a],lc[b],lc[c],lc[d],l,mid,k);
    else return query(rc[a],rc[b],rc[c],rc[d],mid+1,r,k-temp);
}
void dfs(int p,int fat) {
    dep[p]=dep[fat]+1;
    anc[p][0]=fat;
    for (int i=1;i<=18;++i)
        anc[p][i]=anc[anc[p][i-1]][i-1];
    modify(root[p],root[fat],minx,maxx,a[p]);
    for (int i=head[p];~i;i=e[i].nxt) {
        int v=e[i].v;
        if (v!=fat) dfs(v,p);
    }
}
inline int lca(int a,int b) {
    if (dep[a]<dep[b]) a^=b^=a^=b;
    for (int i=18;i>=0;--i)
        if (anc[a][i]&&dep[a]-(1<<i)>=dep[b]) a=anc[a][i];
    if (a==b) return a;
    for (int i=18;i>=0;--i)
        if (anc[a][i]&&anc[b][i]&&anc[a][i]^anc[b][i])
            a=anc[a][i],b=anc[b][i];
    return anc[a][0];
}
int main() {
//  freopen("bzoj 3123.in","r",stdin);
    memset(sum,0,sizeof(sum));
    memset(dep,0,sizeof(dep));
    memset(head,-1,sizeof(head));
    memset(root,0,sizeof(root));
    sb=read();//in fact, only one case...
    n=read(),m=read(),q=read();
    for (register int i=1;i<=n;++i)
        a[i]=read(),maxx=max(maxx,a[i]),minx=min(minx,a[i]),siz[i]=1,fa[i]=i;
//  sort(b+1,b+n+1);
//  int nn=unique(b+1,b+n+1)-b-1;
//  for (register int i=1;i<=n;++i) a[i]=lower_bound(b+1,b+nn+1,a[i])-b;
    for (register int i=0;i<m;++i) {
        int u=read(),v=read();
        adde(u,v);
        u=find(u),v=find(v);
        if (u^v) {
            if (siz[u]<siz[v]) u^=v^=u^=v;
            fa[v]=u,siz[u]+=siz[v];
        }
    }
    for (register int i=1;i<=n;++i)
        if (!dep[i]) dfs(i,0);
    int ans=0;
    for (register int i=0;i<q;++i) {
        char ss;
        while (!isalpha(ss=getchar()));
        if (ss^'L') {
            int x=read()^ans,y=read()^ans,k=read()^ans;
            int ff=lca(x,y);
            ans=query(root[x],root[y],root[ff],root[anc[ff][0]],minx,maxx,k);
            printf("%d\n",ans);
        }
        else {
            int x=read()^ans,y=read()^ans;
            int xx=find(x),yy=find(y);
            if (siz[xx]<siz[yy]) xx^=yy^=xx^=yy,x^=y^=x^=y;
            fa[yy]=xx,siz[xx]+=siz[yy];
            adde(x,y);
            dfs(y,x);
        }
    }
    return 0;
}