[bzoj3123]森林

首先对于询问操作可以使用可持久化线段树来维护，对于连边操作对于两颗树中选取较小的树暴力练到另一个点上，点数可以用并查集然后只修改根的点数即可。

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define mid (l+r>>1)
 4 #define N 100001
 5 struct ji{
 6     int nex,to;
 7 }edge[N<<1];
 8 int E,V,n,m,t,x,y,z,ans,a[N],ff[N],sum[N],f[21][N],head[N],sz[20*N],b[N],son[2][20*N],s[N],r[N],vis[N];
 9 char s1[11];
10 void add(int x,int y){
11     edge[E].nex=head[x];
12     edge[E].to=y;
13     head[x]=E++;
14 }
15 int find(int k){
16     if (k==ff[k])return k;
17     return ff[k]=find(ff[k]);
18 }
19 int ne(int &k){
20     if (!k)k=++V;
21     return k;
22 }
23 int lca(int x,int y){
24     if (s[x]<s[y])swap(x,y);
25     for(int i=20;i>=0;i--)
26         if (s[f[i][x]]>=s[y])x=f[i][x];
27     if (x==y)return x;
28     for(int i=20;i>=0;i--)
29         if (f[i][x]!=f[i][y]){
30             x=f[i][x];
31             y=f[i][y];
32         }
33     return f[0][x];
34 }
35 void merge(int x,int y){
36     ff[find(x)]=find(y);
37     sum[find(y)]+=sum[find(x)];
38 }
39 void update(int k1,int k2,int l,int r,int x){
40     sz[k1]=sz[k2]+1;
41     if (l==r)return;
42     int p=(x>mid);
43     son[p^1][k1]=son[p^1][k2];
44     update(ne(son[p][k1]),son[p][k2],l+(mid-l+1)*p,mid+(r-mid)*p,x);
45 }
46 int query(int l,int r,int x,int a,int b,int c,int d){
47     if (l==r)return l;
48     int t=sz[son[0][a]]+sz[son[0][b]]-sz[son[0][c]]-sz[son[0][d]],p=(x>t);
49     return query(l+(mid-l+1)*p,mid+(r-mid)*p,x-p*t,son[p][a],son[p][b],son[p][c],son[p][d]);
50 }
51 void dfs(int k,int fa){
52     vis[k]=1;
53     s[k]=s[f[0][k]=fa]+1;
54     for(int i=1;i<=20;i++)f[i][k]=f[i-1][f[i-1][k]];
55     update(ne(r[k]),r[fa],1,m,a[k]);
56     for(int i=head[k];i!=-1;i=edge[i].nex)
57         if (edge[i].to!=fa)dfs(edge[i].to,k);
58 }
59 int main(){
60     scanf("%*d%d%d%d",&n,&m,&t);
61     memset(head,-1,sizeof(head));
62     for(int i=1;i<=n;i++){
63         scanf("%d",&a[i]);
64         ff[i]=i;
65         sum[i]=1;
66     }
67     for(int i=1;i<=m;i++){
68         scanf("%d%d",&x,&y);
69         add(x,y);
70         add(y,x);
71         merge(x,y);
72     }
73     memcpy(b,a,sizeof(b));
74     sort(b+1,b+n+1);
75     m=unique(b+1,b+n+1)-b-1;
76     for(int i=1;i<=n;i++)a[i]=lower_bound(b+1,b+m+1,a[i])-b;
77     for(int i=1;i<=n;i++)
78         if (!vis[i])dfs(i,0);
79     while (t--){
80         scanf("%s%d%d",s1,&x,&y);
81         x^=ans;
82         y^=ans;
83         if (s1[0]=='Q'){
84             scanf("%d",&z);
85             printf("%d\n",ans=b[query(1,m,z^ans,r[x],r[y],r[lca(x,y)],r[f[0][lca(x,y)]])]);
86         }
87         if (s1[0]=='L'){
88             if (sum[find(x)]>sum[find(y)])swap(x,y);
89             add(x,y);
90             add(y,x);
91             merge(x,y);
92             dfs(x,y);
93         }
94     }
95 }