RMQ小结

2024-05-03 19:32
文章标签 小结 rmq

本文主要是介绍RMQ小结,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

                                  RMQ小结

区间求值得算法主要有三种:

一、处理O(N),查询O(N*Q)   (朴素)

二、处理O(log N),查询(QlogN) (线段树)

三、处理O(nlogn),查询(1)    (RMQ)

而我们主要来讲一下,O(nlogn)-O(1)的RMQ算法。而RMQ算法的实现又有多种算法,我就选了一种性价比最高的讲解。就是代码容易,时间也可以的算法。因为正宗的RMQ的时间是O(N)-O(1),我主要说的是ST算法O(NlogN)-O(1)的算法。

RMQ的本质思想是动态规划。度娘如是解释:

ST算法(Sparse Table),以求最大值为例,设d[i,j]表示[i,i+2^j-1]这个区间内的最大值,那么在询问到[a,b]区间的最大值时答案就是max(d[a,k], d[b-2^k+1,k]),其中k是满足2^k<=b-a+1(即长度)的最大的k,即k=[ln(b-a+1)/ln(2)]。
d的求法可以用动态规划,d[i, j]=max(d[i, j-1],d[i+2^(j-1), j-1])。

看一下ST算法是怎么实现的(以最大值为例):
   首先是预处理,用一个DP解决。设a是要求 区间 最值的 数列 ,f[i,j]表示从第i个数起连续2^j个数中的最大值。例如数列3 2 4 5 6 8 1 2 9 7 ,f[1,0]表示第1个数起,长度为2^0=1的最大值,其实就是3这个数。f[1,2]=5,f[1,3]=8,f[2,0]=2,f[2,1]=4……从这里可以看出f[i,0]其实就等于a[i]。这样,DP的状态、初值都已经有了,剩下的就是 状态转移方程 。我们把f[i,j](j≥1)平均分成两段(因为j≥1时,f[i,j]一定是偶数个数字),从i到i+2^(j-1)-1为一段,i+2^(j-1)到i+2^j-1为一段(长度都为2^(j-1))。用上例说明,当i=1,j=3时就是3,2,4,5 和6,8,1,2这两段。f就是这两段的最大值中的最大值。于是我们得到了动规方程F[i,j]=max(F[i,j-1],F[i+2^(j-1),j-1])。
接下来是得出最值,也许你想不到计算出f有什么用处,一般要想计算max还是要O(logn),甚至O(n)。但有一个很好的办法,做到了O(1)。还是分开来。如在上例中我们要求区间[2,8]的最大值,就要把它分成[2,5]和[5,8]两个区间,因为这两个区间的最大值我们可以直接由f[2,2]和f[5,2]得到。扩展到一般情况,就是把区间[l,r]分成两个长度为2^n的区间(保证有f对应)。直接给出表达式:
k:=trunc(ln(r-l+1)/ln(2));
ans:=max(F[l,k],F[r-2^k+1,k]);
这样就计算了从l开始,长度为2^k的区间和从r-2^k+1开始长度为2^k的区间的最大值(表达式比较繁琐,细节问题如加1减1需要仔细考虑),二者中的较大者就是整个区间[l,r]上的最大值。

#include <iostream>
#include <cstdio>
#include <cmath>
#include<algorithm>
#define MN 50005
using namespace std;
int mi[MN][17],mx[MN][17],w[MN];
int n,q;
void rmqinit()
{int i,j,m;for(i=1;i<=n;i++){mi[i][0]=mx[i][0]=w[i];}int m = int(trunc(log2(n)));for(i=1;i<=m;i++){for(j=n;j>=1;j--){mx[j][i]=mx[j][i-1];if(j+(1<<(i-1))<=n)mx[j][i]=max(mx[j][i],mx[j+(1<<(i-1))][i-1]);mi[j][i]=mi[j][i-1];if(j+(1<<(i-1)<=n))mi[j][i]=min(mi[j][i],mi[j+(1<<(i-1))][i-1]);}}
}int rmqmin(int l,int r)
{int m=int(trunc(log2(r-l+1)));return min(mi[l][m],mi[r-(1<<m)+1][m]);
}int rmqmax(int l,int r)
{int m=int(trunc(log2(r-l+1)));return max(mx[l][m],mx[r-(1<<m)+1][m]);
}int main()
{cin>>n>>q;for(int i=1;i<=n;i++)scanf("%d",&w[i]);//cin>>w[i];rmqinit();int l,r;for(int i=1;i<=q;i++){scanf("%d%d",&l,&r);//cin>>l>>r;printf("%d %d\n",rmqmax(l,r),rmqmin(l,r));//cout<<rmqmax(l,r)<<" "<<rmqmin(l,r)<<endl;}while(cin>>n)return 0;
}




练练手.HDU3486
题目链接: Click Here~

