倍增 ST表

引入

直接通过一道经典的例题来引入st表以及倍增的概念。

落谷P3865 ST表

描述

给定一个长度为N的数列,和M次询问,求出每一次询问的区间内数字的最大值。

输入格式

第一行包含两个整数N,M(1 \leq N \leq 10^5,1 \leq M \leq 2 \times 10^6),分别表示数列的长度和询问的个数。

第二行包含N个整数(记为 a_i(0 \leq a_i \leq 10^9)),依次表示数列的第i项。

接下来M行,每行包含两个整数l_i, r_i,表示查询的区间为[l_i, r_i]

输出格式

输出包含M行,每行一个整数,依次表示每一次询问的结果。

输入样例

8 8
9 3 1 7 5 6 0 8
1 6
1 5
2 7
2 6
1 8
4 8
3 7
1 8

输出样例

9
9
7
7
9
8
7
9

思路

如果用传统的思路去写这道题,那么无非就是每次都对区间进行一次遍历,那么复杂度将高达N\times M,那一定会T掉,所以我们期望是先进行一遍预处理,然后每次通过O(logn)或者是O(1)的复杂度进行查询,从而挤进去。

倍增

这边,则需要引入一个倍增的思想,利用倍增进行预处理,则可大大节省查询所需的时间,也就是查询可以达到O(1)的复杂度。

为了达到O(1)的查询时间复杂度,那么我们就需要预处理一个表,蕴含所有区间的结果。

而倍增的独特之处在于,他在预处理的时候时间复杂度仅仅只需要O(nlogn),而查询也只需要O(1)

定义

为了预处理出来一个表,那么我们先定义一下这个表的含义。

我们定义一个二维数组:st[i][j]是从i开始,一直到i+2^j这个区间的最大(小)值,也就是st[i][j]=max \{ [i,i+2^j) \}

那么这个有什么用呢,先不着急,我们先想办法预处理出来这个表。

预处理

设我们的原数组是a[],那么首先很显然的是st[i][0] = a[i];

之后我们在进行递推,这边有一个公式st[i][j] = max\{st[i][j - 1],st[i + 2^{j-1}][j - 1]\}

乍一看可能不是那么好理解,其实根据定义再结合一下几个直观的数字便可很容易明白他在做什么。


我们首先需要明白的是,一个区间内,不管是最小值还是最大值,该区间的最大(小)值都等于其所有子区间的最大(小)值的最大(小)值,即满足max\{range\}=max\{max\{subrange_1\},max\{subrange_2\},...,max\{subrange_n\}\}

那么上述公式也很好理解了,根据定义st[i][j]是从i开始到i+2^j这个区间内的最大(小)值,所以他可以裂开成两个子区间,为了能够递推,所以我们裂开平分出两个长度为2^{j-1}的子区间,很显然,这两个区间分别为[i,i+2^{j-1})[i+2^{j-1},i + 2^{j-1} + 2^{j-1}),也即[i+2^{j-1},i+2^j),因此这两个区间恰好把原区间[i,i + 2^j)给分开了。

图示

剩下的便是在代码上的处理了,代码不难理解,直接上吧。

#define rep(i,a,b) for (int i = a;i <= b;i ++)

const int MAXN = 100010;

int st[MAXN][20];
int a[MAXN];

void init() {
    // 定义 st[i][j] 是从i开始,到i + 2^j这一段,即[i,i + 2^j]这一段中的最大/小值
    rep(i,1,n) st[i][0] = a[i];

    for (int j = 1;(1 << j) <= n;j ++) { // 遍历所有的j,j是一个很小的数字,最大值=log2(n)
        rep(i,1,n - (1 << j) + 1) { // 在[1,n]区间范围内,确定j的情况下,把所有的i都遍历求值一遍
            st[i][j] = max(st[i][j - 1], st[i + (1 << (j - 1))][j - 1]); // 套公式
        }
    }
}

查询

查询这边也相应着有一个公式,若需要查询的区间左端下标是l,右端是r,则有x=\lfloor log_2(r-l+1)\rfloor ans=max\{st[l][x],st[r - 2^x + 1][x]\}

解释

我们想查询的区间是[l,r],也就是[l,r+1),st表的定义是st[i][j]=max\{[i,i + 2^j)\}

