倍增与ST表

总说

一、算法介绍

1.1 ST表——O（nlogn + m）

通常用于快速求出一个子区间的最大值，最小值等等，这里以最大值为例

缺点：ST表是静态的，不支持修改操作。

ST表本质上是基于倍增思想的动态规划。

算法介绍

我们规定：

f[i][j]表示：从i开始向后，长度为2^j的区间中的最大值

如果我们想求f[i][j]，我们可以把区间均分为2部分

2部分区间长度为2^(j - 1)，第二个区间的起点就是：i + 2 ^(j - 1)

则转移方程就为：

f[i][j] = max(f[i][j - 1], f[i + (1 << (j - 1))][j - 1]);

我们预处理完之后，我们改如何查询？如图

假设我们要求[L, R]区间的最大值，区间长度为len，则我们能找到一个最大的k，使得：

即：

则有

如图：

中间虽然有重复，但是由于我们求的最大值，就算重复比较但是并不会改变结果

具体到转移方程如下：

至此，我们就学习完了。

中间对于k的求解，我们可以用cmath中的log(x)函数，

代码表示如下：

int k = log(len) / log(2);

起点有n个，对于每个起点要预处理logn次，有m次询问，所以总时间复杂度为O（nlogn + m）

代码模板

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>

using namespace std;

const int N = 2e5 + 10;
const int M = 18;

int n, m;
int f[N][M], w[N];//f[i][j]表示：从i开始 长度为2^j的区间 中的最大值

void init() // 初始化st表
{   
    //类似区间dp，第一维是区间长度，第二维是起点
    for(int j = 0; j < M; j++) //2^0 到2^M-1 枚举长度
        for(int i = 1; i + (1 << j) - 1 <= n; i++) //枚举起点,右端点不能超出界限
        {
            if(!j) //如果长度为1
                f[i][j] = w[i];
            else //将区间一分为2，长度均为2^(j-1)，取2个区间的最大值
                f[i][j] = max(f[i][j - 1], f[i + (1 << (j - 1))][j - 1]);    
        }
}

int query(int l, int r)
{
    int len = r - l + 1; //区间长度
    int k = log(len) / log(2); //由2^k <= len 得出符合条件的最大的k
    return max(f[l][k], f[r - (1 << k) + 1][k]);
}

int main()
{
    scanf("%d",&n);
    for(int i = 1; i <= n; i++)
        scanf("%d", &w[i]);
    
    init(); //初始化
    scanf("%d",&m);
    while(m--)
    {
        int l, r;
        scanf("%d %d",&l, &r);
        printf("%d\n", query(l, r));
    }
    return 0;
}

二、题目练习

2.1 ST表练习

2.1.1 区间最大值——板子题

题目

题目链接：1270. 数列区间最大值 - AcWing题库

思路

板子题，不再说明

代码模板

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>

using namespace std;

const int N = 1e5 + 10;
const int M = 18;

int n, m;
int f[N][M], w[N];

void init()
{
    for(int j = 0; j < M; j++)
    {
        for(int i = 1; i + (1 << j) - 1 <= n; i++)
        {
            if(!j)
                f[i][j] = w[i];
            else
                f[i][j] = max(f[i][j - 1], f[i + (1 <<(j - 1))][j - 1]);
        }
    }
}

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

int main()
{
    scanf("%d %d", &n, &m);
    for(int i = 1; i <= n; i++)
        scanf("%d", &w[i]);
    init();
    while(m--)
    {
        int l, r;
        scanf("%d %d", &l, &r);
        printf("%d\n", query(l, r));
    }
	return 0;
}