RMQ，对于区间最值的一些想法

东辛木

已于 2023-07-10 17:43:47 修改

阅读量60

点赞数

分类专栏：学习日记文章标签：算法动态规划

于 2023-07-10 17:43:45 首次发布

本文链接：https://blog.csdn.net/m0_72610805/article/details/131641153

版权

学习日记专栏收录该内容

16 篇文章 0 订阅

订阅专栏

最近集训，啥都不会，心里有点烦，依稀记得初见st表板子题（p3865)的时候,直接朴素搜索，不出意外的RE了很多点，心灰意冷遂弃之。做任何事，最忌半途而废。且亡羊补牢，从ST表，进一步了解RMQ问题吧。

ST表

oi-wiki:ST表

本质是DP，是利用倍增思想来缩短时间的。而倍增就体现在他数组的定义中：对于f[i][j]，指的是在序列的第i项，向后 $2^j$ 个元素所包含序列间的最大值。

例如数列3 2 4 5 6 8 1 2 9 7 ,st[1，0]表示第1个数起，长度为2^0 =1的最大值，其实就是3这个数。st[1，2]=5，st[1，3]=8，st[2，0]=2，st[2，1]=4……从这里可以看出f[i,0]其实就等于a[i]。这样，DP的状态、初值都已经有了，剩下的就是状态转移方程。

把f[i，j]平均分成两段（因为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[i，j]就是这两段的最大值中的最大值。再解释一下，从第i个位置开始，往后数2^(j-1)的区间内的最值为m1，则st[i][j-1]表示的就是m1；那第二段区间的起点即为2^(j-1)+1（注意！要+1！）依旧往后数2^(j-1）个数，该区间内最值为m2，则st[i][j]的值为max(m1,m2)。于是我们得到了DP方程st[i,j]=max(f[i，j-1],f[i+2^(j-1)，j-1]）.

我们令一个整型k = (log(r-l+1)/log(2))

动态l = r -2^k + 1

计算最值，把区间[l，r]分成两个长度为2^k的区间（保证有f[i，j]对应）。

ans = max(F[l，k]，F[r-2^k+1，k])
例题，洛谷P3865,下为思路，题目数据点有点怪，后五个T了，复制别人题解过去也是TLE

#include <bits/stdc++.h>
typedef long long ll;
const int N = 1e6;
using namespace std;
int f[100100][40],a,x,LC,n,m,p,len,l,r;
int main()
{
    cin>>n>>m;
    for(int i=1;i<=n;i++)
    {	
        cin>>a;
        f[i][0]=a;
    }
    LC=(int)(log(n)/log(2));

    for(int j=1;j<=LC;j++)
        for (int i=1;i<=n-(1<<j)+1;i++)
    		f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
    for(int i=1;i<=m;i++) 
    {
        cin>>l>>r;
  		p=(int)(log(r-l+1)/log(2));
		cout<<max(f[l][p],f[r-(1<<p)+1][p])<<"\n";
    }
    return 0;
}