分块简单入门

分块简单入门

按

树状数组虽超级快，但是很僵硬，不灵活；线段树虽快，但是却不直观，代码量大；所以，速度较慢但直观形象、代码量小的分块大法不实为线段树的替代品。

网络上关于分块的教程不知道为什么很少，虽然早有hzwer大神的分块九讲，但是还是少了入门级详解教程。此篇将分为三个阶段，保证初学者在有意识地阅读后基本掌握分块。

1.简单的入门题

问题引入

给定长度为N(N<=1e5)的正数数列A，然后输入Q(Q<=1e5)行操作命令，指令形如Q l r，表示统计数列中第l~r个数中的最大值

思路

首先想到的简单暴力遍历在这个数量级是会超时，这道题我们可以用分块来做。

先将长度为n，从1开始标号的数组a分为 sqrt(n)块（此时时间复杂度最小）维护，那么易得以下推论

  • a[i]处于(i-1)/sqrt(n)+1块
  • 第i块的左端点为(i-1)*sqrt(n)+1，右端点为min(i*sqrt(n), n)

以上推论便于理解后续实现，可以自己动手模拟验证一下

我们可以先预处理出每一块的最大值并存入数组m[]中，于是对于命令Q l r，我们可以将查询区间[l,r]分为两个部分——所有的整块（整区间）部分和开头和结尾不足一整区间的部分。然后对于整区间我们直接取已维护好了的区间最大值m[i]进行比较，而对于开头和结尾不足一整区间的部分，我们朴素暴力遍历，最终合并答案，得解。

这便是分块算法最基本的思路，概括起来就是大段维护，局部朴素

实现

说明：

blo为区间大小，pos[i]表示a[i]元素位于第pos[i]块，其他变量名同上

例码：

#include <iostream>
#include <cmath>
#define MAXN 100001
#define MAX(a,b) ((a)>(b)?(a):(b))//宏
#define MIN(a,b) ((a)<(b)?(a):(b))
using namespace std;
int blo,n,q,a[MAXN],m[MAXN],pos[MAXN];
int query(int l, int r){
    int ans=0;
    for(register int i=l;i<=MIN(r, pos[l]*blo);i++)//统计左区间（或者查询区间）
        ans=MAX(ans, a[i]);
    for(register int i=pos[l]+1;i<=pos[r]-1;i++)//统计整块区间
        ans=MAX(ans, m[i]);
    if(pos[l]!=pos[r])//如果存在右区间才遍历，防止重复计算
        for(register int i=(pos[r]-1)*blo+1;i<=r;i++)//统计右区间
            ans=MAX(ans, a[i]);
    return ans;
}
int main(){
    ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);//输入输出加速
    cin>>n;
    blo=sqrt(n);
    for(int i=1;i<=n;i++)
        cin>>a[i],pos[i]=(i-1)/blo+1,m[pos[i]]=MAX(m[pos[i]], a[i]);//初始化
    cin>>q;
    char type;
    int l,r;
    while(q--){
        cin>>type>>l>>r;
        cout<<query(l,r)<<endl;
    }
    return 0;
}

2.较难的区间修改

问题引入

给定长度为N(N<=1e5)的数列A，然后输入Q(Q<=1e5)行操作命令

  • 第一类指令形如C l r d，表示把数列中第l~r个数都加d
  • 第二类指令形如Q l r，表示统计数列中第l~r个数的和

思路

同样如上，将长度为n，从1开始标号的数组a分为 sqrt(n)块维护

然后对于操作命令C l r d，[l,r]开头和结尾不足一整区间的部分暴力操作并维护区间之和sum[]；而[l,r]包含的所有整区间则不用暴力操作，而是使用一个add[]数组（类似于线段树的lazy[]延迟标记）记录操作值。

对于命令Q l r也就同理，不足一整区间的部分暴力相加并再加上add[]标记乘以当前区间的个数（因为为了分块算法的效率，操作整块区间时并没有真的修改区间内元素所有值，而是直接用延迟标记add[]暂时记录下来，所以当后面再单独操作取值时，就必须下放延迟标记）；而整块部分则直接sum[]相加即可

总体时间复杂度为

$$

O((N+Q)*\sqrt[]{N})

$$

实现

例码：

