RMQ(Range Maximum/Minimum Query)

taoyiwei17_HNCS

已于 2024-08-14 16:17:16 修改

阅读量434

点赞数 6

文章标签：算法

于 2024-08-14 16:16:27 首次发布

本文链接：https://blog.csdn.net/taoyiwei17_QFOI/article/details/141194062

版权

RMQ(Range Maximum/Minimum Query)

定义

指的是在一个序列中多次进行静态区间求最值。

算法原理

倍增思想。

实现

分为两部分：预处理和询问（S 组算法里最长的预处理）。

first part：预处理

维护 dp_{i,j}dpi,j 表示以下标 ii 为起点，跨度为 2^j2j 的最大值。

状态转移：dp_{i,j}=\max(dp_{i,j-1},dp_{i+2^{j-1},j-1})dpi,j=max(dpi,j−1,dpi+2j−1,j−1)

初始状态：dp_{i,0}=a_idpi,0=ai。

维护 lg_ilgi 表示 ii 取以 22 为底的对数向下取整。lg_i=-1lgi=−1，lg_i=lg_{i>>1}+1lgi=lgi>>1+1。

second part：询问

对于 [l,r][l,r] 的询问，一定能找到 22 的幂超过 r-l+1r−l+1 的一半。

求覆盖区间一半以上的幂，p=lg_{r-l+1}p=lgr−l+1。

\max(dp_{l,p},dp_{r-2^p+1,p})max(dpl,p,dpr−2p+1,p) 即为答案。

代码

#include<bits/stdc++.h>
#define int long long
using namespace std;
int dp[100005][20], a[100005], n, m, lg[100005];
void init()
{
	lg[0]=-1;
	for(int i=1;i<=n;i++)
	{
		dp[i][0]=a[i];
		lg[i]=lg[i>>1]+1;
	}
	for(int j=1;(1<<j)<=n;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]);
		}
	}
}
signed main()
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	{
		cin>>a[i];
	}
	init();
	while(m--)
	{
		int l, r;
		cin>>l>>r;
		int mi=lg[r-l+1];
		cout<<max(dp[l][mi],dp[r-(1<<mi)+1][mi])<<"\n";
	}
}