数据结构进阶之分块

入坑信奥的L同学

已于 2022-11-10 22:04:14 修改

阅读量157

点赞数

分类专栏：笔记文章标签：数据结构分块

于 2022-10-10 21:46:09 首次发布

本文链接：https://blog.csdn.net/m0_55341679/article/details/126504401

版权

分块算法在解决不满足“区间可加性”问题时展现出优势，是线段树和树状数组的有效补充。文章通过四个例题详细讲解了分块算法的应用，包括区间和、众数维护、磁力块问题和莫队算法，展示了分块在不同场景下的解决方案和优化技巧。

摘要由CSDN通过智能技术生成

分块简介

分块算法被称作“优美的暴力”，是维护法与暴力法的完美结合。对于线段树、树状数组所解决不了的不满足“区间可加性”的问题时，分块就显得及其优秀，而对于它们都可以解决的题目，分块的码长居中，思维难度居中，有时也可以使得问题简单很多，从而脱颖而出。
例题0：AcWing 243.一个简单的整数问题2
是的，这题有三种方法可以做：树状数组、线段树、分块。分块的基本思想是：大段维护，局部朴素。也就是将数组划分为多段，对于完全包含的段进行维护，而对于不满一段的通过暴力枚举进行计算。这题用分块可以比较轻松地解决，块内存储区间和即可，代码如下：

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
typedef long long ll;
const int N=1e5+5,M=320;
ll sum[M],add[M];
int L[M],a[N],R[M],cnt=0;
int pos[N];
//pos用于记录每个点所处的块的编号；sum记录块中数的和；
//add是块的标记，表示整个块所有数都发生了add[i]大小的偏移
void modify(int l,int r,int d){
   
    int p=pos[l],q=pos[r];
    if(p==q){
   
        for(int i=l;i<=r;i++) a[i]+=d;
        sum[p]+=d*(r-l+1);
    }
    else{
   
        for(int i=p+1;i<=q-1;i++) add[i]+=d;
        for(int i=l;i<=R[p];i++) a[i]+=d;
        sum[p]+=(R[p]-l+1)*d;
        for(int i=L[q];i<=r;i++) a[i]+=d;
        sum[q]+=(r-L[q]+1)*d;
    }
}
ll query(int l,int r){
   
    int p=pos[l],q=pos[r];
    ll res=0;
    if(p==q){
   
        for(int i=l;i<=r;i++) res+=a[i];
        res+=add[p]*(r-l+1);
    }
    else{
   
        for(int i=p+1;i<=q-1;i++) res+=sum[i]+add[i]*(R[i]-L[i]+1);
        for(int i=l;i<=R[p];i++) res+=a[i];
        res+=(R[p]-l+1)*add[p];
        for(int i=L[q];i<=r;i++) res+=a[i];
        res+=(r-L[q]+1)*add[q];
    }
    return res;
}
int main(){
   
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    int t=(int)sqrt(n);
    for(int i=1;i<=n;i+=t){
   
        L[++cnt]=i;
        R[cnt]=t+i-1;
    }
    if(R[cnt]<n){
   //多余部分
        cnt++;
        L[cnt]=R[cnt-1]+1;
        R[cnt]=n;
    }
    for(int i=1;i<=cnt;i++){
   
        for(int j=L[i];j<=R[i];j++){
   
            pos[j]=i;
            sum[i]+=a[j];
        }
    }
    while(m--){
   
        char op;
        int l,r;
        cin>>op>>l>>r;
        if(op=='C'){