一直觉得分块是一个很高端的东西…一直没敢碰,现在才知道分块就是一种稍微优美一些的暴力,所以没有学过分块的同学不要害怕啦…
首先,分块是什么意思呢,顾名思义就是把要处理的东西进行分块,分成一块一块的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