bzoj 2733: [HNOI2012]永無鄉 (線段樹合併)[省選計劃系列]

2733: [HNOI2012]永無鄉

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

和bzoj3545差不多，3545我寫的splay的啓發式：

bzoj 3545: [ONTAK2010]Peaks (splay啓發式合併)

這題我用的線段樹合併

代碼：

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<algorithm>
#define N 300010
#define M 2000010
using namespace std;
inline int read()
{
    int x=0,f=1;char ch;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,q;
int root[N],s[M],ls[M],rs[M],sz,id[N];
int f[N],v[N];
int find(int x){return x==f[x]?x:f[x]=find(f[x]);}
void insert(int &u,int l,int r,int val)
{
    if(!u)u=++sz;s[u]++;
    if(l==r)return;int mid=(l+r)>>1;
    if(val<=mid)insert(ls[u],l,mid,val);
    else insert(rs[u],mid+1,r,val);
}
int ask(int u,int l,int r,int x,int y)
{
    if(!u)return 0;int mid=(l+r)>>1;
    if(x<=l&&y>=r)return s[u];
    if(y<=mid)return ask(ls[u],l,mid,x,y);
    else if(x>mid)return ask(rs[u],mid+1,r,x,y);
    else return ask(ls[u],l,mid,x,mid)+ask(rs[u],mid+1,r,mid+1,y);
}
int solve(int a,int rk)
{
    int ta=find(a);
    int rt=root[ta];
    if(s[rt]<rk)return 0;
    int l=1,r=n,ret;
    while(l<=r)
    {
		int mid=(l+r)>>1;
		if(ask(rt,1,n,1,mid)>=rk)
		    r=mid-1,ret=mid;
		else l=mid+1;
    }return ret;
}
int merge(int x,int y)
{
    if(!x&&!y)return 0;
    if(!x)return y;
    if(!y)return x;
    ls[x]=merge(ls[x],ls[y]);
    rs[x]=merge(rs[x],rs[y]);
    s[x]=s[ls[x]]+s[rs[x]];
    return x;
}

int main()
{
    n=read(),m=read();
    for(int i=1;i<=n;i++)
		v[i]=read(),f[i]=i,id[v[i]]=i;
    for(int i=1;i<=m;i++)
    {
		int x=read(),y=read();
		int fx=find(x),fy=find(y);
		if(fx==fy)continue;f[fx]=fy;
    }
    for(int i=1;i<=n;i++)
		insert(root[find(i)],1,n,v[i]);
    q=read();char ch[10];id[0]=-1;
    for(int i=1;i<=q;i++)
    {
		scanf("%s",ch);
		int x=read(),y=read();
		if(ch[0]=='Q')
		    printf("%d\n",id[solve(x,y)]);
		else if(ch[0]=='B')
		{
		    int fx=find(x),fy=find(y);
		    if(fx==fy)continue;
		    f[fx]=fy;
		    root[fy]=merge(root[fx],root[fy]);
		}
    }return 0;
}