题目分析:
题目重述,题目说给出一个老板要面试N个人,且这N个人,每个人都有自己的能力值。但是他没有那么多的时间。因此,他决定将这N个人按原来报道的顺序将他们分成m组。而剩下的人由于来的晚,所以,就直接不给他们面试的机会了。题目问,分别从m组中挑选出N/m个人。且boss希望员工在m组中每组的最高能力值的总和不小于k。而要你求出至少要分成几组?

算法分析:
朴素的算法,求出组数是从1~N的遍历,在求出这1~N组中每组的人的最大能力值。这明显超时。但是哪里可以优化呢?我们发现程序中只有两处核心;一,求组数。二,求每组中人的最大能力值。只要解决了这个两个问题,当然就可以解决了超时的问题。
如何,优化组数呢?我们可以想最好的情况是没跟人的能力值都一样的,那么此时就是组数的下界。那么,我们就可以求出N个人中最大能力值(maxv)是多少。而最小组数就是k/maxv了(想想为什么?)而上界就是N组.‘
如何,优化每组中人的最大能力值呢?显然就是我们开篇讲的RMQ.

然后呢,这道题的数据有问题。可以用二分,但是本质上着N个人的能力值并不存在任何关系,因此根本就不能用二分。还有题目描述也很恶心。到底是保证组数最小呢?还是人数最小呢?这要自己猜吗?Are you kidding me?

