RMQ ---- ST（Sparse Table）算法

【概述】

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

一般解决这类区间问题可以用线段树来做，比较通用。但是线段树的代码量有点多。这里介绍另一种算法，ST（Sparse Table）算法。

ST（Sparse Table）算法是一个非常有名的在线处理RMQ问题的算法，它可以在O(nlogn)时间内进行预处理，然后在O(1)时间内回答每个查询。

【构建】  
      该算法是基于倍增思想设计的在线算法，之所以称之为倍增，是因为算法记录从每个元素开始的连续的长度为2^k的区间中元素的最值。
      类似于动态规划，F[i, j]表示从第i个数起连续2^j个数中 ( 区间为[i,i+2^j-1] )的最值。状态转移方程F[i, j]=max/min（F[i，j-1], F[i + 2^j-1，j-1]），转移时需要注意别让超界了。
      这个状态转移方程不难想到,一张图就能说明：
【查询】

对于一个查询区间 [i,j] ,只要找到一个或者多个2的整数倍长度的区间覆盖[i,j] ，取这些区间最值的最值就是答案了。

如何把[i,j]覆盖完整？一种办法是把区间的长度按照二进制分成多个2的整数倍区间，显然这些区间是不重叠的，这样求多次最值就能得到答案。不过这种发放增加了算法常数，一次查询可能就要求几十次最值。

还有种更好的方法，原理是：为了减少分割出的区间数量，允许区间重叠，这样所有的情况下最多只要两个区间就好了，见下图：

   只要求出k就好了，k=(int)((log(j-i+1.0)/log(2.0)))。

【代码】

以下代码基于数组下表从0开始。

初始化：

 //n为元数的个数
 //bitn为n的二进制位数，取下整(int)(log(n)/log(2))
 for (int i=; i<n; ++i)
    f[i][]=input[i];
 for (int j=; j<bitn; ++j)
    for (int i=; i<n; ++i)
    {
        if (i+(<<(j-))>=n) break;
        f[i][j]=max(f[i][j-],f[i+(<<(j-))][j-]);
     }

查询：

 int query(int s,int e)  //查询区间[s,e]的最值
 {
     int k=(int)((log(e-s+1.0)/log(2.0)));
     return max(f[s][k],f[e-(<<k)+][k]);
 }