分块的概念
分块算法用于是区间信息维护与查询的算法
分块的思想是把数据分成若干块,为了时间复杂度的均摊,块的大小通常设置为 n \sqrt n n,以块为单位维护区间的值
对于单点更新,时间复杂度 O ( 1 ) O(1) O(1)
区间更新,若区间不跨越块,直接暴力处理,时间复杂度 O ( n ) O(\sqrt n) O(n),若跨越块,则不能凑成一个整块的暴力处理,整块的直接更新,时间复杂度 O ( n ) O(\sqrt n) O(n)
区间查询同理
分块的特点:分块算法比较简洁,能够维护一些树状数组和线段树无法维护的值,但效率较低,但 1 0 5 10^5 105 级别足矣
疑问:最后一块不足整块的区间会不会影响查询?
答:不会,我们可以虚拟地把后端补齐
例题
Acwing 243 传送门
本题就是区间加法更新,然后区间求和查询,这里的做法是分块,时间复杂度
O
(
m
n
)
O(m\sqrt n)
O(mn)
#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
#include<string>
#include<set>
#include<map>
#include<unordered_map>
#include<queue>
#define mez(x) memset(x,0,sizeof x)
#define mei(x) memset(x,0x3f,sizeof x)
#define lowbit(x) -x&x
#define inf 0x3f3f3f3f
#define INF 0x7fffffff
#define cu(x) cout<<x<<"--"<<endl
#define ex exit(0)
#define pn puts("")
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
typedef vector<int> vc;
inline int read()
{
int x=0,f=1;char ch=getchar();
while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
while (ch>='0'&&ch<='9'){x=x*10+ch-48;ch=getchar();}
return x*f;
}
const int N=1e5+10,M=350;
int n,m,siz; //siz为分块大小
ll tag[M],sum[M]; //tag为分块的懒标记,sum为分块的和(不包含懒标记)
int val[N];
inline int ht(int& x) //数组下标对应的块下标
{
return (x-1)/siz; //为什么-1,由于数组下标从1开始
}
void update(int& l,int& r,int& w)
{
int i,j,k;
if(ht(l)==ht(r)) //若区间位于一个块中,则暴力求解
{
for(i=l;i<=r;++i)
val[i]+=w,sum[ht(l)]+=w;
}
else
{
//首先暴力收缩两端
i=l,j=r;
while(ht(i)==ht(l)) val[i]+=w,sum[ht(i)]+=w,++i;
while(ht(j)==ht(r)) val[j]+=w,sum[ht(j)]+=w,--j;
for(k=ht(i);k<=ht(j);++k)
tag[k]+=w,sum[k]+=(ll)siz*(ll)w;
}
}
ll query(int& l,int& r)
{
int i,j,k;
ll ans=0;
if(ht(l)==ht(r))
{
for(i=l;i<=r;++i)
ans+=val[i]+tag[ht(i)];
}
else
{
i=l,j=r;
while(ht(i)==ht(l)) ans+=val[i]+tag[ht(i)],++i;
while(ht(j)==ht(r)) ans+=val[j]+tag[ht(j)],--j;
for(k=ht(i);k<=ht(j);++k)
ans+=sum[k];
}
return ans;
}
int main()
{
int i,j,k,l,r,w;
n=read(),m=read();
siz=sqrt(n);
for(i=1;i<=n;++i)
val[i]=read(),sum[ht(i)]+=val[i];
char op[2];
while(m--)
{
scanf("%s",op);
if(op[0]=='C')
{
l=read(),r=read(),w=read();
update(l,r,w);
}
else
{
l=read(),r=read();
printf("%lld\n",query(l,r));
}
}
return 0;
}