那么上述公式中st[l][x] = max\{[l,l + 2^x)\} st[r - 2^x + 1][x]=max\{[r - 2^x + 1,r+1)\}

x=log_2(r-l+1),那么上述公式描述的均是max\{[l,r+1)\}

而实际上x=\lfloor log_2(r-l+1)\rfloor,而2^{\lfloor log_2a\rfloor}>\frac a2
所以答案ans只需要求[l,l+2^x)[r - 2^x + 1,r+1)这两段的区间最大值即可。

查询代码

int query(int l, int r)
{
    int x = log2(r - l + 1);
    return max(st[l][x],st[r - (1 << x) + 1][x]);
}

st模版

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <vector>
#include <list>
#include <set>
#include <utility> // pair
#include <map>
#include <iostream>
#include <sstream>
#include <algorithm> // sort
#include <string>
#include <stack>
#include <queue>
#include <fstream>
#include <bitset>

using namespace std;

#define ll long long
#define lll __int128
#define uchar unsigned char
#define ushort unsigned short
#define uint unsigned int
#define ulong unsigned long
#define ull unsigned long long

#define INT_INF 0x7fffffff

#define pi acos(-1)

#define mx(a,b) (a) > (b) ? (a) : (b)
#define mn(a,b) (a) < (b) ? (a) : (b)
#define mem(a,b) memset(a,b,sizeof(a))
#define fre(a) freopen(a,"r",stdin)

#define cio ios::sync_with_stdio(false); // Do not use it with "scanf" and other c input!
#define pb push_back
#define rep(i,a,b) for (int i = a;i <= b;i ++)
#define pre(i,a,b) for (int i = a;i >= b;i --)
#define REP(i,a,b) for (int i = a;i < b;i ++)

#define read(a,s,n) rep(i,s,n) scanf("%d",a + i);
#define READ(a,s,n) REP(i,s,n) scanf("%d",a + i);

#define read_ll(a,s,n) rep(i,s,n) scanf("%lld",a + i);
#define READ_ll(a,s,n) REP(i,s,n) scanf("%lld",a + i);

#define _T_(T) int T;scanf("%d",&T);while (T --)
#define _E_(T) while (~T)

#define Tprint(a,s,e) {int f=1;REP(i,s,e){if(f)f=0;else printf(" ");printf("%d",a[i]);}}


#define endl '\n'

#define itn int
#define nit int
#define inr int
#define mian main
#define ednl endl
#define fro for
#define fir for
#define reutrn return
#define retunr return
#define reutnr return


/* header_useful_h */

inline int read_int()
{
    int x = 0,f = 1;
    char ch = getchar();
    while (!isdigit(ch)) {
        if (ch == '-') f = -1;
        ch = getchar();
    }
    while (isdigit(ch))
    {
        x = x * 10 + ch - 48;
        ch = getchar();
    }
    return x * f;
}

const int MAXN = 100010;

int st[MAXN][20];
int a[MAXN];

int n,m;

void init() {
    // 定义 st[i][j] 是从i开始,到i + 2^j这一段,即[i,i + 2^j]这一段中的最大/小值
    rep(i,1,n) st[i][0] = a[i];

    for (int j = 1;(1 << j) <= n;j ++) { // 遍历所有的j,j是一个很小的数字,最大值=log2(n)
        rep(i,1,n - (1 << j) + 1) { // 在[1,n]区间范围内,确定j的情况下,把所有的i都遍历求值一遍
            st[i][j] = max(st[i][j - 1], st[i + (1 << (j - 1))][j - 1]); // 套公式
        }
    }
}

int query(int l, int r)
{
    int x = log2(r - l + 1);
    return max(st[l][x],st[r - (1 << x) + 1][x]);
}

int main()
{
    scanf("%d%d",&n,&m);
    rep(i,1,n) {
        a[i] = read_int();
    }
    init();
    int l,r;
    while (m --) {
        scanf("%d%d",&l,&r);
        printf("%d\n",query(l, r));
    }
    return 0;
}
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 202,056评论 5 474
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 84,842评论 2 378
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 148,938评论 0 335
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,296评论 1 272
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,292评论 5 363
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,413评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,824评论 3 393
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,493评论 0 256
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,686评论 1 295
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,502评论 2 318
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,553评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,281评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,820评论 3 305
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,873评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,109评论 1 258
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,699评论 2 348
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,257评论 2 341

推荐阅读更多精彩内容