浅谈分块算法

一直觉得分块是一个很高端的东西…一直没敢碰，现在才知道分块就是一种稍微优美一些的暴力，所以没有学过分块的同学不要害怕啦…

---

首先，分块是什么意思呢，顾名思义就是把要处理的东西进行分块，分成一块一块的233，举个很简单的例子，对于一个数列 size（a{ }）=5，我们可以把前2个分到一起，再两个分到一起，最后单下来一个，为什么要这样处理呢？这样处理的好处又是什么呢？
我们也可以这样思考，如果我们把一个数列，当该数列的长度为n的时候，我们以根号n为一段，分出来的段数不超过根号n，如果我们要进行区间的处理，比如加法减法等，可以对于修改区间[ L , R ]可以把其中框起来的块（一块是根号n的大小）直接打上标记，由于每一块的长度不大于根号n，所以对于两边没有框起来的部分，我们直接暴力地进行更新，这样操作次数是最多2倍根号n的，而中间的标记是O（1）处理的，这就是为什么该类算法是根号级别的原因

---

思路看起来很清晰吧，好像很简单的样子，那么我们马上就来试试吧

---

以黄学长的分块学习顺序为例
（本文的所有题目，均来自黄学长，www.hzwer.com）

我们先进行区间修改单点查询操作

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#define MAXN 100000
using namespace std;
int a[MAXN],add[MAXN],b[MAXN],len,n,m;//add为标记
void modify(int l,int r,int ad){
    for(register int i=l;i<=min(b[l]*len,r);i++) a[i]+=ad;
    if(b[l]!=b[r]){
  //注意此处的min操作，可能会在小数据卡到你
        for(register int i=(b[r]-1)*len+1;i<=r;i++) a[i]+=ad; 
    }
    for(register int i=b[l]+1;i<=b[r]-1;i++)add[i]+=ad;
}
int main(){
    scanf("%d%d",&n,&m);
    len=sqrt(n);
    for(register int i=1;i<=n;i++)scanf("%d",&a[i]);
    for(register int i=1;i<=n;i++)b[i]=(i-1)/len+1;
    for(register int i=1;i<=m;i++){
        int temp;
        scanf("%d",&temp);
        if(temp==1){
            int l,r,ad;
            scanf("%d%d%d",&l,&r,&ad);
            modify(l,r,ad);
        }else{
            int loc;
            scanf("%d",&loc);
            printf("%d\n",a[loc]+add[b[loc]]);
        }//该点的值和该店所在的区间累加的标记
    }
    return 0;
}

---

接下来是区间修改，区间查询小于某值的数有多少个

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#include<cmath> 
#define MAXN 1000000+10
#define MINN 5000
using namespace std;
vector<int>v[MINN];
int len,a[MAXN],b[MAXN],add[MINN];
void reset(int zone){
    v[zone].clear();
    for(register int i=(zone-1)*len+1;i<=zone*len;i++) v[zone].push_back(a[i]);
    sort(v[zone].begin(),v[zone].end());
}
void modify(int l,int r,int ad){
    for(register int i=l;i<=min(b[l]*len,r);i++) a[i]+=ad;
    reset(b[l]);
    if(b[l]!=b[r]){
        for(register int i=(b[r]-1)*len+1;i<=r;i++) a[i]+=ad;
        reset(b[r]);
    }
    for(register int i=b[l]+1;i<=b[r]-1;i++) add[i]+=ad;
}
int query(int l,int r,int k){
    int cnt=0;
    for(register int i=l;i<=min(b[l]*len,r);i++){
        if(a[i]+add[b[i]]<k)cnt++;
    }
    if(b[l]!=b[r]){
        for(register int i=(b[r]-1)*len+1;i<=r;i++){
            if(a[i]+add[b[i]]<k)cnt++;
        }
    }
    for(register int i=b[l]+1;i<=b[r]-1;i++){
        int x=k-add[i];
        cnt+=lower_bound(v[i].begin(),v[i].end(),x)-v[i].begin();
    }
    return cnt;
}
int main(){
    int n,m;
    scanf("%d%d",&n,&m);
    len=sqrt(n);
    for(register int i=1;i<=n;i++) scanf("%d",&a[i]);
    for(register int i=1;i<=n;i++){
        b[i]=(i-1)/len+1;
        v[b[i]].push_back(a[i]);
    }
    for(register int i=1;i<=b[n];i++) sort(v[i].begin(),v[i].end());
    for(register