ST算法
来源:互联网 发布:后端开发是指 数据库 编辑:程序博客网 时间:2024/04/27 23:00
原文链接:http://kmplayer.iteye.com/blog/575725
ST算法(Sparse Table):它是一种动态规划的方法。
以最小值为例。a为所寻找的数组.
用一个二维数组f(i,j)记录区间[i,i+2^j-1](持续2^j个)区间中的最小值。其中f[i,0] = a[i];
所以,对于任意的一组(i,j),f(i,j) = min{f(i,j-1),f(i+2^(j-1),j-1)}来使用动态规划计算出来。
这个算法的高明之处不是在于这个动态规划的建立,而是它的查询:它的查询效率是O(1).
假设我们要求区间[m,n]中a的最小值,找到一个数k使得2^k<n-m+1.
这样,可以把这个区间分成两个部分:[m,m+2^k-1]和[n-2^k+1,n].我们发现,这两个区间是已经初始化好的.
前面的区间是f(m,k),后面的区间是f(n-2^k+1,k).
这样,只要看这两个区间的最小值,就可以知道整个区间的最小值!
#include<iostream>#include<cmath>#include<algorithm>using namespace std;#define M 100010#define MAXN 500#define MAXM 500int dp[M][18];/**一维RMQ ST算法*构造RMQ数组 makermq(int n,int b[]) O(nlog(n))的算法复杂度*dp[i][j] 表示从i到i+2^j -1中最小的一个值(从i开始持续2^j个数)*dp[i][j]=min{dp[i][j-1],dp[i+2^(j-1)][j-1]}*查询RMQ rmq(int s,int v)*将s-v 分成两个2^k的区间*即 k=(int)log2(s-v+1)*查询结果应该为 min(dp[s][k],dp[v-2^k+1][k])*/void makermq(int n,int b[]){ int i,j; for(i=0;i<n;i++) dp[i][0]=b[i]; for(j=1;(1<<j)<=n;j++) for(i=0;i+(1<<j)-1<n;i++) dp[i][j]=min(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);}int rmq(int s,int v){ int k=(int)(log((v-s+1)*1.0)/log(2.0)); return min(dp[s][k],dp[v-(1<<k)+1][k]);}void makeRmqIndex(int n,int b[]) //返回最小值对应的下标{ int i,j; for(i=0;i<n;i++) dp[i][0]=i; for(j=1;(1<<j)<=n;j++) for(i=0;i+(1<<j)-1<n;i++) dp[i][j]=b[dp[i][j-1]] < b[dp[i+(1<<(j-1))][j-1]]? dp[i][j-1]:dp[i+(1<<(j-1))][j-1];}int rmqIndex(int s,int v,int b[]){ int k=(int)(log((v-s+1)*1.0)/log(2.0)); return b[dp[s][k]]<b[dp[v-(1<<k)+1][k]]? dp[s][k]:dp[v-(1<<k)+1][k];}int main(){ int a[]={3,4,5,7,8,9,0,3,4,5}; //返回下标 makeRmqIndex(sizeof(a)/sizeof(a[0]),a); cout<<rmqIndex(0,9,a)<<endl; cout<<rmqIndex(4,9,a)<<endl; //返回最小值 makermq(sizeof(a)/sizeof(a[0]),a); cout<<rmq(0,9)<<endl; cout<<rmq(4,9)<<endl; return 0;}
- ST算法
- ST算法
- ST算法
- st算法
- ST算法
- ST算法
- ST算法
- ST算法
- ST算法
- ST算法
- RMQ的ST算法
- RMQ的ST算法
- rmq-st算法
- RMQ问题ST算法
- Uva11235 RMQ ST算法
- RMQ问题 ST算法
- poj3264(ST算法)
- RMS-ST算法
- ACM计算几何题目推荐
- MFC实现父窗口不透明,子窗口半透明效果
- 机房收费系统总结之2——配置文件(初尝)
- 【转载】ibatis 属性详细分析
- HTTPS是如何保证连接安全:每位Web开发者都应知道的
- ST算法
- 【职场版】《漫漫邮子路(一)--走出你身处的“盒子”——朱斌》
- STM32 使用DMA处理ADC 学习笔记!!
- 24.从零开始学习C语言--二进制全部为0的意义
- 有多少条病狗?
- 短训学习录(三)——体会到上导教班的意义
- 我眼里的帝都和她哺育的孩子们
- DevExpress图表控件ChartControl
- 调用哪个方法?