BZOJ 3123: [Sdoi2013]森林【可持久化值域线段树+启发式合并

裸的可持久化值域线段树+启发式（贪心又暴力地）合并

……合并的时候忘记修改深度了……不知道为什么现在稍微长一点的代码就不想调了【好颓啊QwQ

没什么好说的，一切都是套路……

求第k大/第k小很明显满足加减的性质于是

要求链上的值 其实只用记录每个节点到根的数据……然后进行加加减减【根本不用链剖【神犇们说倍增LCA太丑于是都去学LCT求LCA了我就…………跪着

#include<bits/stdc++.h>
#define MAXN 80057
#define MAX2 20000057
using namespace std;	int n,m,t;
int MX;

int cnt[MAX2],cnt_node=0;
int leftson[MAX2],rightson[MAX2];

void insert(int pre,int &now,int l,int r,int key){
cnt[++cnt_node]=cnt[now];
now=cnt_node;
++cnt[now];
if(l==r)	return ;
int mid=(l+r)>>1;
leftson[now]=leftson[pre] , rightson[now]=rightson[pre];
if(key<=mid)	insert(leftson[pre],leftson[now],l,mid,key);
else	insert(rightson[pre],rightson[now],mid+1,r,key);
}

int inqry(int r1,int r2,int r3,int r4,		int l,int r,int k){
if(l==r)	return l;
int mid=(l+r)>>1;
int dt_l=cnt[leftson[r1]]+cnt[leftson[r2]]
-cnt[leftson[r3]]-cnt[leftson[r4]];
if(k<=dt_l)	return inqry(leftson[r1],leftson[r2],leftson[r3],leftson[r4],l,mid,k);
else	return inqry(rightson[r1],rightson[r2],rightson[r3],rightson[r4],mid+1,r,k-dt_l);
}

//=================================================
struct t1{
int to,nxt;
}edge[MAXN<<1];	int cnt_edge=0;
int fst[MAXN];
void addedge(int x,int y){
edge[++cnt_edge].to=y;
edge[cnt_edge].nxt=fst[x];
fst[x]=cnt_edge;
}
int rec[MAXN];

int root[MAXN];

int anc[MAXN][17];
void jump(int now){
for(int i=1;i<17;++i)
anc[now][i]=anc[ anc[now][i-1] ][i-1];
}

int root_num[MAXN],siz[MAXN],dpt[MAXN];
void dfs(int now){
siz[now]=1;
jump(now);
insert(root[ anc[now][0] ],root[now],1,MX,rec[now]);
for(int tmp=fst[now];tmp;tmp=edge[tmp].nxt){
if(edge[tmp].to==anc[now][0])	continue;
int ch=edge[tmp].to;
anc[ch][0]=now;
root_num[ch]=root_num[now];
dpt[ch]=dpt[now]+1;

dfs(ch);
siz[now]+=siz[ch];
}
}

int lca(int x,int y){
if(dpt[x]<dpt[y])	swap(x,y);
int dlt=dpt[x]-dpt[y];
for(int i=16;~i;--i)
if(dlt&(1<<i))	x=anc[x][i];
if(x==y)	return x;
for(int i=16;~i;--i){
if(anc[x][i]!=anc[y][i])
x=anc[x][i], y=anc[y][i];
}
if(x==y)	return x;
return anc[x][0];
}

//============================================================
int work_Q(int x,int y,int k){
int LCA=lca(x,y);
int F_lca=anc[LCA][0];
return inqry(root[x],root[y],root[LCA],root[ anc[LCA][0] ],1,MX,k);
}

void work_L(int x,int y){
addedge(x,y),addedge(y,x);
anc[y][0]=x;
siz[root_num[x]]+=siz[root_num[y]];
root_num[y]=root_num[x];
dpt[y]=dpt[x]+1;
dfs(y);
}

char opt[5];
int read_x,read_y,read_k;
int lastans=0;
int main(){
scanf("%*d%d%d%d",&n,&m,&t);
for(int i=1;i<=n;++i)	scanf("%d",rec+i),MX=max(MX,rec[i]);
for(int i=1;i<=m;++i)
scanf("%d%d",&read_x,&read_y),
addedge(read_x,read_y), addedge(read_y,read_x);

for(int i=1;i<=n;++i)
if(!root_num[i])
root_num[i]=i,dpt[i]=1, dfs(i);

while(t--){
scanf("%s",opt);
if(opt[0]=='Q'){
scanf("%d%d%d",&read_x,&read_y,&read_k);
read_x^=lastans,read_y^=lastans,read_k^=lastans;
printf("%d\n",lastans=work_Q(read_x,read_y,read_k));
}
else{
scanf("%d%d",&read_x,&read_y);
read_x^=lastans,read_y^=lastans;
if(siz[root_num[read_x]]<siz[root_num[read_y]])	swap(read_x,read_y);
work_L(read_x,read_y);
}
}

return 0;
}