ST表用以解决RMQ问题(给定一个长度为N的区间,M个询问,每次询问Li到Ri这段区间元素的最大值/最小值)。复杂度nlog(n)
原理: 预处理:
· 一种利用dp求解区间最值的倍增算法。区间在增加时,每次并不是增加一个长度,而是使用倍增的思想,每次增加2^i个长度。
·设二维数组st[i][j],表示需要查询的数组的从下标i到向后延伸2^j - 1长度(即下标i+2^j - 1)的最值。
·在求解st[i,j]时,ST算法是先对长度为2^j的区间[i,i + 2^j - 1]分成两等份,每份长度均为2^(j - 1)。之后在分别求解这两个区间的 最值st[i,j - 1]和st[i + 2^(j - 1),j - 1]。,最后在结合这两个区间的最值,求出整个区间的最值。 则得到状态方程:
int a[1010];//原始输入数组
int st[1010][20];//st表
void init(int n)
{
for (int i = 0; i < n; i++)//预处理st[][]
st[i][0] = a[i];
for (int j = 1; (1 << j) <= n; j++)
{
for (int i = 0; i + (1 << j) - 1 < n; i++)
st[i][j] = min(st[i][j - 1],st[i + (1 << (j - 1))][j - 1]);
// st[i][j] = min(st[i][j - 1],st[i + 2^(j - 1))][j - 1])
}
}
区间查询:需要查询的区间为[l,r],则需要找到两个相等的覆盖这个闭区间的最小幂区间。
这两个区间可以重复,因为两个区间是否相交对区间最值没有影响。根据上述思想,可以把待查询区间[l,r]分成两个小区间[l,l + 2^k- 1] 和 [r - 2^k + 1,r] ,其又分别对应着st[l,k]和st[r - 2^k + 1,k],此时为了求解整个区间的最小值,我们只需求这两个值得最小值即可。
int search(int l, int r)
{
int k = (int)(log((double)(r - l + 1)) / log(2.0));
return min(st[l][k],st[r - (1 << k) + 1][k]);
}