loj6278 数列分块入门题2

题意：支持区间加，询问区间中元素排名

维护两个域。一个域维护原序列，一个域维护快内排序序列。

每次修改后更新快内排序序列。

修改时O(sqrt(n)log(sqrt(n)))

询问时O(sqrt(n)log(sqrt(n)))

大概是这个量级吧

 #include <bits/stdc++.h>
 using namespace std;

 int bl[],bla,a[],a0[],b[],bi[],n;
 int t1,t2,t3,t4;

 void init() {
     bla=sqrt(n);
     for (int i=;i<=n;i++)
         bl[i]=(i-)/bla+,
         bi[i]=(i-)%bla+;
     for (int l=;l<=n;l+=bla) {
         for(int i=bl[l]*bla-bla+;i<=bl[l]*bla;i++)
             a[i]=a0[i];
         sort(a+bl[l]*bla-bla+,a+bl[l]*bla+);
     }
 }

 void add(int l,int r,int c){
     if(bl[l]==bl[r]) {
         for(int i=l;i<=r;i++)
             a0[i]+=c;
         for(int i=bl[l]*bla-bla+;i<=bl[l]*bla;i++)
             a[i]=a0[i];
         sort(a+bl[l]*bla-bla+,a+bl[l]*bla+);
     }
     else {
         for(int i=l;i<=bl[l]*bla;i++)
             a0[i]+=c;
         for(int i=bl[l]*bla-bla+;i<=bl[l]*bla;i++)
             a[i]=a0[i];
         sort(a+bl[l]*bla-bla+,a+bl[l]*bla+);
         for(int i=bl[r]*bla-bla+;i<=r;i++)
             a0[i]+=c;
         for(int i=bl[r]*bla-bla+;i<=bl[r]*bla;i++)
             a[i]=a0[i];
         sort(a+bl[r]*bla-bla+,a+bl[r]*bla+);
         for(int i=bl[l]+;i<bl[r];i++)
             b[i]+=c;
     }
 }

 int query(int l,int r,long long c){
     int ans=;
     if(bl[l]==bl[r]) {
         for(int i=l;i<=r;i++)
             ans+=(a0[i]<c-b[bl[l]]);
     }
     else {
         for(int i=l;i<=bl[l]*bla;i++)
             ans+=(a0[i]<c-b[bl[l]]);
         for(int i=bl[r]*bla-bla+;i<=r;i++)
             ans+=(a0[i]<c-b[bl[r]]);
         for(int i=bl[l]+;i<bl[r];i++)
             ans+=lower_bound(a+i*bla-bla+,a+i*bla+,c-b[i])-(a+i*bla-bla+);
     }
     return ans;
 }

 void printblock(){
     /*for(int i=1;i<=n/bla+(n%bla)?1:0;i++) {
         printf("block %d : %d\n",i,b[i]);
         printf(" a0: ");
         for(int j=1;j<=bla;j++) printf("%d ",a0[i*bla-bla+j]);
         printf("\n");
         printf(" a:  ");
         for(int j=1;j<=bla;j++) printf("%d ",a[i*bla-bla+j]);
         printf("\n");
     }*/
 }

 int main(){
     ios::sync_with_stdio(false);
     cin>>n;
     for(int i=;i<=n;i++)
         cin>>a0[i];
     init();
     printblock();
     for(int i=;i<=n;i++) {
         cin>>t1>>t2>>t3>>t4;
         if(t1==){
             add(t2,t3,t4);
         }
         else {
             printf("%d\n",query(t2,t3,t4*t4));
         }
         printblock();
     }
 }