学习记录 ST表

稀疏表（ST表）也是一种用于RMQ的数据结构，它利用了倍增与动态规划的思想。如果只存在区间查询，不存在区间查找，则ST表是一种不错的选择。

为了能在O(1)的时间完成查询，可以这样预处理：如果记maxn[i][j]为[i,j]区间上的最大值，易得max[i][j]=maxn([i][j-1],a[i])。但此种方法预处理的复杂度为$O(n^2)$，为了简化，我们可以用倍增思想。记st[i][j]为从i开始连续$2^j$个元素的最大值。则有
$st[i][j]=max(st[i][j-1],st[i+2^{j-1}][j-1])$
最后求解[l,r]上的最大值时，可以将结果看成从l开始向后数f(len)个的最大值到从r开始向前数f(len)个的最大值的最大值。其中f(len)=$\lfloor lo{g}_2(r-l+1)\rfloor$。很显然两个区间有交集，不过这并不影响结果。
ST表详细讲解
代码如下：

#include <cstdio>
#include <algorithm>
const int M=(int)1e5+1;
using namespace std;
int a[M],st[M][17];
int lg2[M];

int main() {
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1; i<=n; i++) {
    	scanf("%d",&a[i]);
    	st[i][0]=a[i];
	}
    for(int i=2; i<=n; i++)    lg2[i]=lg2[i/2]+1;
    //为防止用后面的结果更新前面的，此处应先枚举长度
    for(int j=1; j<=lg2[n]; j++) {
    	for(int i=1; i+(1<<j)-1<=n; i++) {
    		st[i][j]=max(st[i][j-1],st[i+(1<<(j-1))][j-1]);
		}
	}
	while(m--) {
		int l,r;
		scanf("%d%d",&l,&r);
		int lglen=lg2[r-l+1];
		printf("%d\n",max(st[l][lglen],st[r-(1<<lglen)+1][lglen]));
	}
    return 0;
}

其实也用到了区间DP的思想。