#include <iostream>
#include <algorithm>
#include <cmath>
#define MAXN 100010
using namespace std;
int n,q,a[MAXN],blo,sum[MAXN],add[MAXN],pos[MAXN];
void change(int l, int r, int d){
    for(register int i=l;i<=min(pos[l]*blo,r);i++)
        a[i]+=d,sum[pos[l]]+=d;
    for(register int i=pos[l]+1;i<=pos[r]-1;i++)
            add[i]+=d;
    if(pos[l]!=pos[r])//防止重复计算
        for(register int i=(pos[r]-1)*blo+1;i<=r;i++)
            a[i]+=d,sum[pos[r]]+=d;
}
int query(int l, int r){
    int res=0;
    for(register int i=l;i<=min(pos[l]*blo,r);i++)
        res+=a[i]+add[pos[l]];
    for(register int i=pos[l]+1;i<=pos[r]-1;i++)
        res+=sum[i]+add[i]*blo;
    if(pos[l]!=pos[r])//防止重复计算
        for(register int i=(pos[r]-1)*blo+1;i<=r;i++)
            res+=a[i]+add[pos[r]];
    return res;
}
int main(){
    ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);//输入输出加速
    cin>>n;
    blo=(int)sqrt(n);
    for(int i=1;i<=n;i++)
        cin>>a[i],pos[i]=(i-1)/blo+1,sum[pos[i]]+=a[i];//初始化
    cin>>q;
    char type;
    int l,r,d;
    for(int i=0;i<q;i++){
        cin>>type>>l>>r;
        if(type==‘C‘){
            cin>>d;
            change(l,r,d);
        }else
            cout<<query(l,r)<<endl;
    }
    return 0;
}

扩展

其实本题的操作还可以包含区间乘法，即操作涉及区间加法、区间乘法，n个区间和询问。具体思路就是开两个标记数组，分别维护区间加法和区间乘法。值得注意的是，每一次下放标记时，乘法的运算优先级要高于加法，比如当前块加法标记已为m，那么再进行一个乘n的操作后，加法标记应更新为n*m

3.难题

问题引入

[Violet]蒲公英 洛谷 P4168

思路

在线分块大法，首先离散化一下，然后预处理出m[a][b]数组，表示第a块到第b块中蒲公英出现次数最多的编号为m[a][b]；在查询区间[l,r]时，采用分块思想，整区间直接取出与首尾暴力得到的答案进行合并，得解。

实现

（在洛谷上没开O2就RE，开了O2就AC了，也不知道为什么(；′д｀)ゞ）

// luogu-judger-enable-o2
#include <iostream>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <vector>
#include <map>
#define MAXN 50004
#define MIN(a,b) ((a)<(b)?(a):(b))
using namespace std;
int n,q,a[MAXN],num[MAXN],m[505][505][2],blo,cnt[MAXN];
int t[MAXN];
vector <int> po[MAXN];
int getb(int x){//获取a[x]所在块编号
    return (x-1)/blo+1;
}
int getlp(int b){//获取第i块的左端点
    return (b-1)*blo+1;
}
int getrp(int b){//获取第i块的右端点
    return MIN(b*blo, n);
}
int getnum(int l, int r, int x){//获取区间[l,r]中x出现的次数
    return upper_bound(po[x].begin(),po[x].end(),r)-lower_bound(po[x].begin(),po[x].end(),l);
}
void load(int x){//预处理
    memset(cnt, 0, sizeof(cnt));
    int mx=0, mx_num=0;
    for(int i=x;i<=getb(n);i++){
        for(int j=getlp(i);j<=getrp(i);j++){
            cnt[a[j]]++;
            if(cnt[a[j]]>mx||(cnt[a[j]]==mx&&a[j]<mx_num)){//注意如果有若干种蒲公英出现次数相同，则输出种类编号最小的那个
                mx=cnt[a[j]];
                mx_num=a[j];
            }
        }
        m[x][i][0]=mx_num;
        m[x][i][1]=mx;
    }
}
int query(int l, int r){//查询
    int mx_num=m[getb(l)+1][getb(r)-1][0];
    int mx=m[getb(l)+1][getb(r)-1][1];
    for(register int i=l;i<=MIN(r, getrp(getb(l)));i++){
        int t=getnum(l, r, a[i]);
        if(t>mx||(t==mx&&a[i]<mx_num)){
            mx=t;
            mx_num=a[i];
        }
    }
    if(getb(l)!=getb(r))
        for(register int i=getlp(getb(r));i<=r;i++){
            int t=getnum(l, r, a[i]);
            if(t>mx||(t==mx&&a[i]<mx_num)){
                mx=t;
                mx_num=a[i];
            }
        }
    return num[mx_num];
}
int main(){
    ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
    cin>>n>>q;
    blo=sqrt(n);
    for(int i=1;i<=n;i++){
        cin>>a[i],t[i]=a[i];
    }
    sort(t+1, t+1+n);
    int len=unique(t+1, t+1+n)-(t+1);
    for(int i=1;i<=n;i++){
        int temp=lower_bound(t+1, t+1+len, a[i])-t;//从1开始
        num[temp]=a[i],a[i]=temp;
    }
    for(int i=1;i<=n;i++)
        po[a[i]].push_back(i);

    for(int i=1;i<=getb(n);i++) load(i);
    int x=0;
    while(q--){
        int l,r;
        cin>>l>>r;
        l=(l+x-1)%n+1,r=(r+x-1)%n+1;
        if(l>r) swap(l,r);
        x=query(l, r);
        cout<<x<<endl;
    }
    return 0;
}

参考

