ZOJ2112 Dynamic Rankings(整体二分)
今天学习了一个奇技淫巧－－整体二分。关于整体二分的一些理论性的东西，可以参见XRH的《浅谈数据结构题的几个非经典解法》。然后下面是一些个人的心得体会吧，写下来希望加深一下自己的理解，或者如果有人看了或许也有些帮助。
ZOJ2112是一道典型的带修改的区间第k大的问题，有一些树套树等的数据结构可以在线处理这样的问题。但是当题目并不要求在线处理的时候，其实我们可以选择一下整体二分的思想。
个人对整体二分的理解是这样子的，首先对于修改，即把a[xi]=yi（1<=yi<=C）的时候，我们可以把修改的操作划分成两部分，一部分是yi<=mid,另一部分是yi>mid。首先我们忽略掉yi>mid的操作，将yi<=mid以及询问操作按照输入的顺序执行一遍，这样我们就可以知道<=mid的操作对所有询问的贡献，接下来我们就要根据贡献对修改操作划分成两部分，一部分是答案在<=mid里面的，另外一部分是答案在>mid里面的，对于<=mid的，显然我们可以递归求解（因为>mid的操作对那些<=mid是没有影响的），而对于>mid的来说，<=mid的操作是有影响的，但是<=mid所造成的影响已经算过一遍了，所以对于>mid的来说，其实也是独立的。所以如果我们把答案二分的范围的大小设为C，操作的总数设为n，那么划分到<=mid的操作是n1,>mid的操作是n2时（n1+n2=n）
有 T(n,C)=T(n1,C/2)+T(n2,C/2)+O(nlogC) T(n)=nlogC^2...（也不知道复杂度对不对，可以画一个图来YY一下）
好吧，其实我也不知道上面说了什么，感觉真的蛮难用一两段话将整个思路说出来的。其实总的抽象的来说，就是一开始维护的可能答案的区间是[l,r]，
然后我们将所有<=mid的操作作一遍，然后将答案落在[l,mid]的操作（询问和修改）划到一边，将答案落在[mid+1,r]的操作（询问和修改）划到另一边。
#pragma warning(disable:4996)
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
#include <string>
#include <vector>
#include <queue>
using namespace std;
#define maxn 300000
struct Query
{
int x, y; // qt==3 [x,y] qt==1||2 a[x]=y
int qt; // query type 1 for increase 2 for decrease 3 for query
int cur; // the current contribution
int k; // the query is the k-th smalleset
int index;
}q[maxn];
int qtop;
int a[maxn];
int b[maxn];
int bsize;
int n, m;
int tot;
int ans[maxn];
int ansid;
int tmp[maxn];
Query q1[maxn], q2[maxn];
int bit[maxn];
void add(int x, int v)
{
while (x <= n){
bit[x] += v;
x += x&(-x);
}
}
int query(int x)
{
int ret = 0;
while (x > 0){
ret += bit[x];
x -= x&(-x);
}
return ret;
}
void solve(int head, int tail, int l, int r)
{
if (head > tail) return;
if (l == r){
for (int i = head; i <= tail; ++i){
if (q[i].qt == 3) ans[q[i].index] = l;
}
return;
}
int mid = (l + r) >> 1;
// 将所有<=mid的操作作一遍，tmp[i]存的是[head,tail]里<=mid的操作对询问的贡献
for (int i = head; i <= tail; ++i){
if (q[i].qt == 1 && q[i].y <= mid){
add(q[i].x, 1);
}
else if (q[i].qt == 2 && q[i].y <= mid){
add(q[i].x, -1);
}
else{
tmp[i] = query(q[i].y) - query(q[i].x - 1);
}
}
// 将操作撤销一下
for (int i = head; i <= tail; ++i){
if (q[i].qt == 1 && q[i].y <= mid){
add(q[i].x, -1);
}
else if (q[i].qt == 2 && q[i].y <= mid){
add(q[i].x, 1);
}
}
// 将操作划分成两部分
int l1=0, l2 = 0;
for (int i = head; i <= tail; ++i){
if (q[i].qt == 3){
// 如果前面的数加上当前的数>=q[i].k，说明该询问的可行区间在[l,mid]，往左划分
if (q[i].cur + tmp[i] >= q[i].k){
q1[++l1] = q[i];
}
else{
// 否则往右划分，并记下贡献
q[i].cur += tmp[i];
q2[++l2] = q[i];
}
}
else{
if (q[i].y <= mid) q1[++l1] = q[i];
else q2[++l2] = q[i];
}
}
for (int i = 1; i <= l1; ++i) {
q[head + i - 1] = q1[i];
}
for (int i = 1; i <= l2; ++i){
q[head + l1 + i - 1] = q2[i];
}
solve(head, head + l1 - 1, l, mid);
solve(head + l1, tail, mid + 1, r);
}
int main()
{
int T; cin >> T;
while (T--)
{
scanf("%d%d", &n, &m);
qtop = 0;tot = 0;
for (int i = 1; i <= n; ++i){
scanf("%d", a + i);
b[tot++] = a[i];
q[++qtop].qt = 1;
q[qtop].x = i;
q[qtop].y = a[i];
}
char cmd[4];
int xi, yi,ki;
ansid = 0;
for (int i = 0; i < m; ++i){
scanf("%s", cmd);
if (cmd[0] == 'Q'){
scanf("%d%d%d", &xi, &yi, &ki);
q[++qtop].x = xi; q[qtop].y = yi; q[qtop].k = ki;
q[qtop].qt = 3; q[qtop].cur = 0; q[qtop].index = ++ansid;
}
else{
scanf("%d%d", &xi, &yi);
q[++qtop].x = xi; q[qtop].y = a[xi]; q[qtop].qt = 2;
q[++qtop].x = xi; q[qtop].y = yi; q[qtop].qt = 1;
a[xi] = yi;
b[tot++] = yi;
}
}
sort(b, b + tot);
bsize = unique(b, b + tot) - b;
for (int i = 1; i <= qtop; ++i){
if (q[i].qt == 1 || q[i].qt == 2){
q[i].y = lower_bound(b, b + bsize, q[i].y) - b + 1;
}
}
solve(1, qtop, 1, bsize);
for (int i = 1; i <= ansid; ++i){
printf("%d\n", b[ans[i]-1]);
}
}
//system("pause");
return 0;
}