RMQ问题ST算法//洛谷P3865

第一行——咕咕咕（毕竟本来说是要学线段树的。

 

RMQ问题：区间最值查询。即：RMQ(a,i,j)，返回数列a[i]到a[j]中最小or最大元素的下标。询问一次时，跑一次for循环即可，但询问次数多时，就会超时，因此引出ST算法，在nlogn的时间复杂度里预处理，然后O(1)询问。

 

ST算法：

1.用dp[100005][100005]来表示从i起始的连续2^j个数的min或max的值，即dp[1][2]为1到4区间内4个数的max或min值。

2.用dp[i][j] = max（dp[i][j-1]，dp[i+(1<<j)-1][j]）进行维护，其实计算一下就会发现，这两个刚好是i到i+2^j之间的左半和右半区间，完美覆盖。

3.dp[i][0]为i本身，所以初始化为a[i].

4.以L为起点向右查询，以R为终点向左查询，两边比较得最大值是因为cd<=log2(n)，单方向查询不一定全覆盖区间。

 

本代码以洛谷P3865为例（有趣的是，数组我开20RE了，开2000MLE，开200AC美滋滋（手动再见））

需要注意的是，位运算的优先级是低于+or-的（我因为这个也wa过（没想到吧.jpg

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#define ll long long

using namespace std;

int a[100005];
int dp[100010][200];
int n,m;

void init()
{
    for(int i = 1; i <= n; ++i)
        dp[i][0] = a[i];
    for(int j = 1; j <= log2(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]);

}

int query(int l,int r)
{
    int cd = log2(r-l+1);
    return max(dp[l][cd],dp[r-(1<<cd)+1][cd]);
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i = 1; i <= n; ++i)
        scanf("%d",&a[i]);
    init();
    while(m--)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        int ans = query(x,y);
        printf("%d\n",ans);
    }
    return 0;
}

 

好了明天我再更新其他习题（咕咕咕）

欢迎指出错误qwq