原文地址：https://www.cnblogs.com/santiego/p/9452424.html

时间： 2024-10-20 16:27:40

分块简单入门的相关文章

正则表达式简单入门

 正则表达式简单入门    正则表达式在平常编程中有着大量的应用,对于任何一个想学习编程的人来说,正则表达式是一个必须掌握的知识. 废话不多说,下面先对正则表达式做一个简单的入门介绍,在后续的文章中,将会进行详细的介绍.    一.元字符 元字符一共有12个:$ ( ) [ { ? + * . ^ \ | 元字符有特殊的含义,如果要使用其字面值,则必须对其进行转义. 如: \$  \*  \( 等等 二.控制字符或不可打印字符 \a  警报 \e  退出 \f  换页 \n  换行 \r 

程序员，一起玩转GitHub版本控制，超简单入门教程 干货2

本GitHub教程旨在能够帮助大家快速入门学习使用GitHub,进行版本控制.帮助大家摆脱命令行工具,简单快速的使用GitHub. 做全栈攻城狮-写代码也要读书,爱全栈,更爱生活. 更多原创教程请关注头条号.每日更新.也可以添加小编微信:fullstackCourse.一起交流,获取最新全栈教程信息.因为FQ原因,不能下载客户端的同仁,可以关注后回复“GitHub客户端”获取安装软件. 上篇教程:GitHub这么火,程序员你不学学吗? 超简单入门教程 干货 GitHub概念部分出现了一丝纰漏.为

iBatis简单入门教程

iBatis 简介: iBatis 是apache 的一个开源项目,一个O/R Mapping 解决方案,iBatis 最大的特点就是小巧,上手很快.如果不需要太多复杂的功能,iBatis 是能够满足你的要求又足够灵活的最简单的解决方案,现在的iBatis 已经改名为Mybatis 了. 官网为:http://www.mybatis.org/ 搭建iBatis 开发环境: 1 .导入相关的jar 包,ibatis-2.3.0.677.jar .mysql-connector-java-5.1.6

Asp.Net MVC学习总结（一）——Asp.Net MVC简单入门

出处:http://www.cnblogs.com/SeeYouBug/p/6401737.html 一.MVC简单入门 1.1.MVC概念 视图(View) 代表用户交互界面,对于Web应用来说,可以概括为HTML界面,但有可能为XHTML.XML和Applet. 模型(Model) 表示用户对其数据的操作的一个封转.可以分为视图模型(view model)和领域模型(domain models),视图模型就是在视图与控制器之间传输数据的一个封转,而领域模型就是业务逻辑,后台数据模型等的一个集

Java日志系统---Logger之简单入门

Java 中自带的日志系统,今天抽空了解了一点,算是入了门,所以将自己的一些心得记录下来,以备日后查看,有兴趣的朋友,看到此文章,觉得有错误或需要添加的地方,请在下方评论留言,大家可以共同进步,谢谢:) Java中关于日志系统的API,在 java.util.logging 包中,在这个包中,Logger类很重要. Logger类是用来记录 某个级别的日志消息: 级别共分为以下几类,从上倒下,级别依次下降: SEVERE(严重)------级别最高 WARNING(警告) INFO CONFIG

响应式网页设计简单入门（强烈推薦！！！！）

响应式网页设计简单入门 Overview: 构造基本的HTML页面 动态加载样式表 Viewport 字体缩放 侧边栏 导航菜单 图片自适应 其他 总结 说到响应式网页设计(Responsive web design),最近在谷歌加上碰到个奇葩贴子,通过一个原始到无法再简单的网页Motherfucking Website及满屏幕的fuck道出了网页设计的真谛,这孩子不是个激进分子就是个报复社会型的货没错,虽然整篇文章就像是泼妇骂街,但我特么是笑着读完的.. 统计了下全文共用Fuck (包括fuc

Android HttpGet() 请求简单入门实例

HttpClient httpclient = new DefaultHttpClient(); String url = "http://example.com"; List<NameValuePair> params = new ArrayList<NameValuePair>(); params.add( new BasicNameValuePair( "param", "value" ) ); URI uri =

delphi指针简单入门

delphi指针简单入门:         看一个指针用法的例子:     1         var     2             X,   Y:   Integer;       //   X   and   Y   整数类型     3             P:   ^Integer;           //   P   指向整数类型的指针     4         begin     5             X   :=17; //   给   X   赋值     6

windows简单入门（1）

1.1 windows系统的关机 退出所有正在运行的程序,然后点击"关闭计算机"键关机,下面是"关闭计算机"的操作方法:                    图1-2 图1-1 1.2鼠标的操作 ① 鼠标指向一个图标 操作:移动鼠标,将其指针放到屏幕上某一对象上或位置 举例:将鼠标指针从屏幕中的某一处移动到"我的电脑"上,这样就指向了"我的电脑". ② 单击鼠标 操作:将鼠标指向某一对象后,快速按一下鼠标左键,并立即释放,这