ST表
ST表的原理及其实现
ST表类似树状数组,线段树这两种算法,是一种用于解决RMQ(Range Minimum/Maximum Query,即区间最值查询)问题的离线算法
与线段树相比,预处理复杂度同为O(nlogn),查询时间上,ST表为O(1),线段树为O(nlogn)
st表的主体是一个二维数组st[i][j],表示需要查询的数组的从下标i到下标i+2^j - 1的最值,这里以最小值为例
预处理函数:
int a[];//原始输入数组
int st[][];//st表 void init(int n)
{
for (int i = ; i < n; i++)
st[i][] = a[i]; for (int j = ; ( << j) <= n; j++)
{
for (int i = ; i + ( << j) - < n; i++)
st[i][j] = min(st[i][j - ],st[i + ( << (j - ))][j - ]);
}
}
这里首先把从0~n-1的2^0部分进行覆盖,再往下继承
继承这里也很好理解,我们以一个长度为5的数组[5,1,2,3,4]为例
2^0部分覆盖过去自然是5,4,3,2,1
2^1部分的长度为4,从0一直到3,因为从下标为4开始后面只有他自己
st[0][1]是下标为0~1的最小值,自然也就是st[0][0]和st[1][0]的最值
以此往下类推我们可以得出结论:
st[i][j] = min(st[i][j - 1],st[i + 2^(j - 1))][j - 1])
到这里初始化就完成了,注意下标不要越界,如果你对为什么这么处理有困惑的话,请继续看查询
查询函数这里不太好理解
初始化时,每一个状态对应的区间长度都为2^j,由于给出的查询区间长度不一定恰好为2^j,
所以我们要引出一个定理:2^log(a)>a/2 。
因为log(a)表示小于等于a的2的最大几次方。
比如说log(4)=2,log(5)=2,log(6)=2,log(7)=2,log(8)=3,log(9)=3…….
那么我们要查询x到y的最小值。
设len=y-x+1,t=log(len)
根据上面的定理:2^t>len/2
从位置上来说,x+2^t越过了x到y的中间!
因为位置过了一半
所以x到y的最小值可以表示为min(从x往后2^t的最小值,从y往前2^t的最小值)
前面的状态表示为mn[t][x]
设后面(从y往前2^t的最小值)的初始位置是k,
那么k+2^t-1=y,所以k=y-2^t+1
所以后面的状态表示为mn[t][y-2^t+1]
所以x到y的最小值表示为min(mn[t][x],mn[t][y-2^t+1]),所以查询时间复杂度是O(1)
int search(int l, int r)
{
int k = (int)(log((double)(r - l + )) / log(2.0));
return min(st[l][k],st[r - ( << k) + ][k]);
}
示例程序:
#include <iostream>
#include <algorithm> using namespace std; int a[];//原始输入数组
int st[][];//st表 void init(int n)
{
for (int i = ; i < n; i++)
st[i][] = a[i]; for (int i = ; ( << i) <= n; i++)
{
for (int j = ; j + ( << i) - < n; j++)
st[j][i] = min(st[j][i - ],st[j + ( << (i - ))][i - ]);
}
} int search(int l, int r)
{
int k = (int)(log((double)(r - l + )) / log(2.0));
return min(st[l][k],st[r - ( << k) + ][k]);
} int main()
{
int n,m;
while (cin >> n >> m)
{
for (int i = ; i < n; i++)
cin >> a[i]; init(n); while (m--)
{
int l, r;
cin >> l >> r;
cout << search(l,r) << endl;;
}
}
return ;
}
文章内容摘抄于:
ST表的原理及其实现 - 真是啰嗦 - 博客园
ST表算法详解
ST表的更多相关文章
- POJ3693 Maximum repetition substring [后缀数组 ST表]
Maximum repetition substring Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 9458 Acc ...
- 【BZOJ-2006】超级钢琴 ST表 + 堆 (一类经典问题)
2006: [NOI2010]超级钢琴 Time Limit: 20 Sec Memory Limit: 552 MBSubmit: 2473 Solved: 1211[Submit][Statu ...
- 【BZOJ-3956】Count ST表 + 单调栈
3956: Count Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 173 Solved: 99[Submit][Status][Discuss] ...
- 【BZOJ-4569】萌萌哒 ST表 + 并查集
4569: [Scoi2016]萌萌哒 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 459 Solved: 209[Submit][Status] ...
- 【BZOJ-4310】跳蚤 后缀数组 + ST表 + 二分
4310: 跳蚤 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 180 Solved: 83[Submit][Status][Discuss] De ...
- HDU5726 GCD(二分 + ST表)
题目 Source http://acm.hdu.edu.cn/showproblem.php?pid=5726 Description Give you a sequence of N(N≤100, ...
- Hdu 5289-Assignment 贪心,ST表
题目: http://acm.hdu.edu.cn/showproblem.php?pid=5289 Assignment Time Limit: 4000/2000 MS (Java/Others) ...
- Bzoj 2006: [NOI2010]超级钢琴 堆,ST表
2006: [NOI2010]超级钢琴 Time Limit: 20 Sec Memory Limit: 552 MBSubmit: 2222 Solved: 1082[Submit][Statu ...
- ST表poj3264
/* ST表多次查询区间最小值 设 g[j][i] 表示从第 i 个数到第 i + 2 ^ j - 1 个数之间的最小值 类似DP的说 ans[i][j]=min (ans[i][mid],ans ...
- COJ 1003 WZJ的数据结构(三)ST表
WZJ的数据结构(三) 难度级别:B: 运行时间限制:3000ms: 运行空间限制:51200KB: 代码长度限制:2000000B 试题描述 请你设计一个数据结构,完成以下功能: 给定一个大小为N的 ...
随机推荐
- c# ?? 和?
static void Main(string[] args) { double? num1 = null; // ? 说明num1可以为空 ...
- 在ASP.NET MVC 项目中 使用 echarts 画统计图
echarts 官方地址:http://echarts.baidu.com/ 一.根据图中的数据怎么从数据库中获取并组装成对应格式: 从数据库中获取对应数据,然后在项目中引用Newtonsoft.Js ...
- jQuery实现简单的图片淡入淡出效果
整体思路: 1.实现页面布局,设置css样式 2.用jQuery获取需要用到的变量 3.用jQuery为两个按钮绑定事件 一.页面布局: <div class="d1"> ...
- vue 组件通信
组件 组件之间的数据是单向绑定的. 父组件向子组件通信 是通过子组件定义的props属性来实现.通过props定义变量与变量类型和验证方式. props简化定义 在简化定义中,变量是以数组的方式定义. ...
- splay详解(一)
前言 Spaly是基于二叉查找树实现的, 什么是二叉查找树呢?就是一棵树呗:joy: ,但是这棵树满足性质—一个节点的左孩子一定比它小,右孩子一定比它大 比如说 这就是一棵最基本二叉查找树 对于每次插 ...
- SuperMap -WebGL 实现地球的背景透明并显示自定义图片
实现效果如图: 实现代码如下: <!DOCTYPE html> <html lang="en"> <head> <meta charset ...
- 查询SQLSERVER执行过的SQL记录(历史查询记录)
有的时候,需要知道近段时间SQLSERVER执行了什么语句,可以用下面的方法: SELECT TOP 1000 QS.creation_time, SUBSTRING(ST.text, (QS.sta ...
- 如何把activity当成dialog
在工作中有时候需要把activity当成dialog使用,其实做法挺简单的. 1.设置activity的style <style name="DialogTheme" par ...
- gitbook 入门教程之实用插件(新增3个插件)
插件没有什么逻辑顺序,大家可以按照目录大纲直接定位到感兴趣的插件部分阅读即可. 更多插件正在陆续更新中,敬请期待... 最新更新插件 tbfed-pagefooter 版权页脚插件 gitalk 评论 ...
- python读取txt文件最后一行(文件大+文件小)
txt文件小 #coding:utf-8 ''' fname为所读xx.txt文件 输出为:文件第一行和最后一行 ''' fname = 'test.txt' with open(fname, 'r' ...