ST (Sparse Table:稀疏表)算法
1541:【例 1】数列区间最大值
时间限制: 1000 ms 内存限制: 524288 KB
提交数: 600 通过数: 207
【题目描述】
输入一串数字,给你 MM 个询问,每次询问就给你两个数字 X,YX,Y,要求你说出 XX 到 YY 这段区间内的最大数。
【输入】
第一行两个整数 N,MN,M 表示数字的个数和要询问的次数;
接下来一行为 NN 个数;
接下来 MM 行,每行都有两个整数 X,YX,Y。
【输出】
输出共 MM 行,每行输出一个数。
【输入样例】
10 2
3 2 4 5 6 8 1 2 9 7
1 4
3 8
【输出样例】
5
8
【提示】
数据范围与提示:
对于全部数据,1≤N≤105,1≤M≤106,1≤X≤Y≤N1≤N≤105,1≤M≤106,1≤X≤Y≤N。数字不超过 C/C++C/C++ 的 intint 范围。
【来源】http://ybt.ssoier.cn:8088/problem_show.php?pid=1541
解析: 这是指求区间最值的模板问题,这里采用 ST算法来实现,总的时间复杂度O(nlogn+Q)
。
倍增思想:
f[i][j]表示i开始的连续2j 个点的最大值。
则f[i][0]表示i开始连续1个点的最大值即a[i];
f[i][1]表示i开始连续2个点的最大值即a[i]和a[i+1]的最大值;
f[i][2]表示i开始连续4个点的最大值即a[i]~a[i+3]中的最大值;
f[i][3]表示i开始连续8个点的最大值即a[i]~a[i+7]中的最大值;
......
f[i][log(n)/log(2)开始连续n个点的最大值即 a[i]~a[i+n-1];(i+n-1<=n)
f[1][0]=a[1]=3; | f[2][0]=a[2]=2; | f[3][0]=a[3]=4; | f[4][0]=a[4]=5; | f[5][0]=a[5]=6; | f[6][0]=a[6]=8; | f[7][0]=a[7]=1; | f[8][0]=2; | f[9][0]=9; | f[10][0]=7; |
f[1][1]=max(f[1][0],f[2][0])=3 | f[2][1]=max(f[2][0],f[3][0])=4 | f[3][1]=max(f[3][0],f[4][0])=5 | f[4][1]=max(f[4][0],f[5][0])=6 | f[5][1]=max(f[5][0],f[6][0])=8 | f[6][1]=8 | f[7][1]=2 | f[8][1]=9 | f[9][1]=9 | |
f[1][2]=max(f[1][1],f[3][1])=4 | f[2][2]=max(f[2][1],f[4][1])=6 | f[3][2]=max(f[3][1],f[5][1])=8 | f[4][2]=max(f[4][1],f[6][1])=8 | f[5][2]=max(f[5][1],f[7][1])= 8 | f[6][2]=9 | f[7][2]=9 | |||
f[1][3]=max(f[1][2],f[5][2])=8 | f[2][3]=max(f[2][2],f[6][2])=9 | f[3][3]=max(f[3][2],f[7][2])=9 | |||||||
ST:
void st(){ //时间复杂度O(nlogn)
int k=log(n)/log(2);//求深度
for(int i=1;i<=n;i++)f[i][0]=a[i];//初始化
for(int j=1;j<=k;j++)
for(int i=1;i+(1<<j)-1<=n;i++)//长度为2^j的区间[i,i+2^j-1] ,所以当i+(1<<j)-1>n时循环停止
f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
}
RMQ(Range Minimum/Maximum Query)问题 //区间最值查询:
int qwt(int L,int R){//时间复杂度O(1)
int k=log(R-L+1)/log(2);//2k<=R-L+1,但2*2k>R-L+1
return max(f[L][k],f[R-(1<<k)+1][k]);//L+2k-1<=R,R-2k+1>=L,所以区间[L,L+2k-1]和[,R-2k+1,R]是正好相切或者有部分重合,也就是说两个区间是无缝衔接的,这样求两个区间最值的最值就是所求答案。
}
L=1,R=4时 求的是 max(f[1][2],f[1][2]) 即区间【1,4】【1,4】
L=3,R=8时 求的是 max(f[3][2],f[6][2])即区间【3,6】【5,8】,两个区间中5,6是重合的
L=1,R=10时 求的是 max(f[1][3],f[3][10])即区间【1,8】【3,10】,两个区间中3~8是重合的
#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
const int maxn=;
int n,m;
int f[maxn][],a[maxn];
void st(){
int k=log(n)/log();
for(int i=;i<=n;i++)f[i][]=a[i];
for(int j=;j<=k;j++)
for(int i=;i+(<<j)-<=n;i++)//长度为2^j的区间[i,i+2^j-1]
f[i][j]=max(f[i][j-],f[i+(<<(j-))][j-]);
}
int qwt(int L,int R){
int k=log(R-L+)/log();
return max(f[L][k],f[R-(<<k)+][k]);
}
int main(){
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++)
scanf("%d",&a[i]);
st();
for(int i=;i<=m;i++){
int x,y;
scanf("%d%d",&x,&y);
printf("%d\n",qwt(x,y));
}
return ;
}
注意:当输入输出数据的规模达到10的6次方时,就需要用scanf和pritf输入输出。此时用cin和cout是绝对会超时的。
ST (Sparse Table:稀疏表)算法的更多相关文章
- RMQ(ST(Sparse Table))(转载)
1. 概述 RMQ(Range Minimum/Maximum Query),即区间最值查询,是指这样一个问题:对于长度为n的数列A,回答若干询问RMQ(A,i,j)(i,j<=n),返回数列A ...
- 稀疏表(ST / Sparse Table)
RMQ问题: 给定一个序列,每次询问一个区间最小值 / 最大值. 没有修改. //拿区间最大值来举例. memset(ans, -INF, sizeof(ans)); for (int i = 1; ...
- RMQ (Range Minimal Query) 问题 ,稀疏表 ST
RMQ ( 范围最小值查询 ) 问题是一种动态查询问题,它不需要修改元素,但要及时回答出数组 A 在区间 [l, r] 中最小的元素值. RMQ(Range Minimum/Maximum Query ...
- codeforce 359D 二分+ 动态规划(sparse table)
原题链接:http://codeforces.com/problemset/problem/359/D 思路:首先对符合题目的长度(r-l)从0到n-1进行二分查找,对每一个长度进行check,看是否 ...
- 算法学习 - ST表 - 稀疏表 - 解决RMQ问题
2017-08-26 21:44:45 writer:pprp RMQ问题就是区间最大最小值查询问题: 这个SparseTable算法构造一个表,F[i][j] 表示 区间[i, i + 2 ^ j ...
- RMQ ---- ST(Sparse Table)算法
[概述] RMQ(Range Minimum/Maximum Query),即区间最值查询,是指这样一个问题:对于长度为n的数列A,回答若干询问RMQ(A,i,j)(i,j<=n),返 ...
- 基于稀疏表(Sparse Table)的RMQ(区间最值问题)
在RMQ的其他实现方法中,有一种叫做ST的算法比较常见. [构建] dp[i][j]表示的是从i起连续的2j个数xi,xi+1,xi+2,...xi+2j-1( 区间为[i,i+2j-1] )的最值. ...
- 动态规划——稀疏表求解RMQ问题
RMQ (Range Minimum/Maximum Query)问题,即区间最值查询问题,是求解序列中的某一段的最值的问题.如果只需要询问一次,那遍历枚举(复杂度O(n))就是最方便且高效的方法,但 ...
- 一维二维Sparse Table
写在前面: 记录了个人的学习过程,同时方便复习 Sparse Table 有些情况,需要反复读取某个指定范围内的值而不需要修改 逐个判断区间内的每个值显然太浪费时间 我们希望用空间换取时间 ST表就是 ...
随机推荐
- 【OF框架】搭建标准工作环境
前言 统一工作环境,统一工具集,是沟通效率的基础.如同语言一样,使用不同语言的人,需要花更多的精力去理解语言,然后才是理解语言的内容,而使用相同语言的人,对话过程直接进入到内容.对于语言不统一增加的效 ...
- 读取xml时,报错:xml.etree.ElementTree.ParseError: no element found: line 20, column 9
读取xml时,出现报错:xml.etree.ElementTree.ParseError: no element found: line 20, column 9 原因是xml文件格式有问题,可以检查 ...
- Scala环境安装步骤
1.scala解释器本地安装 2.IDEA安装 3.安装IDEA的scala插件 4.创建maven项目 5.安装ScalaSDK
- 牛客练习赛50 D tokitsukaze and Event (最短路,思维)
牛客练习赛50 D tokitsukaze and Event 链接:https://ac.nowcoder.com/acm/contest/1080/D来源:牛客网 时间限制:C/C++ 1秒,其他 ...
- MyBatis3_[tp_41-42-43]-_动态sql_trim_自定义字符串截取_choose分支选择_update的set与if-trim 结合的动态更新
笔记要点出错分析与总结 /** 笔记: * 查询的时候,如果某些条件,没带可能SQL拼装会有问题; * 1.-->给where 后面加上 1=1, 以后的条件都and XXX * 2. < ...
- 学习php doctrine
了解symfony3.3.PHP FIG.Doctrine: 了解angular2.material2:. 熟悉git:了解开源项目:openstack docker ceph等: NoSQL(HBa ...
- SHELL编程基础01
首先shell是在linux下运行的一种环境,它是以shell脚本来运行的,学会了它基本可以解决任何问题,也可以用shell脚本开发. 和java,python的相比,其弱类型的语言没有那么复杂的结构 ...
- 大数据之路week06--day07(完全分布式Hadoop的搭建)
前提工作: 克隆2台虚拟机完成后:新的2台虚拟机,请务必依次修改3台虚拟机的ip地址和主机名称[建议三台主机名称依次叫做:master.node1.node2 ] 上一篇博客 (三台虚拟机都要开机) ...
- java TCP 通信:服务端与客服端
1.首先先来看下基于TCP协议Socket服务端和客户端的通信模型: Socket通信步骤:(简单分为4步) 1.建立服务端ServerSocket和客户端Socket 2.打开连接到Socket的输 ...
- java中的AIO
AIO(异步非阻塞)AIO采用了Proactor模式,AIO与NIO的不同之处在于当AIO在进行读写操作时,不用先等通知,可直接调用相应的read/write方法,这两种方法均为异步的,对于读操作而言 ...