分块入门

分块——一巴掌把题拍扁的算法

就是把序列分成块,每块内部暴力处理,外部打标记的算法

时间往往是n*(m+n/m)很容易得出当m=根号n是最小

所以n<=1e5时往往适用

 

ORZ hzwer

https://loj.ac/problem/6277

题目描述

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

输入格式

第一行输入一个数字 nnn。

第二行输入 nnn 个数字,第 iii 个数字为 aia_ia​i​​,以空格隔开。

接下来输入 nnn 行询问,每行输入四个数字 opt\mathrm{opt}opt、lll、rrr、ccc,以空格隔开。

若 opt=0\mathrm{opt} = 0opt=0,表示将位于 [l,r][l, r][l,r] 的之间的数字都加 ccc。

若 opt=1\mathrm{opt} = 1opt=1,表示询问 ara_ra​r​​ 的值(lll 和 ccc 忽略)。

输出格式

对于每次询问,输出一行一个数字表示答案。

样例

样例输入

4
1 2 2 3
0 1 3 1
1 0 1 0
0 1 2 2
1 0 2 0

样例输出

2
5

数据范围与提示

对于 100% 100\%100% 的数据,1≤n≤50000,−231≤others 1 \leq n \leq 50000, -2^{31} \leq \mathrm{others}1≤n≤50000,−2​31​​≤others、ans≤231−1 \mathrm{ans} \leq 2^{31}-1ans≤2​31​​−1。

 

吐槽一下:从输出描述可以看出这题目出的多随意)

如果加的区间包含块,就标记下,如果不包含,暴力修改

时间O(n根号n),优秀

#include<cstdio>
#include<cmath>
#define ll long long
using namespace std;

int read()
{
	int ret=0; bool f=0;
	char ch=getchar();
	while(ch<'0'||ch>'9') 
	{
		if(ch=='-') f=1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9')
		ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar();
	return f?-ret:ret;
}

const int N=1e5+5;
int id[N],n;

struct A
{
	int num;
	ll a[N],t[N];
	inline void pre()
	{
		num=sqrt(n);
		for(int i=1;i<=n;i++)
			id[i]=(i-1)/num+1,a[i]=read();
	}	
	
	inline void add(int l,int r,int c)
	{
		for(int i=id[l]+1;i<=id[r]-1;i++) t[i]+=c;
		if(id[l]==id[r]) 
		{
			for(int i=l;i<=r;i++) a[i]+=c;
			return;
		}
		for(int i=l;i<=num*id[l];i++) 
			a[i]+=c;
		for(int i=num*(id[r]-1)+1;i<=r;i++)
			a[i]+=c;
	}
	
	inline ll sum(int x) 
	{
		return a[x]+t[id[x]];
	}
}fk;

int main()
{
	n=read();
	fk.pre();
	for(int i=1;i<=n;i++)
	{
		int t=read(),l=read(),r=read(),c=read();
		if(!t) fk.add(l,r,c);
			else printf("%lld\n",fk.sum(r));
	}
	return 0;
}

 

### 数列分块入门第8题的算法实现与解析 数列分块是一种高效的处理区间查询和修改的技术,其核心思想是将数组划分为若干个连续的小块,每一块内的元素可以快速更新或查询。对于数列分块入门第8题,假设问题是涉及区间的加法操作以及最大值查询,则可以通过以下方法来解决。 #### 1. 数据结构设计 为了高效完成区间加法和最大值查询的操作,我们可以维护两个辅助数组: - `block_sum[]`:存储每个块的最大值。 - `lazy_tag[]`:标记每个块是否有延迟更新(即尚未应用到具体元素上的增量)。 这些数据结构的设计使得我们可以在 $ O(\sqrt{n}) $ 的时间复杂度下完成单次操作[^3]。 #### 2. 初始化过程 初始化时,我们需要计算初始状态下的块划分情况,并填充上述辅助数组的内容。以下是具体的代码实现: ```cpp #include <bits/stdc++.h> using namespace std; const int MAXN = 1e5 + 5; int a[MAXN], block_size, block_num; long long block_max[400]; bool lazy_flag[400]; void build_block(int n) { block_size = sqrt(n); block_num = (n + block_size - 1) / block_size; memset(block_max, 0, sizeof(block_max)); memset(lazy_flag, false, sizeof(lazy_flag)); for (int i = 0; i < n; ++i) { int idx = i / block_size; block_max[idx] = max(block_max[idx], (long long)a[i]); } } ``` #### 3. 延迟标记的应用 当执行区间加法时,为了避免逐一遍历整个范围中的每一个元素,引入懒惰传播机制。如果某个整块完全被覆盖在当前操作范围内,则直接对该块打上标签并记录增量;否则逐一访问该块内部受影响的部分。 下面是针对这一逻辑的具体函数定义: ```cpp // Apply the pending update to all elements within specified block. void propagate(int blk_idx, int delta) { if (!lazy_flag[blk_idx]) return; // Update maximum value of this block accordingly. block_max[blk_idx] += delta * block_size; lazy_flag[blk_idx] = false; } // Add 'delta' to range [l, r]. void add_range(int l, int r, int delta, int n) { int start_blk = l / block_size, end_blk = r / block_size; if (start_blk == end_blk) { for (int i = l; i <= min(r, (start_blk + 1) * block_size - 1); ++i) a[i] += delta; // Recalculate new maximum after modification. block_max[start_blk] = 0; for (int i = start_blk * block_size; i < ((start_blk + 1) * block_size && i < n); ++i) block_max[start_blk] = max((long long)a[i], block_max[start_blk]); } else { // Process first incomplete block separately. for (int i = l; i < (start_blk + 1) * block_size; ++i) a[i] += delta; block_max[start_blk] = 0; for (int i = start_blk * block_size; i < ((start_blk + 1) * block_size && i < n); ++i) block_max[start_blk] = max((long long)a[i], block_max[start_blk]); // Fully covered blocks can simply apply tag updates. for (int b = start_blk + 1; b < end_blk; ++b){ block_max[b] += delta * block_size; lazy_flag[b] |= true; } // Handle last partial block similarly as above case. for (int i = end_blk * block_size; i <= r; ++i) a[i] += delta; block_max[end_blk] = 0; for (int i = end_blk * block_size; i < ((end_blk + 1) * block_size && i < n); ++i) block_max[end_blk] = max((long long)a[i], block_max[end_blk]); } } ``` #### 4. 查询最值功能 最后一步是在给定区间内查找最大的数值。这同样依赖于之前构建好的块级信息来进行加速检索。 ```cpp long long query_max(int l, int r) { int start_blk = l / block_size, end_blk = r / block_size; long long result = LLONG_MIN; if (start_blk == end_blk) { for (int i = l; i <= r; ++i) result = max(result, (long long)a[i]); } else { for (int i = l; i < (start_blk + 1) * block_size; ++i) result = max(result, (long long)a[i]); for (int b = start_blk + 1; b < end_blk; ++b) result = max(result, block_max[b]); for (int i = end_blk * block_size; i <= r; ++i) result = max(result, (long long)a[i]); } return result; } ``` 通过以上步骤即可有效应对数列分块相关的题目需求。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值