算法基础——分块

分块的概念

分块算法用于是区间信息维护与查询的算法

分块的思想是把数据分成若干块，为了时间复杂度的均摊，块的大小通常设置为 $\sqrt{n}$，以块为单位维护区间的值

对于单点更新，时间复杂度 $O(1)$

区间更新，若区间不跨越块，直接暴力处理，时间复杂度 $O(\sqrt{n})$，若跨越块，则不能凑成一个整块的暴力处理，整块的直接更新，时间复杂度 $O(\sqrt{n})$

区间查询同理

分块的特点：分块算法比较简洁，能够维护一些树状数组和线段树无法维护的值，但效率较低，但 $10^5$ 级别足矣

疑问：最后一块不足整块的区间会不会影响查询？
答：不会，我们可以虚拟地把后端补齐

例题

Acwing 243 传送门
本题就是区间加法更新，然后区间求和查询，这里的做法是分块，时间复杂度 $O(m\sqrt{n})$

#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;
}