BZOJ 2733: [HNOI2012]永无乡（线段树启发式合并）

2733: [HNOI2012]永无乡

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 3926   Solved: 2093
[ 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

HINT

线段树+并查集，感觉复杂度是q*n*log(n),但是好像启发式合并证明总复杂度只有nlog(n)。。。

#include<set>  
#include<map>     
#include<stack>            
#include<queue>            
#include<vector>    
#include<string> 
#include<time.h>
#include<math.h>            
#include<stdio.h>            
#include<iostream>            
#include<string.h>            
#include<stdlib.h>    
#include<algorithm>   
#include<functional>    
using namespace std;            
#define ll __int64        
#define inf  1000000000       
#define mod 1000000007             
#define maxn  100100
#define lowbit(x) (x&-x)            
#define eps 1e-9
int ls[maxn*18],rs[maxn*18],sum[maxn*18];
int rt[maxn],a[maxn],n,m,k,siz,q,p[maxn],id[maxn];
int find(int x)
{
	if(p[x]==x)
		return x;
	return p[x]=find(p[x]);
}
void insert(int &id,int l,int r,int val)
{
	if(!id)
		id=++siz;
	if(l==r)
	{
		sum[id]=1;
		return;
	}
	int mid=(l+r)/2;
	if(val<=mid)
		insert(ls[id],l,mid,val);
	else
		insert(rs[id],mid+1,r,val);
	sum[id]=sum[ls[id]]+sum[rs[id]];
}
int merge(int x,int y)
{
	if(!x) return y;
	if(!y) return x;
	ls[x]=merge(ls[x],ls[y]);
	rs[x]=merge(rs[x],rs[y]);
	sum[x]=sum[ls[x]]+sum[rs[x]];
	return x;
}
int query(int id,int l,int r,int val)
{
	if(l==r)
		return l;
	int mid=(l+r)/2;
	if(sum[ls[id]]>=val)
		return query(ls[id],l,mid,val);
	else
		return query(rs[id],mid+1,r,val-sum[ls[id]]);
}
int main(void) 
{
	char c;
	int i,j,x,y,t1,t2;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++)
		scanf("%d",&a[i]),p[i]=i;
	for(i=1;i<=m;i++)
	{
		scanf("%d%d",&x,&y);
		t1=find(x);t2=find(y);
		p[t1]=t2;
	}
	for(i=1;i<=n;i++)
	{
		  insert(rt[find(i)],1,n,a[i]);
		  id[a[i]]=i;
	}
	scanf("%d",&q);
	while(q--)
	{
		scanf(" %c%d%d",&c,&x,&y);
		if(c=='B')
		{
			t1=find(x);
			t2=find(y);
			if(t1!=t2)
			{
				p[t1]=t2;
				rt[t2]=merge(rt[t1],rt[t2]);
			}
		}
		else
		{
			t1=find(x);
			if(sum[rt[t1]]<y)
			{
				printf("-1\n");
				continue;
			}
			int t=query(rt[t1],1,n,y);
			printf("%d\n",id[t]);
		}
	}
	return 0;
}