#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;const int N = 2e5 + 5;
int n,dp[N][20];void RMQ_init()
{int m = int(trunc(log2(n)));
//    int m = (int)(log10(n*1.0)/log10(2.0));for(int j = 1;j <= m;++j)for(int i = 1;i+(1<<j)-1 <= n;i++)dp[i][j] = max(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
}
int RMQ(int L,int R)
{
//    int k = (int)(log10(R-L+1)*1.0/log10(2.0));int k = int(trunc(log2(R-L+1)));return max(dp[L][k],dp[R-(1<<k)+1][k]);
}
int main()
{int k;while(scanf("%d%d",&n,&k),(n>=0&&k>=0)){int sum = 0,maxv = 0;for(int i = 1;i <= n;++i){scanf("%d",&dp[i][0]);sum += dp[i][0];maxv = max(maxv,dp[i][0]);}if(sum <= k){puts("-1");continue;}RMQ_init();int ans = n,st = k/maxv;if(!st) st++; //注意出现除0的情况---->Runtime Error!!!!!!!!!!!!!!for(int i = st;i < n;++i)  //枚举组数{int seg,s = 0;seg = n/i;            //可以分成几个区间for(int j = 1;j <= i;j++) //在i个区间内,各个区间的起点和终点{s += RMQ((j-1)*seg+1,j*seg);if(s > k)break;}if(s > k){ans = i;break;}}printf("%d\n",ans);}return 0;
}





                         AMagic Lamp



题目链接:Click Here~

题意解析:

    题目要求,给出一个字符串。让你在删除m个字符之后,求剩下的数组成的最小数是多少。(同理,求剩下的数组成的最大数也是同样的解法)

思路解析:

显然,要删除m个数则原字符串比剩下n-m个(原字符串长度为n.则,我们知道当两个数的位数相同的时候,就要比较他的首位。

比如,19992000.因此,我们要尽量的使首位变小。而这只要在一个区间中找出最小的那个就好了。而问题有出现了,要怎么确定去区间呢?显然,我们在找第一位的时候,后面肯定要至少剩下n-m-1个数。且我们假设我们找到的第一位索引为i1。则第二位的区间可以很容易的确定是从i+1为起点后面至少剩下n-m-2个的区间内。同理,后面的都是相同的道理。

    而我们知道区间求极值,我们又可以使用RMQ求得最值。而当运用RMQ的时候我们又可以从逆向思维来考虑。当我们找第一位的时候既然后面要剩下n-m-1个,那不是说明第一位的最值一定是在[0,m]内,第二个在[i1,m+1],ect?这是显然的。所以,我们可以轻松的运用RMQ求出每个区间的最值是多少。



#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;const int N = 1e3 + 5;
int n,pos[N][20];
char str[N];
void RMQ_init()
{for(int i = 0;i <= n;++i) pos[i][0] = i;int m = int(trunc(log2(n)));for(int j = 1;j <= m;++j)for(int i = 0;i +(1<<j)-1 <= n;++i){int s1 = pos[i][j-1];int s2 = pos[i+(1<<(j-1))][j-1];pos[i][j] = str[s1]<=str[s2]?s1:s2;}
}
int RMQ(int L,int R)
{int m = int(trunc(log2(R-L+1)));int s1 = pos[L][m];int s2 = pos[R-(1<<m)+1][m];return str[s1] <= str[s2]?s1:s2;
}
int main()
{int m;while(~scanf("%s%d",str,&m)){n = strlen(str); // printf("str ---> %s   n -- > %d\n",str,n);RMQ_init();
//        char res[N];bool first = true;int index,idx = 0,j = m;
//        memset(res,'\0',sizeof(res));for(int i = 1;i <= n-m;++i){index = RMQ(idx,j);j++;if(first&&str[index]!='0'){first = false;putchar(str[index]);}else if(!first){putchar(str[index]);}idx = index+1;}if(first) putchar('0');putchar('\n');}return 0;
}

 

Max Sum of Max-K-sub-sequence


     也可以运用单调队列。Click Here~




                                




这篇关于RMQ小结的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/957475

相关文章

分布式系统的个人理解小结

分布式系统:分的微小服务,以小而独立的业务为单位,形成子系统。 然后分布式系统中需要有统一的调用,形成大的聚合服务。 同时,微服务群,需要有交流(通讯,注册中心,同步,异步),有管理(监控,调度)。 对外服务,需要有控制的对外开发,安全网关。

Linux环境配置中问题小结

在Linux环境配置中,遇到问题首先猜测: 1、是否是权限问题; 2、软连接是否配置;

long long,_int64使用小结

前言:   在16位环境下,int/unsigned int 占16位,long/unsigned long占32位   在32位环境下,int占32位,unsigned int占16位,long/unsigned long占32位 何时需要使用:   long 和 int 范围是[-2^31,2^31),即-2147483648~2147483647,而unsigned范围是[0,2^32),

密码学读书笔记小结

密码学是保证消息的私密性和完整性以及消息认证的基础。加密算法的选择和密钥的管理是安全机制的效率、性能和可用性的关键。 公钥加密算法: 分发密钥比较容易,但是对大数据量的加密性能较差密钥加密算法: 更适合大批的加密任务混合型加密协议: 例如TLS,先用公钥加密建立一个安全通道,然后使用通道交换密钥,并将此密钥用于后续数据交换。 对分布式系统攻击的分类: 窃听: 未经授权获得消息副本伪装: 在未

Android 源码中jni项目 加载so目录小结

Android 源码中jni项目 加载so目录小结 文章目录 Android 源码中jni项目 加载so目录小结一、前言二、so目录验证测试1、jni so文件错误报错(1)报错1 - 未找到so文件:(2)报错2 - so文件中未找到native方法: 2、验证的几种情况(1)apk下面的 lib/arm64/ 放置正确的so文件(2)apk下面的 lib/arm64/ 放置错误的so文

maven打包成可执行的jar,以及读取配置文件问题小结

文章来源 https://blog.csdn.net/chasonsp/article/details/88852353 折腾的几天,使用maven打包后发现了问题,首先是打包的配置文件读取问题,使用getResource().getPah()会发现在访问jar包的文件时,路径里会有感叹号(杠杠滴~~)是这样的 …jar!.. 经过不断的查找资料及反复验证后,终于找到了可行的方法:

【硬刚大数据】Flink在实时在实时计算平台和实时数仓中的企业级应用小结

欢迎关注博客主页:https://blog.csdn.net/u013411339 欢迎点赞、收藏、留言 ,欢迎留言交流!本文由【王知无】原创,首发于 CSDN博客!本文首发CSDN论坛,未经过官方和本人允许,严禁转载! 本文是对《【硬刚大数据之学习路线篇】从零到大数据专家的学习指南(全面升级版)》的面试部分补充。 大数据领域自 2010 年开始,以 Hadoop、Hive 为代

Kafka常见问题学习路径源码阅读小结 | 写在Kafka3.0发布之际

严格来说,这篇文章也不是今天写的。是之前断断续续写在了几篇文章中。 2021年9月21日,随着Kafka3.0的发布,Kafka在「分布式流处理平台」这个目标上的努力进一步得到加强!Kafka不满足于「消息引擎」的定位,正式基于这样的定位,Kafka 社区于 0.10.0.0 版本正式推出了流处理组件 Kafka Streams,也正是从这个版本开始,Kafka 正式"变身"为分布式的流处理平台

C++中继承及virtual小结

一、继承基础知识 C++中的继承 1.1继承的基本概念 类与类之间的关系 has-A,包含关系,用以描述一个类由多个“部件类”构成,实现has-A关系用类的成员属性表示,即一个类的成员属性是另一个已经定义好的类。 use-A,一个类使用另一个类,通过类之间的成员函数相互联系,定义友元或者通过传递参数的方式来实现。(和组合不同) is-A,即继承关系,关系具有传递性。 继承的特点 子类拥有

【HDU】5436 Transmigration tree【树链剖分+dp+rmq】

题目链接:【HDU】5436 Transmigration tree #pragma comment(linker, "/STACK:16777216")#include <stdio.h>#include <string.h>#include <vector>#include <algorithm>using namespace std ;typedef long long LL ;