ST算法(倍增)

在RMQ问题(区间最值问题)中,著名的ST算法就是倍增的产物。给定一个长度为N的数列A,ST算法能在O(N\log N)时间的预处理后,以O(1)的时间复杂度在线回答“数列A中下标在 l~r之间的数的最大值是多少”这样的区间最值问题。

一个序列的子区间个数显然有O(N^{2})个,根据倍增思想,我们首先在这个规模为O(N^{2})的状态空间里选择一些2的整数次幂的位置作为代表值。

F[i,j]表示数列A中下标在子区间[i,i+2^{j}-1]里的数的最大值,也就是从i开始的2^{j}个数的最大值。递推边界显然是F[i,0]=A[i],即数列A在子区间[i,i]里的最大值。

在递推时,我们把子区间的长度成倍增长,有公式f[i,j]=max(F[i,j-1],F[i+2^{j-1}][j-1]),即长度为2^{j}的子区间的最大值是左右两半长度为2^{j-1}的子区间的最大值中较大的一个。

void ST_perwork(int n){
    for(int i=1;i<=n;i++) f[i][0]=A[i];
    int 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]);
        }
    }
}

当询问任意区间[l,r]的最值时,我们先计算出一个k,满足2^{k},也就是使2的k次幂小于区间长度的前提下最大的k。那么“从l开始的2^{k}个数”和“以r结尾的2^{k}个数”这两段一定覆盖了整个区间[l,r],这两段的最大值分别是F[l,k]F[r-2^{k}+1,k],二者中较大的那个就是整个区间[l,r]的最值。因为求的是最大值,所以这两段只要覆盖区间[l,r]即可,即使有重叠也没关系。

int ST_query(int l,int r){
    int k=log(r-l+1)/log(2);
    return max(f[l][k],f[r-(1<<k)+1][k]);
}

简便起见,我们再代码中使用了cmath库的log函数。该函数效率较高,一般来说对程序的性能影响不大。更严格的讲,为了保证复杂度为O(1),应该O(N)预处理1~N这N种区间长度各自对应的k值,在询问时直接使用。

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值