RMQ问题:
- 给定一个长度为n的区间(1-n),m次询问,每次询问[l,r]区间的 极值。是一种离线方法。
ST表:
- 一种利用dp求解区间最值的倍增算法。
关键:
- 构造f[i][j] 数组代表的是以第i个元素为起点的长度为2^j的区间极值。
预处理:
-
f[i][0] = s[i],长度为1的区间的极值是元素本身。
-
状态转移:每个区间的长度都是2的次方,因此区间可以均分为两个等长区间,该区间的极值就等于均分的两个区间极值的极值。
-
第一区间:从i开始,长度为2^(j-1)–>f[i][j-1]
-
第二区间:从i+2^ (j-1)开始,长度为2^(j-1)–>f[i+(1<<(j-1))][j-1])
-
转移方程:f[i][j] = max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
查询:
- 区间[l,r],长度为k = r-l+1,查询两个区间长度就是log2k,一个区间是[l,l+2^k-1],另一个区间是[l+2 ^k,r]。
代码:
#include<bits/stdc++.h>
using namespace std;
int f[1010][110];
int s[1010];
int n;
void ST() // 预处理区间极值
{
//f[i][j]表示以第i个元素为起点的2^j长度的区间极值
for(int i = 1; i <= n; i++)
{
f[i][0] = s[i]; // 长度为1的区间极值为它本身
}
for(int j = 1; (1<<j) <= n; j++) // 枚举区间2^j的长度(长度小于总区间长度)
{
for(int i = 1; i+(1<<j)-1 <= n; i++) // 枚举起点(起点加区间长度小于总长度)
{
// f[i][j] : s[i]-s[i+2^j-1]
// f[i][j-1] : s[i]-s[i+2^(j-1)-1]
// f[i+(1<<(j-1))][j-1] : s[i+(1<<(j-1))]-s[i+(1<<(j-1))+2^(j-1)]
/*
由上面三个区间可以得出由于区间长度都是2的次方
所以当前区间可以分割为两个长度为2^(j-1)的区间