[bzoj1901][整体二分]Dynamic Rankings

Description

给定一个含有n个数的序列a[1],a[2],a[3]……a[n]，程序必须回答这样的询问：对于给定的i,j,k，在a[i],a[i+1
],a[i+2]……a[j]中第k小的数是多少(1≤k≤j-i+1)，并且，你可以改变一些a[i]的值，改变后，程序还能针对改
变后的a继续回答上面的问题。

Input

第一行有两个正整数n(1≤n≤10000)，m(1≤m≤10000)。 分别表示序列的长度和指令的个数。
第二行有n个数，表示a[1],a[2]……a[n]，这些数都小于10^9。 接下来的m行描述每条指令 每行的格式是下面两种格式中的一种。
Q i j k 或者 C i t Q i j k （i,j,k是数字，1≤i≤j≤n, 1≤k≤j-i+1）
表示询问指令，询问a[i]，a[i+1]……a[j]中第k小的数。 C i t (1≤i≤n，0≤t≤10^9)表示把a[i]改变成为t
m,n≤10000

Output

对于每一次询问，你都需要输出他的答案，每一个输出占单独的一行。

Sample Input

5 3

3 2 1 4 7

Q 1 4 3

C 2 6

Q 2 5 3

Sample Output

3

6

题解

主席树模板题同时也是整体二分模板
大致说一下理解
对于一整个序列操作都需要二分的，我们可以直接二分
使用某种数据结构统计正确值在大于或者小于二分值的两个区间内的操作
重新分成两段区间再进行整体二分
其实就和cdq分治差不多.

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
struct ch{int u,v,op,gg;}q[410000],lq[410000],rq[410000];int cnt;
int n,m,ans[410000],tlen;

int s[410000];
int lowbit(int x){return x&-x;}
void change(int x,int c){while(x<=n)s[x]+=c,x+=lowbit(x);}
int findsum(int x){int ret=0;while(x>=1)ret+=s[x],x-=lowbit(x);return ret;}
//bit
void sol(int lval,int rval,int st,int ed)
{
    if(st>ed)return ;
    if(lval==rval)
    {
        for(int i=st;i<=ed;i++)if(q[i].gg)ans[q[i].gg]=lval;
        return ;
    }
    int mid=(lval+rval)>>1;
    int lt=0,rt=0;
    for(int i=st;i<=ed;i++)
    {
        if(q[i].gg)
        {
            int sum=findsum(q[i].v)-findsum(q[i].u-1);
            if(sum>=q[i].op)lq[++lt]=q[i];
            else q[i].op-=sum,rq[++rt]=q[i];
        }
        else
        {
            if(q[i].v<=mid)change(q[i].u,q[i].op),lq[++lt]=q[i];
            else rq[++rt]=q[i];
        }
    }
    for(int i=ed;i>=st;i--)if(!q[i].gg && q[i].v<=mid)change(q[i].u,-q[i].op);
    for(int i=1;i<=lt;i++)q[i+st-1]=lq[i];
    for(int i=1;i<=rt;i++)q[i+st+lt-1]=rq[i];
    sol(lval,mid,st,st+lt-1);
    sol(mid+1,rval,st+lt,ed);
}
//binary search
int a[410000];
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        cnt++;
        scanf("%d",&q[cnt].v);q[cnt].u=i;q[cnt].op=1;
        a[i]=q[cnt].v;
    }
    for(int i=1;i<=m;i++)
    {
        char ch[10];int u,v,k;
        scanf("%s%d%d",ch+1,&u,&v);
        if(ch[1]=='Q')
        {
            scanf("%d",&k);
            cnt++;q[cnt].u=u;q[cnt].v=v;q[cnt].op=k;q[cnt].gg=++tlen;
        }
        else
        {
            cnt++;q[cnt].u=u;q[cnt].v=a[u];q[cnt].op=-1;
            cnt++;q[cnt].u=u;q[cnt].v=v;q[cnt].op=1;
            a[u]=v;
        }
    }

    //init
    sol(0,1e9,1,cnt);
    for(int i=1;i<=tlen;i++)printf("%d\n",ans[i]);
    return 0;
}