蒟蒻的分块入门

前言

发现分块这个东西已经渗透到经济生活的各个领域（大雾），所以我不得不学些简单的分块技巧，否则会被同机房的jar形绕懵圈的。

解析

首先%%%%hzwer:   http://hzwer.com/8053.html

对于一个有n个元素的数列，我们把它拆碎了，拆成除了最后一块以外，每块都是k的n/k（或许要+1）个块。这样对于某些块，我们可以整块操作，优化时间。

例题：洛谷P3372

区间操作：对于每一个区间，一般是由前面一个不完整的块，后面一个不完整的块，和中间一些完整的块组成的。所以对于前后不完整的块，暴力修改，并且记录每个块的元素和，中间完整的块，用lazy标记标记整个块里每个元素都要加几，就可以了！

代码

#include<iostream>#include<algorithm>#include<string>#include<vector>#include<cstdlib>#include<cstdio>#include<cstring>#include<queue>#include<map>#include<cmath>#include<climits>using namespace std;int n,m,sqn;int bj[100005];long long v[100005],ful[405],laz[405];void upda(int l,int r,long long num){int i,j,lim,kl;lim=min(bj[l]*sqn,r);for(i=l;i<=lim;i++){v[i]+=num;ful[bj[i]]+=num;}if(bj[l]!=bj[r]){kl=(bj[r]-1)*sqn+1;for(i=kl;i<=r;i++){v[i]+=num;ful[bj[i]]+=num;}}for(i=bj[l]+1;i<=bj[r]-1;i++)laz[i]+=num;//懒惰标记}long long find(int l,int r){int i,j,lim,kl;long long ans=0;lim=min(bj[l]*sqn,r);for(i=l;i<=lim;i++)ans+=v[i]+laz[bj[i]];if(bj[l]!=bj[r]){kl=(bj[r]-1)*sqn+1;for(j=kl;j<=r;j++)ans+=v[j]+laz[bj[j]];}for(i=bj[l]+1;i<=bj[r]-1;i++)ans+=ful[i]+sqn*laz[i];return ans;}int main(){   int i,j,k,x,y,z;   long long num;   scanf("%d%d",&n,&m);sqn=sqrt(n);   for(i=1;i<=n;i++){   scanf("%lld",&v[i]);bj[i]=(i-1)/sqn+1;//标记在哪一块   ful[bj[i]]+=v[i];//整块的元素和   }   for(i=1;i<=m;i++){   scanf("%d%d%d",&x,&y,&z);   if(x==1){   scanf("%lld",&num);upda(y,z,num);   }   else {printf("%lld\n",find(y,z));}   }    return 0; }