[BZOJ 2733][HNOI2012]永无乡

并没有什么, Splay的启发式合并


[HNOI2012]永无乡

【题目描述】


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


输入文件第一行是用空格隔开的两个正整数 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 
 
【输出格式】


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


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 
  
【样例输出】


-1
  2
  5
  1
  2


#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define maxn 100010
using namespace std;

//*---------------------------------------------------------------*//
int c[maxn][2], size[maxn], fa[maxn], val[maxn];

void pushup(int x){
	size[x] = size[c[x][0]] + size[c[x][1]] + 1;
}

void rotate(int p, int x){
	int mark = p == c[x][1], y = c[p][mark ^ 1], z = fa[x];
	if(c[z][0] == x)c[z][0] = p;
	if(c[z][1] == x)c[z][1] = p;
	if(y)fa[y] = x;
	fa[p] = z;c[p][mark ^ 1] = x;
	fa[x] = p;c[x][mark] = y;
	pushup(x);
}

void splay(int p){
	while(fa[p]){
		int x = fa[p], y = fa[x];
		if(y == 0)rotate(p, x);
		else if(p == c[x][0] ^ x == c[y][0])
		    rotate(p, x), rotate(p, y);
		else rotate(x, y), rotate(p, x);
	}
	pushup(p);
}

void insert(int p, int x){
	//splay(x);
	int now = x, ret = now;
	while(now){
		ret = now;
		if(val[p] < val[now])
		    now = c[now][0];
		else now = c[now][1];
	}
	c[ret][val[p] > val[ret]] = p;
	fa[p] = ret;
	c[p][0] = c[p][1] = 0;
	size[p] = 1;
	splay(p);
}

int Qry(int p, int k){
	splay(p);
	if(size[p] < k)return -1;
	int now = p;
	while(true){
		if(size[c[now][0]] + 1 == k)
			return now;
		if(size[c[now][0]] >= k)
		    now = c[now][0];
		else{
			k -= size[c[now][0]] + 1;
			now = c[now][1];
		}
	}
}

int que[maxn];

void Union(int x, int y){
	splay(x), splay(y);
	if(size[x] > size[y])
	    swap(x, y);
	int head = 0, tail = 0;
	que[tail ++] = x;
	while(head != tail){
		int u = que[head ++];
		if(c[u][0])que[tail ++] = c[u][0];
		if(c[u][1])que[tail ++] = c[u][1];
  		insert(u, y);
	}
}
//*---------------------------------------------------------------*//

int n, m;

int main(){
    freopen("bzoj_2733.in","r",stdin);
	freopen("bzoj_2733.out","w",stdout);
	scanf("%d%d", &n, &m);
	for(int i = 1; i <= n; i ++)
	    scanf("%d", &val[i]);
	int u, v;
	for(int i = 1; i <= n; i ++)
	    size[i] = 1;

	for(int i = 1; i <= m; i ++){
		scanf("%d%d", &u, &v);
		if(u && v)
			Union(u, v);
	}

	scanf("%d", &m);
	while(m --){
		char QAQ = getchar();
		for(; QAQ < '!'; QAQ = getchar());
		scanf("%d%d", &u, &v);
		if(QAQ == 'Q')
			printf("%d\n", Qry(u, v));
		else
			Union(u, v);
	}
	return 0;
}



  

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值