数列分块入门

数列分块就是把数组分成等大的区域便于统一操作。可以说是优雅的暴力了。

具体思路源自

http://www.cnblogs.com/CHerish_OI/category/1176577.html

1.loj6277  区间修改+单点查询    

        将数组分成sqrt（n）长度的子数列，用pos数组记录每个点所属区间的标号，add数组记录每个分块所增加的数值。当要对区间[ l , r ]进行操作（+c）时，分开处理两头的【因为可能需要更改的范围并不能恰好充满区间】，再对中间整块的数列进行分块处理。

2.loj6278  区间加法+询问区间内小于某个元素x的元素的个数

还是一样分成sqrt（n）的长度的区间。加法同上。因为区间同时增加和减少并不影响整体的排序，因此在输入的时候，将分块另存入一个vector中，并预处理排序一下。在算结果是对于一个完整的区间直接使用lower_bound函数就可以了。

3.loj6279  区间加法+求序列中某特定元素的前驱【即小于x的最大值】

分区间和加法同上，求前驱时分开对两端求最大值，再使用set【因为set实际上是基于红黑树实现的】lower_bound直接查找。

这里只给出sol函数

void sol(int l,int r,int c)
{
    int ans=-1;
    for(int i=l;i<=min(pos[l]*block,r);i++)
if(a[i]+ad[pos[i]]<c)ans=max(ans,a[i]+ad[pos[i]]);
    if(pos[l]!=pos[r])
        for(int i=(pos[r]-1)*block+1;i<=r;i++)
            if(a[i]+ad[pos[i]]<c)
                ans=max(ans,a[i]+ad[pos[i]]);//单独处理两端的 
    for(int i=pos[l]+1;i<=pos[r]-1;i++)
    {
        int x=c-ad[i];
        set<int>::iterator it=s[i].lower_bound(x);//前向迭代器 
        if(it==s[i].begin())continue;
        it--;
        ans=max(ans,*it+ad[i]);
    }
    printf("%d\n",ans);

}

4.loj6280  区间加法+区间查询

这个其实可以直接线段树实现的。还是用打标记的思想，sol的时候还是分开处理两端，中间完全覆盖的部分直接用add*数组长度+sum即可。

5.loj6281  区间开方+区间查询

这里有一个小小的优化就是当不断的开方操作进行后，最后可能有0或者1出现，这个时候就不要再对其开方了直接跳过就可以了。

这里给出处理的函数

void reset(int x) {
if (v[x])return;
v[x]=1; 
sum[x]=0;
for (int i=(x-1)*block+1;i<=x*block;i++) 
{
a[i]=sqrt(a[i]);
sum[x]+=a[i];
if(a[i]>1)v[x] = 0;
}
}
void update(int l, int r) {
if (!v[pos[l]]) {
for (int i=l;i<=min(pos[l]*block,r);i++) 
{
sum[pos[l]]-=a[i];
a[i]=sqrt(a[i]);
sum[pos[l]]+=a[i];
}
v[pos[l]]=1;
for (int i=(pos[l]-1)*block+1;i<=pos[l]*block;i++)
if (a[i]>1)v[pos[l]]=0;
}
if (pos[l] == pos[r])return;
if (!v[pos[r]]) 
{
for (int i=findl(pos[r]);i<=r;i++) 
{
sum[pos[r]]-=a[i];
a[i]=sqrt(a[i]);
sum[pos[r]]+=a[i];
}
v[pos[r]]=1;
for (int i=findl(pos[r]);i<=findr(pos[r]);i++)
if (a[i]>0)v[pos[r]] = 0;
}
for (int i=pos[l]+1;i<=pos[r]-1;i++)reset(i);

}

不得不说这个真的饿十分暴力呢。

6.loj6282  单点插入单点查询