[BZOJ]3065: 带插入区间K小值 块状链表

Description
从前有n只跳蚤排成一行做早操，每只跳蚤都有自己的一个弹跳力a[i]。跳蚤国王看着这些跳蚤国欣欣向荣的情景，感到非常高兴。这时跳蚤国王决定理性愉悦一下，查询区间k小值。他每次向它的随从伏特提出这样的问题: 从左往右第x个到第y个跳蚤中，a[i]第k小的值是多少。
这可难不倒伏特，他在脑袋里使用函数式线段树前缀和的方法水掉了跳蚤国王的询问。
这时伏特发现有些跳蚤跳久了弹跳力会有变化，有的会增大，有的会减少。
这可难不倒伏特，他在脑袋里使用树状数组套线段树的方法水掉了跳蚤国王的询问。（orz 主席树）
这时伏特发现有些迟到的跳蚤会插入到这一行的某个位置上，他感到非常生气，因为……他不会做了。
请你帮一帮伏特吧。
快捷版题意：带插入、修改的区间k小值在线查询。

题解：

分块。维护每块的原数组与排序后数组，当某块的size刚好达到两个块的大小时，分裂成两块，其它都是普通分块的套路了。如果要对拍，记得在对拍完之后加上强制在线，别问我是怎么知道的……

代码：

#include<bits/stdc++.h>using namespace std;#define LL long long#define pa pair<int,int>const int block_size=600;const int Maxn=70010;const int inf=2147483647;int read(){    int x=0,f=1;char ch=getchar();    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();    return x*f;}int n,num,next[120],b[120][1210],a[120][1210],size[120];//a为原数组 b为排序后数组 num为块数 next[i]为第i块的下一块序号 void rebuild(int x)//重构第x块 {    for(int i=1;i<=size[x];i++)b[x][i]=a[x][i];    sort(b[x]+1,b[x]+1+size[x]);}int binary(int x,int v)//第x块有多少数比v小 {    int l=1,r=size[x];    while(l<=r)    {        int mid=l+r>>1;        if(b[x][mid]<v)l=mid+1;        else r=mid-1;    }    return l-1;}int t1,t2;void get(int x)//t1表示x所属块的序号 t2表示x是块中的第几个 {    t1=1;t2=x;    while(t2>=size[t1])    {        if(t2==size[t1])break;        if(next[t1]==-1)break;        t2-=size[t1];        t1=next[t1];    }}int query(int L,int R,int k){    int bl,br,tl,tr;    get(L);bl=t1;tl=t2;    get(R);br=t1;tr=t2;    int l=0,r=70000;    while(l<=r)    {        int mid=l+r>>1,tot=0;        if(bl==br)        {            for(int i=tl;i<=tr;i++)            tot+=(a[bl][i]<mid);        }        else        {            for(int i=tl;i<=size[bl];i++)tot+=(a[bl][i]<mid);            for(int i=1;i<=tr;i++)tot+=(a[br][i]<mid);            int t=next[bl];            while(t!=br)            {                tot+=binary(t,mid);                t=next[t];            }        }        if(tot<=k-1)l=mid+1;        else r=mid-1;    }    return l-1;}void modify(int x,int v){    get(x);int t=t1;x=t2;    int i;    for(i=1;i<=size[t];i++)    if(a[t][x]==b[t][i])break;    a[t][x]=b[t][i]=v;    while(i<size[t]&&b[t][i]>b[t][i+1])swap(b[t][i],b[t][i+1]),i++;    while(i>1&&b[t][i]<b[t][i-1])swap(b[t][i],b[t][i-1]),i--;}void insert(int x,int v){    get(x);int t=t1;x=t2;    for(int i=size[t]+1;i>x;i--)a[t][i]=a[t][i-1];    a[t][x]=b[t][++size[t]]=v;    if(size[t]==2*block_size)    {        next[++num]=next[t];next[t]=num;        size[t]=size[num]=block_size;        for(int i=1;i<=block_size;i++)a[num][i]=a[t][i+block_size];        rebuild(num);rebuild(t);    }    else    {        int i=size[t];        while(i>1&&b[t][i]<b[t][i-1])swap(b[t][i],b[t][i-1]),i--;    }}int main(){    memset(next,-1,sizeof(next));    n=read();num=(n-1)/block_size+1;    for(int i=1;i<=n;i++)    {        int t1=(i-1)/block_size+1,t2=(i-1)%block_size+1;        a[t1][t2]=b[t1][t2]=read();        size[t1]++;    }    for(int i=1;i<=num;i++)    {        rebuild(i);        if(i!=num)next[i]=i+1;    }    int zaixian=1;    int ans=0,Q=read();    while(Q--)    {        char str[3];        scanf("%s",str);        int x=read()^(ans*zaixian),y=read()^(ans*zaixian);        if(str[0]=='Q')ans=query(x,y,read()^(ans*zaixian)),printf("%d\n",ans);        else if(str[0]=='M')modify(x,y);        else insert(x,y);    }}