【分块】bzoj3196 Tyvj 1730 二逼平衡树

分块 或 树套树。

在每个块中维护一个有序表，查询时各种二分，全都是分块的经典操作，就不详细说了。

块的大小定为sqrt(n*log2(n))比较快。

 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 #include<cstring>
 using namespace std;
 int sum,sz,l[],r[],num[],a[],b[],x,y,k,n,m,op,tmp[];
 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;
     for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
 }
 inline void Update()
 {
     *lower_bound(b+l[num[x]],b+r[num[x]]+,a[x])=y;
     sort(b+l[num[x]],b+r[num[x]]+);
     a[x]=y;
 }
 inline void Rank()
 {
     int ans=;
     if(num[x]+>=num[y]) {for(int i=x;i<=y;i++) if(a[i]<k) ans++;}
     else
       {
         for(int i=x;i<=r[num[x]];i++) if(a[i]<k) ans++;
         for(int i=l[num[y]];i<=y;i++) if(a[i]<k) ans++;
         for(int i=num[x]+;i<num[y];i++) ans+=lower_bound(b+l[i],b+r[i]+,k)-(b+l[i]);
       }
     printf("%d\n",ans+);
 }
 inline void Kth()
 {
     if(num[x]+>=num[y])
       {
         int en=;
         for(int i=x;i<=y;i++) tmp[++en]=a[i];
         sort(tmp+,tmp+en+);
         printf("%d\n",tmp[k]);
       }
     else
       {
         int L=,R=;
         while(L!=R)
           {
             int mid=L+R>>,cnt=;
             for(int i=x;i<=r[num[x]];i++) if(a[i]<mid) cnt++;
             for(int i=l[num[y]];i<=y;i++) if(a[i]<mid) cnt++; //统计<mid的值数
             for(int i=num[x]+;i<num[y];i++) cnt+=lower_bound(b+l[i],b+r[i]+,mid)-(b+l[i]);
             if(cnt>=k) R=mid;
             else L=mid+;
           }
         printf("%d\n",L-);
       }
 }
 inline void Front()
 {
     int ans=-;int* p;
     if(num[x]+>=num[y]) {for(int i=x;i<=y;i++) if(a[i]>ans && a[i]<k) ans=a[i];}
     else
       {
         for(int i=x;i<=r[num[x]];i++) if(a[i]>ans && a[i]<k) ans=a[i];
         for(int i=l[num[y]];i<=y;i++) if(a[i]>ans && a[i]<k) ans=a[i];
         for(int i=num[x]+;i<num[y];i++)
           {
             p=lower_bound(b+l[i],b+r[i]+,k);
             if(p!=b+l[i] && *(p-)<k && *(p-)>ans) ans=*(p-);
           }
       }
     printf("%d\n",ans);
 }
 inline void Next()
 {
     int ans=;int* p;
     if(num[x]+>=num[y]) {for(int i=x;i<=y;i++) if(a[i]<ans && a[i]>k) ans=a[i];}
     else
       {
         for(int i=x;i<=r[num[x]];i++) if(a[i]<ans && a[i]>k) ans=a[i];
         for(int i=l[num[y]];i<=y;i++) if(a[i]<ans && a[i]>k) ans=a[i];
         for(int i=num[x]+;i<num[y];i++)
           {
             p=upper_bound(b+l[i],b+r[i]+,k);
             if(p!=b+r[i]+ && *p>k && *p<ans) ans=*p;
           }
       }
     printf("%d\n",ans);
 }
 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("%d%d%d",&op,&x,&y);
         if(op==) Update();
         else
           {
               scanf("%d",&k);
             if(op==) Rank();
             else if(op==) Kth();
             else if(op==) Front();
             else Next();
           }
       }
     return ;
 }