蒟蒻的分块入门

前言

发现分块这个东西已经渗透到经济生活的各个领域（大雾），所以我不得不学些简单的分块技巧，否则会被同机房的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; 
}