bzoj 3123 [Sdoi2013]森林

http://www.elijahqi.win/archives/3755
Description

Input
第一行包含一个正整数testcase，表示当前测试数据的测试点编号。保证1≤testcase≤20。
第二行包含三个整数N，M，T，分别表示节点数、初始边数、操作数。第三行包含N个非负整数表示 N个节点上的权值。
接下来 M行，每行包含两个整数x和 y，表示初始的时候，点x和点y 之间有一条无向边, 接下来 T行，每行描述一个操作，格式为“Q x y k”或者“L x y ”，其含义见题目描述部分。

Output
对于每一个第一类操作，输出一个非负整数表示答案。

Sample Input
1

8 4 8

1 1 2 2 3 3 4 4

4 7

1 8

2 4

2 1

Q 8 7 3 Q 3 5 1

Q 10 0 0

L 5 4

L 3 2 L 0 7

Q 9 2 5 Q 6 1 6

Sample Output
2

2

1

4

2
HINT

对于第一个操作 Q 8 7 3，此时 lastans=0，所以真实操作为Q 8^0 7^0 3^0，也即Q 8 7 3。点8到点7的路径上一共有5个点，其权值为4 1 1 2 4。这些权值中，第三小的为 2，输出 2，lastans变为2。对于第二个操作 Q 3 5 1 ，此时lastans=2，所以真实操作为Q 3^2 5^2 1^2 ，也即Q 1 7 3。点1到点7的路径上一共有4个点，其权值为 1 1 2 4 。这些权值中，第三小的为2，输出2，lastans变为 2。之后的操作类似。

路径第k小 主席树搞定

存在合并怎么办？ 启发式合并即可解决 每次dfs一遍 注意更新Fa的时候全部更新即可

#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int N=8e4+10;
struct node{
    int y,next;
}data[N<<1];
int fa1[N],fa[N][20],rt[N],h[N],num,dep[N],Log[N],cnt;
bool mark[N];int size[N],a[N],w[N],T,n,m,q,nn;
inline int find(int x){
    while(x!=fa1[x]) x=fa1[x]=fa1[fa1[x]];return x;
}
inline void merge1(int x,int y){
    int xx=find(x),yy=find(y);if (xx==yy) return;
    fa1[yy]=xx;size[xx]+=size[yy];
}
inline void insert1(int x,int y){
    merge1(x,y);
    data[++num].y=y;data[num].next=h[x];h[x]=num;
    data[++num].y=x;data[num].next=h[y];h[y]=num;
}
inline int lca(int x,int y){
    if (dep[x]<dep[y]) swap(x,y);int dis=dep[x]-dep[y];
    for (int i=Log[dis];~i;--i) if ((1<<i)&dis) x=fa[x][i];
    if(x==y) return x;
    for (int i=Log[dep[y]];~i;--i)
        if (fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
    return fa[x][0];
}
struct node1{
    int left,right,v;
}tree[N*400];
inline void ins(int &x,int l,int r,int p){
    tree[++cnt]=tree[x];x=cnt;++tree[x].v;
    if (l==r) return;int mid=(l+r)>>1;
    if (p<=mid) ins(tree[x].left,l,mid,p);
    else ins(tree[x].right,mid+1,r,p);
}
inline int query(int rt1,int rt2,int rt3,int rt4,int l,int r,int k){
    if (l==r) return l;int mid=(l+r)>>1;
    int l1=tree[rt1].left,l2=tree[rt2].left,l3=tree[rt3].left,l4=tree[rt4].left,
    tmp=tree[l1].v+tree[l2].v-tree[l3].v-tree[l4].v;
    if (tmp>=k) return query(l1,l2,l3,l4,l,mid,k);
    else return query(tree[rt1].right,tree[rt2].right,tree[rt3].right,tree[rt4].right,mid+1,r,k-tmp);
}
inline void dfs(int x,int f){
    dep[x]=dep[f]+1;fa[x][0]=f;mark[x]=1;
    for (int i=1;i<=Log[n];++i) fa[x][i]=fa[fa[x][i-1]][i-1];
    rt[x]=rt[f];ins(rt[x],1,nn,a[x]);
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y;if(y==fa[x][0]) continue;
        dfs(y,x);
    }
}
inline void merge(int x,int y){
    int xx=find(x),yy=find(y);
    fa1[yy]=xx;size[xx]+=size[yy];
    dfs(y,x);
}
int main(){
    freopen("bzoj3123.in","r",stdin);
    T=read();Log[0]=-1;
    n=read();m=read();q=read();
    for (int i=1;i<=n;++i) Log[i]=Log[i>>1]+1;
    for (int i=1;i<=n;++i) fa1[i]=i,size[i]=1;
    for (int i=1;i<=n;++i) a[i]=w[i]=read();int last_ans=0;
    sort(w+1,w+n+1);nn=unique(w+1,w+n+1)-w-1;memset(mark,0,sizeof(mark));
    for (int i=1;i<=n;++i) a[i]=lower_bound(w+1,w+nn+1,a[i])-w;
    for (int i=1;i<=m;++i) insert1(read(),read());
    for (int i=1;i<=n;++i) if (!mark[i]) dfs(i,0);
    while(q--){
        char op=gc();while(op!='L'&&op!='Q') op=gc();
        if (op=='Q'){
            int x=read()^last_ans,y=read()^last_ans,k=read()^last_ans;
            int f=lca(x,y),t=fa[f][0];
            last_ans=w[query(rt[x],rt[y],rt[f],rt[t],1,nn,k)];
            printf("%d\n",last_ans);
        }else{
            int x=read()^last_ans,y=read()^last_ans;
            if (size[find(x)]>size[find(y)]) merge(x,y);
            else merge(y,x);insert1(x,y);
        }
    }
    return 0;
}