RMQ模版
2015-09-13 12:23
232 查看
#include<iostream> #include<cmath> #include<algorithm> using namespace std; #define M 100010 #define MAXN 500 #define MAXM 500 int 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;
相关文章推荐
- LintCode Wood Cut
- 第二周项目3(2) 汉诺塔
- 蘑菇街的IM 网络层
- Socket(一)
- 生活感悟
- Android应用的基本组件介绍
- VEX Context(语言环境)
- 网络协议
- Head First Python 第一章 初始Python
- 2015-09-13
- 穿越两次PIX8.0防火墙并两次静态NAT的FTP测试
- 穿越两次PIX8.0防火墙并两次静态NAT的FTP测试
- C#跨线程调用
- 话说Mac终端下 遇到的千奇百怪的报错及解决方法
- web即时通讯2--基于Spring websocket达到web聊天室
- 从特征描述符到深度学习:计算机视觉发展20年
- C++的函数隐藏、覆盖和重载
- 快速幂取模算法
- 项目管理需要一个在线团队协作平台
- C 双向链表