loj6278 数列分块入门题2

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

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

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

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

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

大概是这个量级吧

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

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

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

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+1;i<=bl[l]*bla;i++)
            a[i]=a0[i];
        sort(a+bl[l]*bla-bla+1,a+bl[l]*bla+1);
    }
    else {
        for(int i=l;i<=bl[l]*bla;i++)
            a0[i]+=c;
        for(int i=bl[l]*bla-bla+1;i<=bl[l]*bla;i++)
            a[i]=a0[i];
        sort(a+bl[l]*bla-bla+1,a+bl[l]*bla+1);
        for(int i=bl[r]*bla-bla+1;i<=r;i++)
            a0[i]+=c;
        for(int i=bl[r]*bla-bla+1;i<=bl[r]*bla;i++)
            a[i]=a0[i];
        sort(a+bl[r]*bla-bla+1,a+bl[r]*bla+1);
        for(int i=bl[l]+1;i<bl[r];i++)
            b[i]+=c;
    }
}

int query(int l,int r,long long c){
    int ans=0;
    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+1;i<=r;i++)
            ans+=(a0[i]<c-b[bl[r]]);
        for(int i=bl[l]+1;i<bl[r];i++)
            ans+=lower_bound(a+i*bla-bla+1,a+i*bla+1,c-b[i])-(a+i*bla-bla+1);
    }
    return ans;
}

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

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