[st算法]对于RMQst算法的一些理解

RMQ算法简介

RMQ（Range Minimum/Maximum Query），即区间最值查询，是指这样一个问题：对于长度为n的数列A，回答若干询问RMQ（A,i,j）(i,j<=n)，返回数列A中下标在i，j之间的最小/大值。

Part 1:概述

搜索区间最值的方法较为常用的是线段树和RMQ，线段树更加灵活多变，但是写起来比较复杂，RMQ相对来说更加方便书写。RMQ的st算法和其他st算法一样都需要一个预处理的过程，这个过程的时间复杂度$O（Nlog_2{N}）$ ,查询的复杂度$O(1)$。综合来看还是一种较为快捷的方法。

Part 2:分析

RMQ的st算法是基于DP预处理的，假设DP数组为$dp[maxn][20]$ ,原数组为$num[maxn]$;
DP状态中$dp(i,j)$ 表示的是第i个数字起连续$2^j$ 个数中的最值。

例如num数列：1 7 6 3 1 2 5 6 3 8 6
dp(1,0) = max(1)；
dp(1,1) = max(1,7);
dp(1,2) = max(1,7,6,3);

所以DP状态中的初值dp(i,0)也就是第i个数值本身。

状态和初始值确定之后然后就是找到转移方程，我们把dp(i,j)分成两个部分（由上述分析可以得到dp(i,j)中一共有$2^j$ 个数字，一定为偶数，所以每一个部分就是$2^{j-1}$ 个数字）$[i ,i + 2^{j-1}-1]$ 为一个区间,$[i+2^{j-1},i+2^j-1]$ 为第二个区间

上例num数列：1 7 6 3 1 2 5 6 3 8 6
i=1,j=3时两个区间分别是[1,7,6,3]和[1,2,5,6]

那么dp(i,j)就是这个两个区间中的最大值的较大的那个值(最小值同理)

Part 3:处理过程

有上述分析我们可以得到下面的c++代码：

void RMQ(int num) //预处理->O(nlogn)  
{  
    for(int j = 1; j < 20; ++j)  
        for(int i = 1; i <= num; ++i)  
            if(i + (1 << j) - 1 <= num){  
    max_dp[i][j]=max(max_dp[i][j-1],max_dp[i + (1<<(j-1))][j-1]);  
    min_dp[i][j]=min(min_dp[i][j-1],min_dp[i + (1<<(j-1))][j-1]);  
    }  
}  

对于查询来说也就是查询区间$[i,j]$，按照上述的其区间长度为$j-i+1$
去$k = log_2{(j-i+1)}$ 则$RMQ(A, i, j)=max{dp(i , k), dp( j - 2 ^ k + 1, k)}$
所以当我们查询的时候代码如下：

int Mininum(int L,int H){
    int k=(int)(log((double)H-L+1)/log(2.0));
    return min(DPMIN[L][k],DPMIN[H-(1<<k)+1][k]);
}
int Maxinum(int L,int H){
    int k=(int)(log((double)H-L+1)/log(2.0));
    return max(DPMAX[L][k],DPMAX[H-(1<<k)+1][k]);
}

Part 4:例子

51nod1174区间中最大数

#include<bits/stdc++.h>
#define MAXN 100010
#define MAXE 100010
#define ll long long
#define eps 1e-6
#define PI acos(-1)
#define FILEIN freopen("in.txt","r",stdin)
#define FILEOUT freopen("out.txt","w",stdout)
#define CLR(x) memset(x,0,sizeof(x))
#define MEM(a,x) memset(a,x,sizeof(a))
#define PB(x) push_back(x)
#define PF(x) push_front(x)
void sf(int &x){scanf("%d",&x);}
void sf(char &x){scanf("%c",&x);}
void sf(double &x){scanf("%lf",&x);}
void sf(char *x){scanf("%s",x);}
void sf(ll &x){scanf("%I64d",&x);}
void pf(int x){printf("%d",x);}
void pf(char x){printf("%c",x);}
void pf(double x){printf("%f",x);}
void pf(char *x){printf("%s",x);}
void pf(ll x){printf("%I64d",x);}
void pf(){printf("\n");}
const ll mod = 1e9+7;
const int maxn = 100010;
const int INF = 0x7ffffff;
const long long LINF = 0x7fffffffffffffff;
using namespace std;
ll q_pow(ll x,ll y){ll res=1;while(y){if(y&1) res=(res*x)%mod; y>>=1;x=(x*x)%mod;}return res;}
ll inv(ll x){return q_pow(x,mod-2);}
int Read(){int x=0,F=1;char C=getchar();while(C<'0'||C>'9'){if(C=='-')F=-F;C=getchar();}while(C>='0'&&C<='9'){x=x*10-'0'+C,C=getchar();}return x*F;}
int Trans(){char c;sf(c);int res;res=c-'0';return res;}

class RMQ{
public:
    void rmq(int n,int *num)  {
        int temp=(int)(log((double)n)/log(2.0));
        for(int i=0;i<n;i++)  DPMIN[i][0]=DPMAX[i][0]=num[i];
        for(int j=1;j<=temp;j++)
            for(int i=0;i<n;i++)
                if(i+(1<<(j-1))<n) DPMIN[i][j]=min(DPMIN[i][j-1],DPMIN[i+(1<<(j-1))][j-1]);
        for(int j=1;j<=temp;j++)
            for(int i=0;i<n;i++)
                if(i+(1<<(j-1))<n) DPMAX[i][j]=max(DPMAX[i][j-1],DPMAX[i+(1<<(j-1))][j-1]);
    }
    int Mininum(int L,int H){
        int k=(int)(log((double)H-L+1)/log(2.0));
        return min(DPMIN[L][k],DPMIN[H-(1<<k)+1][k]);
    }
    int Maxinum(int L,int H){
        int k=(int)(log((double)H-L+1)/log(2.0));
        return max(DPMAX[L][k],DPMAX[H-(1<<k)+1][k]);
    }
    void Init(){CLR(DPMIN);CLR(DPMAX);}
private:
    int DPMAX[10010][20];
    int DPMIN[10010][20];
};

int num[10010],n,m;
int main(){
    RMQ R;
    R.Init();
    cin>>n;
    for(int i=0;i<n;i++){sf(num[i]);}
    R.rmq(n,num);
    cin>>m;
    while(m--){int l,r;sf(l),sf(r);pf(R.Maxinum(l,r));pf();}
    return 0;
}

Part 5:注意点

….