ST表,稀疏表,用于求解经典的RMQ问题。即区间最值问题。
Problem:
给定n个数和q个询问,对于给定的每个询问有l,r,求区间[l,r]的最大值。.
Solution:
主要思想是倍增和区间dp。
状态:dp[i][j] 为闭区间[i,i+2^j-1]的最值。
这个状态与转移方程的关系很大,即闭区间的范围涉及到了转移方程的简便性。
转移方程:dp[i][j]=max(dp[i][j-1],dp[i+2^(j-1)][j-1])。
这是显然的,但这里有个细节:第一个项的范围为[i,i+2^(j-1)-1],即第i个数到第i+2^(j-1)个数的前一个数,而第二个项的范围是[i+2^(j-1),i+2^j-1]。这里容易弄混,所以导致无法理解整个方程,或者写错。
询问:[l,r]区间的最大值。
令g=log2(r-l+1)向下取整。则区间的最值就是max(dp[l][g],dp[r-2^g+1][g])。
好的,我们开始对这些方程做一番分析。假设现在有5个数,为如下情况。这些分别是,dp[1][0],dp[2][0]....dp[5][0]
接着,我们可以统计到红色区域的最大值。这些分别是dp[1][1],dp[2][1]...dp[4][1]
最后,我们再统计到dp[1][2],dp[2][2],也就是1~1+4-1和2~2+4-1的最大值。
直到整个求完。(5后面的点因为皆为0,所以在取区间最大值的时候可以直接忽略)
有了这些后,如果要求l~r的最大值,我们需要让两个区间覆盖这个区间。假设我们询问1~5的最大值。
按道理来讲,我们假设2^x=5-1+1,然后dp[1][x]不就完事了?然而很可惜,x可能为小数。因此,x我们应当向下取整x',那么显然这样子取dp[1][x']会比dp[1][x]少一段。
于是我们再采用dp[5-2^x'+1][x']的方式来弥补。什么意思呢?
像这样,在两个互相被覆盖的区间里取一个最大值,是不是完美的解决了这个问题呢?而这个x就等于log2(x),x'即为log2(x)向下取整。
你可能会说这样有没有可能会没有全部覆盖,这很显然不可能。因为l+2^x'<r-2^x' 等价于 2^(x'+1)<r-l。而显然这个x'+1>=x,所以2^(x‘+1)>2^x=r-l+1>r-l。因此原不等式绝对不成立。
以上,ST表是一个很好的倍增思想入门。在LCA中也用到了与ST表非常类似的倍增思想。
ST表很简单。请注意常数,很容易就会了。Luogu P3865即为模板题。
下面是代码:
#include<bits/stdc++.h>
using namespace std;
template <typename T> T gn(T &x){
x=;
T plus=;
char ch=getchar();
while(ch<''||ch>'')plus=(ch=='-'?-:),ch=getchar();
while(ch>=''&&ch<='')x=x*+ch-'',ch=getchar();
return x*plus;
}
const int N=1e5+;
int a[N][];
template <typename T> inline int mmax(T &a,T &b){return a<b?b:a;}
inline int dpow(int a,int x){
int ret=;
while(x){
if(x&)ret*=a;
a*=a;
x>>=;
}
return ret;
} int main(){
int n,m;
gn(n),gn(m);
for(int i=;i<=n;i++)
gn(a[i][]);
int plus=;
for(int j=;j<;j++){
for(int i=;i+plus<=n;i++)
a[i][j]=mmax(a[i][j-],a[i+plus][j-]);
plus*=;
}
int na,nb,x;
for(int i=;i<m;i++){
gn(na),gn(nb);
x=log2(nb-na+);
printf("%d\n",mmax(a[na][x],a[nb-dpow(,x)+][x]));
}
return ;
}