分块入门思想

这篇博客探讨了数列分块思想,主要应用于处理长度为n的数列,通过将数列均分为根号n个块,采用线段树的lazy懒惰标记策略进行区间操作。文章提到了两种分块方法:函数法和数组法,并指出函数法常数较大,推荐使用数组法。还详细介绍了分块函数和操作函数的实现,并提供了一个例题来说明如何在实际问题中应用这些概念。
摘要由CSDN通过智能技术生成

主要是对于分块思想的一些个人理解

主要思想

最主要的思想是把一个长度为n数列,均等的分成根号n个块,在对区间进行整体操作的时候,遇到完整的块,就直接对整个块进行操作,需要用到线段树lazy懒惰标记的区间处理的办法,对于两端不完整的块,直接去暴力遍历它。

对于数列分块的分块处理

主要是有两种方法每一种是函数法每一种是数组法
我自己之前一直用的是函数法,因为函数比较直观,但是直到某天在写一道牵涉分块的题目的时候,被卡常了,隔壁队伍的好兄弟告诉我,函数法常数很大,最好要用数组。

分块函数

第一种

len=sqrt(n)
int get(int x)
{
    return (x - 1) / len + 1;
}//函数分块

第二种

blo=sqrt(n)
for(int i=1;i<=n;i++)
{
	bl[i]=(i-1)/blo+1;
}

操作函数

主要就是三个for循环,第一个for循环操作第一个不完整的块,直接暴力遍历,第二个操作最后一个不完整的块,第三个循环是遍历中间完整的块,查询函数也是同理。中间牵涉区间修改的话还要有lazy标记

void modify(int a,int b,int c)
{
    for(int i=a;i<=min(bl[a]*blo,b);i++)
        v[i]+=c,sum[bl[a]]+=c;;
    if(bl[a]!=bl[b])
        for(int i=(bl[b]-1)*blo+1;i<=b;i++)
            v[i]+=c,sum[bl[b]]+=c;
    for(int i=bl[a]+1;i<=bl[b]-1;i++)
        atag[i]+=c;
}

例题

数列分块入门一

#include<map>
#include<set>
#include<cmath>
#include<stack>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define mod 998244353
#define pi acos(-1)
#define inf 0x7fffffff
#define ll long long
using namespace std;
ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,blo;
int bl[50005];
ll v[50005],atag[50005],sum[50005];
void add(int a,int b,int c)
{
    for(int i=a;i<=min(bl[a]*blo,b);i++)
        v[i]+=c,sum[bl[a]]+=c;;
    if(bl[a]!=bl[b])
        for(int i=(bl[b]-1)*blo+1;i<=b;i++)
            v[i]+=c,sum[bl[b]]+=c;
    for(int i=bl[a]+1;i<=bl[b]-1;i++)
        atag[i]+=c;
}
ll query(int a,int b)
{
    ll ans=0;
    for(int i=a;i<=min(bl[a]*blo,b);i++)
        ans+=v[i]+atag[bl[a]];
    if(bl[a]!=bl[b])
        for(int i=(bl[b]-1)*blo+1;i<=b;i++)
            ans+=v[i]+atag[bl[b]];
    for(int i=bl[a]+1;i<=bl[b]-1;i++)
        ans+=sum[i]+blo*atag[i];
    return ans;
}
int main()
{
    n=read();blo=sqrt(n);
    for(int i=1;i<=n;i++)v[i]=read();
    for(int i=1;i<=n;i++)
    {
        bl[i]=(i-1)/blo+1;
        sum[bl[i]]+=v[i];
    }
    for(int i=1;i<=n;i++)
    {
        int f=read(),a=read(),b=read(),c=read();
        if(f==0)add(a,b,c);
        if(f==1)
            printf("%d\n",query(a,b)%(c+1));
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值