关于分块
关于分块,其主要思想是大局维护,局部朴素;即吧一个序列分成若干个小块,对任意区间进行操作的时候:如果区间内包含完整块,就进行快速的整体修改;若所修改区间内的一部分是某个个块中的一部分,则直接进行暴力修改就行。根据有关数学知识,可以证明分成 n \sqrt n n的时间效率最高,一般的时间复杂度则是 O ( n n ) O(n\sqrt n) O(nn)
分块入门1
给出一个长为n的数列,以及n个操作,操作涉及区间加法,单点查值。
这道题就是分块一个灰常简单的模板了,具体可以这么做:
- 设修改的区间为 [ l , r ] [l,r] [l,r],此时在分块中预处理了每一个区间的左端点和右端点,当前枚举为第i个区间。
- 若区间 i i i在 [ l , r ] [l,r] [l,r]以内,则直接用数组 a d d [ i ] add[i] add[i]进行区间的整体操作,表示第 i i i个区间的所有数的变化量,只需要累加一次。
- 若区间i的一部分在 [ l , r ] [l,r] [l,r]以内,暴力修改即可。直接在a数组上修改。
- 对于查询的点x,输出的值为: a [ i ] + a d d [ p o s i ] a[i]+add[pos_i] a[i]+add[posi], p o s i pos_i posi表示 i 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 n + 2 ∗ n l o g n ) O(n\sqrt n+2*\sqrt n \ log\sqrt n) O(nn+2∗n logn)
代码如下:
#include<bits/stdc++.h>
#define find lower_bound
using namespace std;
int n,m;
int a[100000];
int L[100000];
int R[100000];
int pos[100000];
int sum[100000];
vector<int>num[100000];
void reset(int x)
{
num[x].clear();
for (int i=L[x];i<=R[x];++i)
num[x].push_back(a[i]);
sort(num[x].begin(),num[x].end());
}
void change(int l,int r,int v)
{
for (int i=1;i<=m;++i)
{
if (l>R[i] || r<L[i]) continue;
if (l<=L[i] && r>=R[i]) sum[i]+=v;
else
{
for (int j=max(l,L[i]);j<=min(R[i],r);++j)
a[j]+=v;
reset(i);
}
}
}
void print(int l,int r,int k)
{
int ans=0;
for (int i=1;i<=m;++i)
{
if (l>R[i] || r<L[i]) continue;
if (l<=L[i] && r>=R[i])
ans += find (num[i].begin(), num[i].end(), k*k-sum[i]) - num[i].begin();
else
for (int j=max(l,L[i]);j<=min(R[i],r);++j)
if (a[j]+sum[pos[j]]<k*k) ans++;
}
printf("%d\n",ans);
}
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