COJ 0986 WZJ的数据结构（负十四） 区间动态k大

题解：哈哈哈我过了！！！主席树+树状数组写起来还真是hentai啊。。。

在这里必须分享我的一个沙茶错！！！看这段代码：

void get(int x,int d){
    if(!d)for(lt[ltot=]=root[x];x;x-=x&-x)if(c[x])lt[++ltot]=c[x];
    else  for(rt[rtot=]=root[x];x;x-=x&-x)if(c[x])rt[++rtot]=c[x];return;
}

似乎挺正常是吧？而且缩进也很优美对不对！！！

可是！！！！！！！！窝萌发现它是错的！！！！！！！！！！

因为下面那个else跟上面的if配对了！！！！！！所以这段程序其实是这样的：

void get(int x,int d){
    if(!d)for(lt[ltot=]=root[x];x;x-=x&-x){
        if(c[x])lt[++ltot]=c[x];
        else{
            for(rt[rtot=]=root[x];x;x-=x&-x)
            if(c[x])rt[++rtot]=c[x];
        }
    }return;
}

我要疯了！！！！！！！！！！！！！！！！！！！

 #include<iostream>
 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 #include<queue>
 #include<cstring>
 #define PAU putchar(' ')
 #define ENT putchar('\n')
 #define CH for(int d=0;d<2;d++)if(ch[d])
 #define lson ls(x),y->ch[0],L,M
 #define rson rs(x),y->ch[1],M+1,R
 using namespace std;
 const int maxn=+,maxnode=+,inf=-1u>>,vl=,vr=;
 struct node{
     node*ch[];int siz;node(){siz=;}
 }pol[maxnode],*nodecnt=pol,*root[maxn],*c[maxn],*lt[maxn],*rt[maxn];int ltot,rtot;
 node*ls(node*x){return x?x->ch[]:x;}
 node*rs(node*x){return x?x->ch[]:x;}
 int sz(node*x){return x?x->siz:;}
 int n,m,A[maxn],cv;
 void build(int v,int pos,node*x,node*&y,int L=vl,int R=vr){
     y=nodecnt++;y->siz=sz(x)+v;if(L==R)return;int M=L+R>>;
     if(pos<=M)y->ch[]=rs(x),build(v,pos,lson);else y->ch[]=ls(x),build(v,pos,rson);return;
 }
 void update(int x,int v){
     for(int w=x;w<=n;w+=w&-w)build(-,A[x],c[w],c[w]);
     for(int w=x;w<=n;w+=w&-w)build(,A[x]=v,c[w],c[w]);return;
 }
 void get(int x,int d){
     if(!d){for(lt[ltot=]=root[x];x;x-=x&-x)if(c[x])lt[++ltot]=c[x];}
     else  {for(rt[rtot=]=root[x];x;x-=x&-x)if(c[x])rt[++rtot]=c[x];}return;
 }
 inline int read(){
     int x=,sig=;char ch=getchar();
     for(;!isdigit(ch);ch=getchar())if(ch=='-')sig=;
     for(;isdigit(ch);ch=getchar())x=*x+ch-'';
     return sig?x:-x;
 }
 inline void write(int x){
     if(x==){putchar('');return;}if(x<)putchar('-'),x=-x;
     int len=,buf[];while(x)buf[len++]=x%,x/=;
     for(int i=len-;i>=;i--)putchar(buf[i]+'');return;
 }
 void init(){
     n=read();m=read();
     for(int i=;i<=n;i++)build(,A[i]=read(),root[i-],root[i]);
     return;
 }
 void work(){
     int x,y,v;
     while(m--){
         if(!read())x=read(),v=read(),update(x,v);
         else{x=read();y=read();v=read();
             get(x-,);get(y,);int L=vl,R=vr,M,kth;
             while(L<R){M=L+R>>;kth=;//printf("%d %d %d %d\n",L,R,ltot,rtot);
                 for(int i=;i<=ltot;i++)kth-=sz(ls(lt[i]));
                 for(int i=;i<=rtot;i++)kth+=sz(ls(rt[i]));
                 if(kth>=v){R=M;
                     for(int i=;i<=ltot;i++)lt[i]=ls(lt[i]);
                     for(int i=;i<=rtot;i++)rt[i]=ls(rt[i]);
                 }else{L=M+;v-=kth;
                     for(int i=;i<=ltot;i++)lt[i]=rs(lt[i]);
                     for(int i=;i<=rtot;i++)rt[i]=rs(rt[i]);
                 }
             }write(L);ENT;
         }
     }
     return;
 }
 void print(){
     return;
 }
 int main(){init();work();print();return ;}