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

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

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

void get(int x,int d){
    if(!d)for(lt[ltot=1]=root[x];x;x-=x&-x)if(c[x])lt[++ltot]=c[x];
    else  for(rt[rtot=1]=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=1]=root[x];x;x-=x&-x){
        if(c[x])lt[++ltot]=c[x];
        else{
            for(rt[rtot=1]=root[x];x;x-=x&-x)
            if(c[x])rt[++rtot]=c[x];
        }
    }return;
}

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

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cmath>
 4 #include<algorithm>
 5 #include<queue>
 6 #include<cstring>
 7 #define PAU putchar(‘ ‘)
 8 #define ENT putchar(‘\n‘)
 9 #define CH for(int d=0;d<2;d++)if(ch[d])
10 #define lson ls(x),y->ch[0],L,M
11 #define rson rs(x),y->ch[1],M+1,R
12 using namespace std;
13 const int maxn=100000+10,maxnode=20000000+10,inf=-1u>>1,vl=1,vr=100000;
14 struct node{
15     node*ch[2];int siz;node(){siz=0;}
16 }pol[maxnode],*nodecnt=pol,*root[maxn],*c[maxn],*lt[maxn],*rt[maxn];int ltot,rtot;
17 node*ls(node*x){return x?x->ch[0]:x;}
18 node*rs(node*x){return x?x->ch[1]:x;}
19 int sz(node*x){return x?x->siz:0;}
20 int n,m,A[maxn],cv;
21 void build(int v,int pos,node*x,node*&y,int L=vl,int R=vr){
22     y=nodecnt++;y->siz=sz(x)+v;if(L==R)return;int M=L+R>>1;
23     if(pos<=M)y->ch[1]=rs(x),build(v,pos,lson);else y->ch[0]=ls(x),build(v,pos,rson);return;
24 }
25 void update(int x,int v){
26     for(int w=x;w<=n;w+=w&-w)build(-1,A[x],c[w],c[w]);
27     for(int w=x;w<=n;w+=w&-w)build(1,A[x]=v,c[w],c[w]);return;
28 }
29 void get(int x,int d){
30     if(!d){for(lt[ltot=1]=root[x];x;x-=x&-x)if(c[x])lt[++ltot]=c[x];}
31     else  {for(rt[rtot=1]=root[x];x;x-=x&-x)if(c[x])rt[++rtot]=c[x];}return;
32 }
33 inline int read(){
34     int x=0,sig=1;char ch=getchar();
35     for(;!isdigit(ch);ch=getchar())if(ch==‘-‘)sig=0;
36     for(;isdigit(ch);ch=getchar())x=10*x+ch-‘0‘;
37     return sig?x:-x;
38 }
39 inline void write(int x){
40     if(x==0){putchar(‘0‘);return;}if(x<0)putchar(‘-‘),x=-x;
41     int len=0,buf[15];while(x)buf[len++]=x%10,x/=10;
42     for(int i=len-1;i>=0;i--)putchar(buf[i]+‘0‘);return;
43 }
44 void init(){
45     n=read();m=read();
46     for(int i=1;i<=n;i++)build(1,A[i]=read(),root[i-1],root[i]);
47     return;
48 }
49 void work(){
50     int x,y,v;
51     while(m--){
52         if(!read())x=read(),v=read(),update(x,v);
53         else{x=read();y=read();v=read();
54             get(x-1,0);get(y,1);int L=vl,R=vr,M,kth;
55             while(L<R){M=L+R>>1;kth=0;//printf("%d %d %d %d\n",L,R,ltot,rtot);
56                 for(int i=1;i<=ltot;i++)kth-=sz(ls(lt[i]));
57                 for(int i=1;i<=rtot;i++)kth+=sz(ls(rt[i]));
58                 if(kth>=v){R=M;
59                     for(int i=1;i<=ltot;i++)lt[i]=ls(lt[i]);
60                     for(int i=1;i<=rtot;i++)rt[i]=ls(rt[i]);
61                 }else{L=M+1;v-=kth;
62                     for(int i=1;i<=ltot;i++)lt[i]=rs(lt[i]);
63                     for(int i=1;i<=rtot;i++)rt[i]=rs(rt[i]);
64                 }
65             }write(L);ENT;
66         }
67     }
68     return;
69 }
70 void print(){
71     return;
72 }
73 int main(){init();work();print();return 0;}