数列分块2

题意：给出一个长为n的数列，以及n个操作，操作涉及区间加法，询问区间内小于某个值x的元素个数。 

有了上一题的经验，我们可以发现，数列简单分块问题实际上有三项东西要我们思考：

对于每次区间操作：

1.不完整的块 的O(√n)个元素怎么处理？

2.O(√n)个 整块 怎么处理？

3.要预处理什么信息（复杂度不能超过后面的操作）？

 

我们先来思考只有询问操作的情况，不完整的块枚举统计即可；而要在每个整块内寻找小于一个值的元素数，于是我们不得不要求块内元素是有序的，这样就能使用二分法对块内查询，需要预处理时每块做一遍排序，复杂度O(nlogn)，每次查询在√n个块内二分，以及暴力2√n个元素，总复杂度O(nlogn + n√nlog√n)。

可以通过均值不等式计算出更优的分块大小，就不展开讨论了

 

那么区间加怎么办呢？

套用第一题的方法，维护一个加法标记，略有区别的地方在于，不完整的块修改后可能会使得该块内数字乱序，所以头尾两个不完整块需要重新排序，复杂度分析略。

在加法标记下的询问操作，块外还是暴力，查询小于（x – 加法标记）的元素个数，块内用（x – 加法标记）作为二分的值即可。

#include <bits/stdc++.h>
using namespace std;
const int maxn = 50005;
typedef long long LL;
vector<LL> ve[maxn];
LL a[maxn],laz[maxn];
int bl[maxn],st[maxn],ed[maxn];
int block,n;
void force( int l,int r,LL v ){
    for( int i =l ;i <= r;i++ ) a[i] += v;
}
void reset( int x ){
    ve[x].clear();
    for( int i = st[x];i <= ed[x];i++ ) ve[x].push_back(a[i]);
    sort( ve[x].begin(),ve[x].end() );
}
void update( int l,int r,LL v ){
    for( int i = l;i <= min( ed[bl[l]],r );i++ )a[i] += v;
    reset( bl[l] );
    if( bl[l] != bl[r] ){
        for( int i = st[bl[r]];i <= r;i++ ) a[i] += v;
        reset( bl[r] );
    }
    for( int i = bl[l]+1;i < bl[r];i++ ){
        laz[i] += v;
    }
}
int query( int l,int r,LL v ){
    int res = 0;
    for( int i = l;i <= min( ed[bl[l]],r );i++ ){
        if( a[i] + laz[ bl[l] ] < v ) res++;
    }
    if( bl[l] !=bl[r] ){
        for( int i = st[ bl[r] ];i <= r;i++ ){
            if( a[i] + laz[ bl[r] ] < v ) res++;
        }
    }
    for( int i = bl[l]+1;i < bl[r];i++ ){
        LL x = v-laz[i];
        res += lower_bound( ve[i].begin(),ve[i].end(),x ) -ve[i].begin();
    }
    return res;
}
int main(){
    //freopen( "123.txt","r",stdin );
    //freopen("456.txt","w",stdout);
    int op,l,r,c;
    scanf("%d",&n);
    block = sqrt(n);
    for( int i = 1;i <= n;i++ )scanf("%lld",&a[i]);
    for( int i = 1;i <= n;i++ ){
        bl[i] = (i-1)/block+1;
        ve[ bl[i] ].push_back( a[i] );
        if( bl[i] != bl[i-1] ){
            st[ bl[i] ] = i;
            ed[ bl[i-1] ] = i-1;
        }
    }
    ed[ bl[n] ] = n;
    for( int i = 1;i <= bl[n];i++ ){
        //sort( a+st[i],a+ed[i]+1 );
        sort( ve[i].begin(),ve[i].end() );
    }
    for( int i = 1;i <= n;i++ ){
        scanf("%d%d%d%d",&op,&l,&r,&c);
        l = max( l,l );r = min( r,n );
        if( l > r ) swap(l,r);
        if( op == 0 ){
            update( l,r,c );
        }else{
            int cnt = query( l,r,1LL*c*c );
            printf("%d\n",cnt);
        }
    }
    return 0;
}