【算法·分块】分块入门①-⑥

关于分块

关于分块，其主要思想是大局维护，局部朴素；即吧一个序列分成若干个小块，对任意区间进行操作的时候：如果区间内包含完整块，就进行快速的整体修改；若所修改区间内的一部分是某个个块中的一部分，则直接进行暴力修改就行。根据有关数学知识，可以证明分成$\sqrt{n}$的时间效率最高，一般的时间复杂度则是$O(n\sqrt{n})$

具体内容，我们可以通过六道分块习题进行理解。

分块入门1

给出一个长为n的数列，以及n个操作，操作涉及区间加法，单点查值。

这道题就是分块一个灰常简单的模板了，具体可以这么做：

• 设修改的区间为$[l,r]$，此时在分块中预处理了每一个区间的左端点和右端点，当前枚举为第i个区间。
• 若区间$i$在$[l,r]$以内，则直接用数组$add[i]$进行区间的整体操作，表示第$i$个区间的所有数的变化量，只需要累加一次。
• 若区间i的一部分在$[l,r]$以内，暴力修改即可。直接在a数组上修改。
• 对于查询的点x，输出的值为：$a[i]+add[po{s}_i]$,$po{s}_i$表示$i$所属的块的编号。
  代码如下：

（注意：代码可能与上述讲解的数组含义有所出入）

例如：上面的add是sum，而下面的add是单点变化量，等同于在原数组上做修改。

#include<bits/stdc++.h>
using namespace std;
int n,m;
int a[100000];
int L[100000];
int R[100000];
int add[100000];
int sum[100000];
int num[100000];
void Plus(int l,int r,int c)
{
   
	for (int i=1;i<=m;++i)
	{
   
		if (l>R[i] || r<L[i]) continue;
		if (l<=L[i] && r>=R[i]) sum[i]+=c;
		else 
		{
   
			for (int j=max(L[i],l);j<=min(R[i],r);++j) add[j]+=c;
		}
	}
}
void ask(int x)
{
   
	printf("%d\n",a[x]+sum[num[x]]+add[x]);
}
int main(void)
{
   
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout);
	scanf("%d",&n);
	for (int i=1;i<=n;++i) scanf("%d",a+i);
	m=sqrt(n);
	for (int i=1;i<=m;++i)
	{
   
		L[i]=sqrt(n)*(i-1)+1;
		R[i]=sqrt(n)*i;
	}	
	if (R[m]<n) 
	{
   
		m++;
		L[m]=R[m-1]+1;
		R[m]=n;
	}
	for (int i=1;i<=m;++i)
	    for (int j=L[i];j<=R[i];++j)
	        num[j]=i;
	for (int i=1,opt,l,r,c;i<=n;++i)
	{
   
		scanf("%d %d %d %d",&opt,&l,&r,&c);
		if (opt == 0) Plus(l,r,c);
		if (opt == 1) ask(r);
	}
	return 0;
}

分块入门2

给出一个长为n的数列，以及n个操作，操作涉及区间加法，询问区间内小于某个值x的元素个数。

区间加法的方法同分块1，我们来考虑如何求小于x的元素个数。

对于每一个块，我们可以在一开始分别进行排序，若不进行修改操作就可以直接用二分查找得到答案。

对于区间修改来说，如果是整块的修改则没有关系，因为仍然保证有序性。

如果是两端零零散散的修改，在修改完之后再对这些块分别进行排序。

这个就可以用过二分查找或者lowerbound函数在log的复杂度内解决问题了。

因为分成的块是logn块，时间复杂度是$O(n\sqrt{n}+2\ast$