【分块】bzoj1901 Zju2112 Dynamic Rankings

区间k大，分块大法好，每个区间内存储一个有序表。

二分答案，统计在区间内小于二分到的答案的值的个数，在每个整块内二分、零散的暴力即可。

还是说∵有二分操作，∴每个块的大小定为sqrt(n*log2(n))比较快呢。

 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 using namespace std;
 int n,a[],num[],qx,qy,k,p,v,m,sz,sum,l[],r[],b[],tmp[];
 char s[];
 void makeblock()
 {
     memcpy(b,a,sizeof(a));
     sz=sqrt((double)n*log2(n));
     for(sum=;sum*sz<n;sum++)
       {
         l[sum]=(sum-)*sz+;
         r[sum]=sum*sz;
         for(int i=l[sum];i<=r[sum];i++)
           num[i]=sum;
         sort(b+l[sum],b+r[sum]+);
       }
     l[sum]=sz*(sum-)+;
     r[sum]=n;
     sort(b+l[sum],b+r[sum]+);
     for(int i=l[sum];i<=r[sum];i++)
       num[i]=sum;
 }
 void query(int L,int R,int k)
 {
     if(num[L]+>=num[R])
       {
         int en=;
         for(int i=L;i<=R;i++) tmp[++en]=a[i];
         sort(tmp+,tmp+en+);
         printf("%d\n",tmp[k]);
       }
     else
       {
         int x=,y=;
         while(x!=y)
           {
             int mid=x+y>>,cnt=;
             for(int i=L;i<=r[num[L]];i++) if(a[i]<mid) cnt++;
             for(int i=l[num[R]];i<=R;i++) if(a[i]<mid) cnt++; //统计<mid的值数
             for(int i=num[L]+;i<num[R];i++) cnt+=( lower_bound(b+l[i],b+r[i]+,mid) - (b+l[i]) );
             if(cnt>=k) y=mid;
             else x=mid+;
           }
         printf("%d\n",x-);
       }
 }
 void update(int p,int v)
 {
     *lower_bound(b+l[num[p]],b+r[num[p]]+,a[p])=v;
     sort(b+l[num[p]],b+r[num[p]]+);
     a[p]=v;
 }
 int main()
 {
     scanf("%d%d",&n,&m);
     for(int i=;i<=n;i++) scanf("%d",&a[i]);
     makeblock();
     for(int i=;i<=m;i++)
       {
         scanf("%s",s);
         if(s[]=='Q') {scanf("%d%d%d",&qx,&qy,&k);query(qx,qy,k);}
         else {scanf("%d%d",&p,&v);update(p,v);}
       }
     return ;
 }