RMQ_ST算法

 

RMQ算法是求给定区间的最值问题，例如给你一段区间A[6]={3,5,2,1,6,7} 求区间[2,4]内的最大数，看的出来最大数为6

但是有时区间很大，大概有几万，并且询问的次数达到几十万次，普通的搜索算法是肯定会超时的。ST算法，全称为Sparse_Table算法，每次查询的时间复杂度可以达到O(1)，这样十几万次的查询就不再话下了！

 

实现的原理（easy!）

它的思想有一部分动态规划的思想，千万不要着急，很简单的动态规划，和递推差不多。假如我们初始的一段序列存储在num[]数组中（以后都是这个一维数组存储要查询的数据）。

我们规定二维的一个数组f[i][j]的值为区间[i,i+2^j-1]内那个最大数，例如num[6]={3,5,2,1,6,7}那么

f[0][0]表示区间[0,0+2^0-1]即区间[0,0];

f[2][2]表示区间[2,2+2^2-1]即区间[2,5];

f[1][0]表示区间[1,1+2^0-1]即区间[1,1];

f[2][0]表示区间[2,2+2^0-1]即区间[2,2];

我们可以看出f[i][0];其实就对应了num[i]的值

注意：f数组的纵坐标不用太大因为因为我们必须保证 i+2^j<=n(n是数组的元素个数)

这样先给出状态转移的方程f[i][j]=min(f[i][j-1], f[i+2^(j-1)][j-1]);

证明如下：

f[i][j]表示区间 [i,i+2^j-1]中最小的元素；

f[i][j-1]表示区间[i,i+2^(j-1)-1]中最小的元素；

f[i+2^(j-1)][j-1]表示区间[i+2^(j-1),i+2^(j-1)+2^(j-1)-1]中最小元素；

看的出两个区间刚好是把区间[i,i+2^j-1]划分成数量相等的两部分，因此分别从这两部分中求出最小值即可。

 

初始化完成f数组的之后，就要到转到利用f数组怎样查询的问题，我们肯定是利用f数组中已有的值！假如我们现在要查询的区间是[m,n],求出这个区间的最小值，对于这个区间[m,n]假如我们能把它分成几个区间，并且这几个区间就是

f[][]保存的值，然后从这几个最小值中选出最小值，然而实际上我们只要分成两个区间即可，我们知道求区间[m,m+2^k-1]中的最小值就是f[m][k]但是这个m+2^k-1必须小于n因此我们尽可能使得m+2^k-1接近于n 即求出满足

2^k<=n-m+1的最大值，那么剩下区间[m+2^k-1,n]难道就不管了么？当然不是，你看n-(2^k-1)>=m,因此我们可以求出区间[n-(2^k-1),n]的最小值，刚好是f[n-(2^k-1)][k]里面的值，这两个区间肯定是有重复区间这是明显的，这样查询就很简单了 。如下图    至于k怎求出，这是一个数学问题 k=(int)(log(doube(n-m+1))/log(2.00))

 

练一下啊？连接http://poj.org/problem?id=3264

 

注意：本题的题目中的数组从1，我们从0..

#include <stdio.h>
#include <iostream>
#include <string.h>
#include <math.h>
using namespace std;

#define mmax(a,b) ((a)>(b)?(a):(b))
#define mmin(a,b) ((a)>(b)?(b):(a))
#define maxn 50050
int num[maxn];
int f1[maxn][32];
int f2[maxn][32];

int N,Q;

void Init()
{
    int k=(int)(log(N*1.0)/log(2.0));
    int i,j,m;
    for(i=0;i<N;i++)
    {
        f2[i][0]=f1[i][0]=num[i];
    }
    for(j=1;j<=k;j++)
    {
        for(i=0;i+(1<<j)-1<N;i++)
        {
            m=i+(1<<(j-1));
            f1[i][j]=mmax(f1[i][j-1],f1[m][j-1]);
            f2[i][j]=mmin(f2[i][j-1],f2[m][j-1]);
        }
    }
}
void RMQ(int i,int j)
{
    int k=(int)(log(double(j-i+1))/log(2.0)),t1,t2;
    t1=mmax(f1[i][k],f1[j-(1<<k)+1][k]);
    t2=mmin(f2[i][k],f2[j-(1<<k)+1][k]);
    printf("%d\n",t1-t2);
}

int main()
{
    int A,B;
   while( scanf("%d%d",&N,&Q)!=EOF)
   {
       memset(f1,0,sizeof(f1));
       memset(f2,0,sizeof(f2));
       for(int i=0;i<N;i++)
        scanf("%d",&num[i]);
       Init();
       for(int i=0;i<Q;i++)
       {

           scanf("%d%d",&A,&B);
           RMQ(A-1,B-1);
       }
   }
    return 0;
}