【ZOJ 2112 Dynamic Rankings】动态区间第k大：树状数组+主席树

思路纯属个人YY，如有错误请各位dalao不吝赐教。Orz
萌新第一次接触这么阔怕的东西，瑟瑟发抖啊嘤嘤嘤。。。
ZOJ 2112 Dynamic Rankings
题目大意就是，求区间第k大，支持单点修改。

首先，区间第k大，那肯定是要用传说中的主席树了。
（刚刚接触主席树的萌新持续瑟瑟发抖）

---

首先想一下，如果是整体第k大，应该怎么算呢？sort
可以开一棵权值线段树，对于每个叶子 i，表示在整个区间里，i 这个数出现了多少次，然后向上维护和~
这样查询的时候，二分就阔以了~

---

那区间和整体该怎么转换呢？当然是搞一堆权值线段树，分别储存各个前缀区间的整体第k大，然后加加减减就好啦~
这样做的话，很多节点都是重复的（毕竟每次都只加进去1个数字），所以可以使用各种指针乱指，使用之前已经建好的节点，避免浪费啦~
大概就是这个亚子~

[随便从度娘那里偷的图]

---

现在我们阔以查询区间第k大了，但是如果加上修改的话，该怎么办呢？~
把某一个数修改了（我们假设把a改成b），就需要修改主席树里以$min（a,b）$为区间右端点开始往后的所有树，$nlogn$的巨大复杂度显然无法接受。

但是觉不觉得目前为止我们的各种操作都很眼熟呢？前缀和相减求区间和，要提升修改的效率。。好像有个$logn$查询，$logn$修改的东西……树状数组！
[再放一张偷的图，from star-city dalao]
维护的时候使用树状数组来维护，每个树状数组树状之前的节点（什么鬼）都看作一个前缀区间和。实际上树状数组就长这个亚子：
// 偷图from FlashHu dalao

---

那现在该怎么查询呢？
插入之前是两棵树一起向下跳，求差来获得$[l,r]$，二分获得第k大。树状数组维护之后的话当然是$logn$棵树和$logn$棵树一起跳啦~

---

代码是抄的kuangbin dalao的板子
自己乱加了一些注释，希望有一丢丢用处。。
要记得空间开40倍啊！//可怕

OTZ

// 单点修改+查询区间第k大
//树状数组套主席树，ZOJ2112 
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 60010;
const int M = 2500010;//主席树本体的节点数，O(nlogn + mlogn),要开40倍 
int n,q,m,tot;
int a[N],//原始数组 
t[N];//所有权值 
int T[N],//n个权值的权值线段树的树根 
lson[M],
rson[M],
c[M];//主席树本体的各个节点 
int S[N];//树状数组的数组 
struct Query{
	int kind;
	int l,r,k;
}query[10010];

void Init_hash(int k){
	sort(t,t+k);
	m = unique(t,t+k) -t;
} 

int hash(int x){
	return lower_bound(t,t+m,x) -t;
}

int build(int l,int r){
	int root = tot ++;
	c[root] = 0;
	if(l != r){
		int mid = (l+r)>>1;
		lson[root] = build(l,mid);
		rson[root] = build(mid+1,r);
	}
	return root;
}

int Insert(int root,int pos,int val){
	int newroot = tot ++,tmp = newroot;
	int l = 0,r = m-1;
	c[newroot] = c[root] + val;
	while(l < r){
		int mid = (l + r) >> 1;
		if(pos <= mid){
			lson[newroot] = tot ++;
			rson[newroot] = rson[root];
			newroot = lson[newroot];
			root = lson[root];
			r = mid;
		}
		else{
			rson[newroot] = tot ++;
			lson[newroot] = lson[root];
			newroot = rson[newroot];
			root = rson[root];
			l = mid +1; 
		}
		c[newroot] = c[root] + val;
	}
	return tmp;
}

int lowbit(int x){
	return x&(-x);
}

int use[N];
void add(int x,int pos,int val){
	while(x <= n){
		S[x] = Insert(S[x],pos,val);
		x += lowbit(x);
	}
}

int sum(int x){
	int ret = 0;
	while(x > 0){
		ret += c[lson[use[x]]];
		x -= lowbit(x);
	}
	return ret;
}

int Query(int left,int right,int k){
	int left_root = T[left-1];
	int right_root = T[right];
	int l = 0,r = m-1;
	for(int i = left-1; i; i -= lowbit(i)) use[i] = S[i];
	for(int i = right; i; i -= lowbit(i)) use[i] = S[i];
	while(l < r){ //二分找第k大 
		int mid = (l+r) >> 1;
		int tmp = sum(right) - sum(left-1) + c[lson[right_root]] - c[lson[left_root]];
		if(tmp >= k){
			r = mid;
			for(int i=  left-1; i; i -= lowbit(i))
				use[i] = lson[use[i]];
			for(int i = right; i; i-= lowbit(i))
				use[i] = lson[use[i]];
			left_root = lson[left_root];
			right_root = lson[right_root];
		}
		else{
			l = mid+1;
			k -= tmp;
			for(int i=  left-1; i; i -= lowbit(i))
				use[i] = rson[use[i]];
			for(int i = right; i; i-= lowbit(i))
				use[i] = rson[use[i]];
			left_root = rson[left_root];
			right_root = rson[right_root];
		}
	}
	return l;
}

void Modify(int x,int p,int d){
	while(x <= n){
		S[x] = Insert(S[x],p,d);
		x += lowbit(x);
	}
}

int main(){
	int Tcase;
	scanf("%d",&Tcase);
	while(Tcase --){
		scanf("%d%d",&n,&q); //n个点，q次询问 
		tot = 0;
		m = 0;
		for(int i = 1;i <= n;i ++){
			scanf("%d",&a[i]);
			t[m ++] = a[i];
		}
		char op[10];
		for(int i = 0;i < q;i ++){
			scanf("%s",op);
			if(op[0] == 'Q'){
				query[i].kind = 0;
				scanf("%d%d%d",&query[i].l,&query[i].r,&query[i].k);
			}
			else{
				query[i].kind = 1;
				scanf("%d%d",&query[i].l,&query[i].r);
				t[m ++] = query[i].r;//先保存下来询问 主要就是因为这个！ 
			}
		}
		Init_hash(m);//离散化，方便建立权值线段树 
		T[0] = build(0,m-1); //先建一棵空树
		for(int i = 1; i <= n;i ++) //主席树~统计1~i区间每个值出现的次数 
			T[i] = Insert(T[i-1],hash(a[i]),1); 
		for(int i = 1;i <= n;i ++){
			S[i] = T[0];
		}
		for(int i = 0;i < q;i ++){
			if(query[i].kind == 0)
				printf("%d\n",t[Query(query[i].l,query[i].r,query[i].k)]); //l~r第k大 
			else{
				Modify(query[i].l,hash(a[query[i].l]),-1); //修改l处元素，改为r 
				Modify(query[i].l,hash(query[i].r),1);
				a[query[i].l] = query[i].r;
			}
		}
	}
	return 0;
}