st表-RMQ问题解法

RMQ

RMQ 问题是指:对于长度为 n 的数列 A,回答若干询问 RMQ (A , i , j ) ( i , j ≤ n),返回数列A中下标在 i , j 里的最小(大)值,也就是说,RMQ问题是指求区间最值的问题。

二:st算法

ST算法是动态规划以及倍增的思想相结合。只适用于静态区间求最值,如果是动态的,需要用线段树

2.1 基本思想

首先我们定义f[i][j]表示以下标i为起点,连续2^j个数中的最大值。eg:f[2][2]表示:从a[2]-a[5]的最大值。f[i][0]也就表示为a[i]了。

2.2 通过倍增的方法获取值

如我们想要取得a[2]-a[14]中间的最大值的话,如果我们能够二分的获得里面的值就好了比如a[2]-a[9],a[9]-a[14],或者说a[2]-a[9],a[7]-a[14]的话就好了。这样我们取最大值就可以得到a[2]-a[14]的最大值了。
居然f[i][j]中的j代表2^j,又因为我们需要二分,很自然的想到了对数运算了。
st表-RMQ问题解法

那么我们可以得到查询的时候的方法:

//查询a[l]-a[r]的最大值
int query(int l,int r){
    int k = log2(r-l+1);
    return max(f[l][k],f[r-(1<<k)+1][k]);
}

2.3预处理

为了让我们能够实现2.2中的算法,我们需要对f[i][j]进行处理,使其支持区间查询操作。

//j的大小可以根据实际数据的范围来选择,一般数据为1e5,那么21位的完全够用
for(int j=1;j<=21;j++){
	//这里的n代表a[]的总长度,因此当前操作区间长度也就是i+(1<<j)-1
	//可以参考2.2中为什么+1,这就是为什么-1了。
        for(int i=1;i+(1<<j)-1<<n;i++)
            f[i][j] = max(f[i][j-1],f[i+1<<(j-1)][j-1]);
}

模板总结

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
using namespace std;

const int N = 1e5+10;
int f[N][21];
inline int read()
{
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}

int query(int l,int r){
    int k = log2(r-l+1);
    return max(f[l][k],f[r-(1<<k)+1][k]);
}


int main(){
    int n,m;
    n=read();
    m=read();
    for(int i=1;i<=n;i++){
        f[i][0]=read();
    }

    for(int j=1;j<=21;j++){
        for(int i=1;i+(1<<j)-1<<n;i++)
            f[i][j] = max(f[i][j-1],f[i+1<<(j-1)][j-1]);
    }


    for(int i=1;i<=m;i++){
        int l=read();
        int r=read();
        printf("%d\n",query(l,r));
    }


    return 0;
}
上一篇:矿视界译文:Chia网络是如何运作的?


下一篇:chia/xch rpc形式充提生成地址