bzoj 2733(splay启发式合并)

2733: [HNOI2012]永无乡

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 2566   Solved: 1353
[ Submit][ Status][ Discuss]

Description

永无乡包含 n 座岛，编号从 1 到 n，每座岛都有自己的独一无二的重要度，按照重要度可 以将这 n 座岛排名，名次用 1 到 n 来表示。某些岛之间由巨大的桥连接，通过桥可以从一个岛 到达另一个岛。如果从岛 a 出发经过若干座（含 0 座）桥可以到达岛 b，则称岛 a 和岛 b 是连 通的。现在有两种操作：B x y 表示在岛 x 与岛 y 之间修建一座新桥。Q x k 表示询问当前与岛 x连通的所有岛中第 k 重要的是哪座岛，即所有与岛 x 连通的岛中重要度排名第 k 小的岛是哪 座，请你输出那个岛的编号。 
 

Input

输入文件第一行是用空格隔开的两个正整数 n 和 m，分别 表示岛的个数以及一开始存在的桥数。接下来的一行是用空格隔开的 n 个数，依次描述从岛 1 到岛 n 的重要度排名。随后的 m 行每行是用空格隔开的两个正整数 ai 和 bi，表示一开始就存 在一座连接岛 ai 和岛 bi 的桥。后面剩下的部分描述操作，该部分的第一行是一个正整数 q， 表示一共有 q 个操作，接下来的 q 行依次描述每个操作，操作的格式如上所述，以大写字母 Q 或B 开始，后面跟两个不超过 n 的正整数，字母与数字以及两个数字之间用空格隔开。 对于 20%的数据 n≤1000,q≤1000 
 
对于 100%的数据 n≤100000,m≤n，q≤300000 
 

Output

对于每个 Q x k 操作都要依次输出一行，其中包含一个整数，表 示所询问岛屿的编号。如果该岛屿不存在，则输出-1。 
 

Sample Input

5 1
4 3 2 5 1
1 2
7
Q 3 2
Q 2 1
B 2 3
B 1 5
Q 2 1
Q 2 4
Q 2 3

Sample Output

-1
2
5
1
2

解题思路：splay+启发式合并。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
int n,m;
int fa[110000],f[110000],size[110000],l[110000],r[110000];
int q[110000],zhi[110000];


inline int read()
{
char y; int x=0,f=1; y=getchar();
while (y<'0' || y>'9') {if (y=='-')f=-1; y=getchar();}
while (y>='0' && y<='9') {x=x*10+int(y)-48; y=getchar();}
return x*f;
}


int find(int now)
 {
  if (f[now]!=now) f[now]=find(f[now]);
  return f[now];
 }


void pushdown(int now)
 {
  size[now]=size[l[now]]+size[r[now]]+1;
 }


void rotate(int x)
 {
  int y=fa[x]; int z=fa[y];
if (z!=0)
{
if (l[z]==y) l[z]=x;else r[z]=x;
 } 
fa[x]=z; fa[y]=x;
if (l[y]==x)
{
fa[r[x]]=y; l[y]=r[x]; r[x]=y; 
}else
 {
  fa[l[x]]=y; r[y]=l[x]; l[x]=y;
 }
pushdown(y); pushdown(x);
 }


void splay(int x)
 {
  while (fa[x])
  {
   int y=fa[x]; int z=fa[y];
   if (z!=0)
    if (l[z]==y ^ l[y]==x) rotate(y);else rotate(x);
   rotate(x);
}
 }


void insert(int z,int &now,int from)
 {
  if (now==0) {now=z;fa[z]=from;return;}
    ++size[now]; 
    if (zhi[z]<zhi[now]) insert(z,l[now],now);
     else insert(z,r[now],now);
 }


void work(int x,int y)
 {
  splay(x); splay(y);
if (size[x]>size[y]) swap(x,y);
int tail=1,head=0; q[tail]=x;
while (head<tail)
{
++head;
if (l[q[head]])
{
++tail; q[tail]=l[q[head]];
}
if (r[q[head]])
{
++tail; q[tail]=r[q[head]];
}
 }
for (int i=1;i<=tail;++i)
{
insert(q[i],y,0); 
splay(q[i]);
}
 }


int find(int sum,int now)
 {
  if (size[l[now]]+1==sum) return now;
  if (size[l[now]]+1<sum) return find(sum-size[l[now]]-1,r[now]);
  else return find(sum,l[now]);
 }


int main()
{
n=read(); m=read();
for (int i=1;i<=n;++i)
{
size[i]=1; zhi[i]=read(); f[i]=i;
}
for (int i=1;i<=m;++i)
{
int x=read(); int y=read();
int u1=find(x); int u2=find(y);
if (u1==u2) continue;
work(x,y);
}
int q=read();
for (int i=1;i<=q;++i)
{
  char c[100];
  scanf("%s",c);
  if (c[0]=='Q')
   {
     int x,k; x=read(); k=read(); splay(x);
     if (size[x]<k) 
      {
      printf("-1\n");
      continue;
  }
     int ans=find(k,x); printf("%d\n",ans);
}else
{
int x,y; x=read(); y=read();
int u1=find(x); int u2=find(y);
if (u1==u2) continue;
work(x,y); 
}
}
}