ST算法分析及代码

ST旨在利用倍增思想处理区间最值问题。相比于线段树和树状数组的优势在于只需经过O(nlogn)的预处理即可实现O(1)在线的查询l-r的最值。

分为两个部分，一：预处理，用f[i][j]表示从[i,i+2^j-1]的最值，初始化f[i,0]=a[i]。递推过程如下：将区间[i,i+2^j-1]分成两个已知区间[i,i+2^(j-1)-1]和[i+2^(j-1),i+2^j-1]表示为f[i][j-1],f[i+(1<<(j-1))][j-1]。

二：查询，给出l和r后确定一个k，使2的k次幂小于区间长度的前提下最大，那么从l开始的2^k个数和以k结尾的k个数，这两段一定覆盖了整个区间，这两段的最大值为f[l][k],f[r-(1<<k)+1][k]。

代码如下

#include <bits/stdc++.h>
using namespace std;
const int N=100001;
int n,m,a[N],f[N][20],t,k,l,r;
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
	scanf("%d",&a[i]);
	for(int i=1;i<=n;i++)
	f[i][0]=a[i];
	t=log(n)/log(2);
	for(int j=1;j<=t;j++)
	for(int i=1;i<=n-(1<<j)+1;i++)
	f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
	
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&l,&r);
		k=log(r-l+1)/log(2);
		printf("%d\n",max(f[l][k],f[r-(1<<k)+1][k]));
	}
}