RMQ五大神器：ST表，树状数组，线段树，主席树，莫队

最新推荐文章于 2024-03-10 11:08:00 发布
最新推荐文章于 2024-03-10 11:08:00 发布
阅读量465 收藏
点赞数 1
分类专栏： 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/codeisking/article/details/50839526
版权

一、ST表

题目一 http://poj.org/problem?id=3264

#include<cstdio>
#include<cmath>
#include<algorithm>
#define N 100010
using namespace std;
int st1[N][30];
int st2[N][30];
int n,m,a,b;
void ST(int n){
    for(int j=1;j<=30;j++)
        for(int i=1;i+(1<<j)-1<=n;i++){
            st1[i][j]=max(st1[i][j-1],st1[i+(1<<(j-1))][j-1]);
            st2[i][j]=min(st2[i][j-1],st2[i+(1<<(j-1))][j-1]);
        }
}
int MX(int l,int r){
    int k=(int)(log(r-l+1.0)/log(2.0));
    return max(st1[l][k],st1[r-(1<<k)+1][k]);
}
int MN(int l,int r){
    int k=(int)(log(r-l+1.0)/log(2.0));
    return min(st2[l][k],st2[r-(1<<k)+1][k]);
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%d",&st1[i][0]);
        st2[i][0]=st1[i][0];
    }
    ST(n);
    while(m--){
        scanf("%d%d",&a,&b);
        printf("%d\n",MX(a,b)-MN(a,b));
    }
    return 0;
}

二、树状数组

 

 

题目一 http://ybt.ssoier.cn:8088/problem_show.php?pid=1535

【例 1】数列操作
【题目描述】
给定n个数列,规定有两种操作，一是修改某个元素，二是求子数列[a,b]的连续和。数列元素个数最多10万个，询问操作最多10万次。

【输入】
第一行2个整数n,m(n表示输入n个数，m表示m操作）
第二行n个整数
接下来m行，每行三个数k,a,b(k=0,表示求子数列[a,b]的连续和；k=1,表示第a个数加b)。

【输出】
若干行,表示k=0时，对应子数列[a,b]连续和。

【输入样例】
10 5
1 2 3 4 5 6 7 8 9 10
1 1 5
0 1 3
0 4 8
1 7 5
0 4 8
【输出样例】
11
30
35
 

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10;
int n,m,k,a,b,C[N];
int lowbit(int x){return x&-x;}
void add(int x,int y){while(x<=n){C[x]+=y;x+=lowbit(x);}}
int ask(int x){int z=0;while(x){z+=C[x];x-=lowbit(x);}return z;}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){scanf("%d",&a);add(i,a);}
    while(m--){
        scanf("%d%d%d",&k,&a,&b);
        if(k==1)add(a,b);
        else printf("%d\n",ask(b)-ask(a-1));
    }
    return 0;
}

题目二 http://poj.org/problem?id=2352

#include<stdio.h>
#define N 32010
int n,c[N],ans[N];
int lowbit(int x){return x & -x;}
void add(int x){while(x<N){c[x]++;x+=lowbit(x);}}
int ask(int x){
    int s=0;
    while(x){s+=c[x];x-=lowbit(x);}
    return s;
}
int main(){
    scanf("%d",&n);
    for(int i=1,x,y;i<=n;i++){
        scanf("%d%d",&x,&y);
        ans[ask(x+1)]++;
        add(x+1);
    }
    for(int i=0;i<n;++i)printf("%d\n",ans[i]);
    return 0;
}

 

 

题目三 http://ybt.ssoier.cn:8088/problem_show.php?pid=1537

【例 3】校门外的树
【题目描述】
校门外有很多树，学校决定在某个时刻在某一段种上一种树，保证任一时刻不会出现两段相同种类的树，现有两种操作：
K=1，读入 l,r 表示在 l 到 r 之间种上一种树，每次操作种的树的种类都不同；
K=2，读入 l,r 表示询问 l 到 r 之间有多少种树。
注意：每个位置都可以重复种树。

【输入】
第一行 n,m 表示道路总长为 n，共有 m 个操作；
接下来 m 行为 m 个操作。

【输出】
对于每个 k=2 输出一个答案。

【输入样例】
5 4
1 1 3
2 2 5
1 2 4
2 3 5
【输出样例】
1
2
【提示】
数据范围与提示：
对于 20% 的数据，1≤n,m≤100；
对于 %60% 的数据，1≤n≤103,1≤m≤5×104 ；
对于 %100% 的数据，1≤n,m≤5×104 ，保证 l,r>0。

#include<bits/stdc++.h>
using namespace std;
const int N=5e4+10;
int n,m,k,l,r,sum1[N],sum2[N];
int lowbit(int x){return x&(-x);}
void add1(int x,int val){ //更新左端点的数目 
    while(x<=N){
        sum1[x]+=val;
        x+=lowbit(x);
    }
}
void add2(int x,int val){  //更新右端点的数目 
    while(x<=N){
        sum2[x]+=val;
        x+=lowbit(x);
    }
}
int ask1(int x){  //查找区间[1,x]中左端点的数目 
    int res=0;
    while(x){
        res+=sum1[x];
        x-=lowbit(x);
    }
    return res;
}
int ask2(int x){  //查找区间[1,x]中右端点的数目 
    int res=0;
    while(x){
        res+=sum2[x];
        x-=lowbit(x);
    }
    return res;
}
int main(){
    scanf("%d%d",&n,&m);
    while(m--){
        scanf("%d%d%d",&k,&l,&r);
        if(k==1){add1(l,1);add2(r,1);}
        else printf("%d\n",ask1(r)-ask2(l-1));
    }
    return 0;
}

题目四 http://acm.hdu.edu.cn/showproblem.php?pid=4325

#include<bits/stdc++.h>
using namespace std;
const int N=100000+10;
struct Node{int b,e;}A[N];
int _,n,m,a[N*4],B[N],z[N],c;

int lowbit(int x){return x&-x;}
void add(int x,int y){for(int i=x;i<=c;i+=lowbit(i))z[i]+=y;}
int ask(int x){int y=0;while(x){y+=z[x];x-=lowbit(x);}return y;}
int main(){
    scanf("%d",&_);
    for(int cas=1;cas<=_;cas++){
        scanf("%d%d",&n,&m);
        memset(z,0,sizeof(z));
        c=0;
        for(int i=1;i<=n;i++){
            scanf("%d%d",&A[i].b,&A[i].e);
            a[++c]=A[i].b;a[++c]=A[i].e;
        }
        for(int i=1;i<=m;i++){scanf("%d",&B[i]);a[++c]=B[i];}
        sort(a+1,a+c+1);//排序
        c=unique(a+1,a+c+1)-a-1;//去重
        for(int i=1;i<=n;i++){
            int x=lower_bound(a+1,a+c+1,A[i].b)-a;
            int y=lower_bound(a+1,a+c+1,A[i].e)-a;
            add(x,1);add(y+1,-1);
        }
        printf("Case #%d:\n",cas);
        for(int i=1;i<=m;i++){
            int t=lower_bound(a+1,a+c+1,B[i])-a;
            printf("%d\n",ask(t));
        }
    }
    return 0;
}

 

三、线段树

题目一 http://ybt.ssoier.cn:8088/problem_show.php?pid=1537

【 例 1】区间和
【题目描述】
给定一数列，规定有两种操作，一是修改某个元素，二是求区间的连续和。

【输入】
输入数据第一行包含两个正整数n,m(n≤100000,m≤500000),以下是m行，
每行有三个正整数k,a,b(k=0或1,a,b≤n).k=0时表示将a处数字加上b,k=1时表示询问区间[a,b]内所有数的和。

【输出】
对于每个询问输出对应的答案。

【输入样例】
10 20
0 1 10
1 1 4
0 6 6
1 4 10
1 8 9
1 4 9
0 10 2
1 1 8
0 2 10
1 3 9
0 7 8
0 3 10
0 1 1
1 3 8
1 6 9
0 5 5
1 1 8
0 4 2
1 2 8
0 1 1
【输出样例】
10
6
0
6
16
6
24
14
50
41
 


#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=400010;
int n,m,k,a,b;ll c[N];
ll query(int k,int l,int r,int x,int y){
    if(r<x||y<l)return 0;
    if(x<=l&&r<=y)return c[k];
    int mid=(l+r)/2;
    return query(2*k,l,mid,x,y)+query(2*k+1,mid+1,r,x,y);
}
void update(int k,int l,int r,int x,int v){
    if(x<l||x>r)return;
    if(l==r){c[k]+=v;return;}
    int mid=(l+r)/2;
    update(2*k,l,mid,x,v);
    update(2*k+1,mid+1,r,x,v);
    c[k]=c[2*k]+c[2*k+1];
}
int main(){
    scanf("%d%d",&n,&m);
    while(m--){
        scanf("%d%d%d",&k,&a,&b);
        if(k)printf("%lld\n",query(1,1,n,a,b));
        else update(1,1,n,a,b);
    }
    return 0;
}

题目二 http://poj.org/problem?id=3468

 

#include<cstdio>
#include<algorithm>
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define LL long long
using namespace std;
 
const int maxn=111111;
LL sum[maxn<<2],add[maxn<<2];
void PushUp(int rt){ sum[rt]=sum[rt<<1]+sum[rt<<1|1];}

void PushDown(int rt,int m){
    if(add[rt]){
        add[rt<<1]+=add[rt];
        add[rt<<1|1]+=add[rt];
        sum[rt<<1]+=add[rt]*(m-(m>>1));
        sum[rt<<1|1]+=add[rt]*(m>>1);
        add[rt]=0;
    }
}

void build(int l,int r,int rt){
    add[rt]=0;
    if(l==r){scanf("%lld",&sum[rt]); return; }
    int m=(l+r)>>1;
    build(lson);
    build(rson);
    PushUp(rt);
}

void update(int L,int R,int c,int l,int r,int rt){
    if(L<=l && r<=R){
        add[rt]+=c;
        sum[rt]+=(LL)c*(r-l+1);
        return ;
    }
    PushDown(rt,r-l+1);
    int m=(l+r)>>1;
    if(L<=m) update(L,R,c,lson);
    if(m<R) update(L,R,c,rson);
    PushUp(rt);
}

LL query(int L,int R,int l,int r,int rt){
    if(L<=l && r<=R) return sum[rt];
    PushDown(rt,r-l+1);
    int m=(l+r)>>1;
    LL ret=0;
    if(L<=m) ret+=query(L,R,lson);
    if(m<R) ret+=query(L,R,rson);
    return ret;
}

int main(){
    int N,Q;
    scanf("%d%d",&N,&Q);
    build(1,N,1);
    while(Q--){
        char op[2];
        int a,b,c;
        scanf("%s",op);
        if(op[0]=='Q'){
            scanf("%d%d",&a,&b);
            printf("%lld\n",query(a,b,1,N,1));
        }
        else{
            scanf("%d%d%d",&a,&b,&c);
            update(a,b,c,1,N,1);
        }
    }
    return 0;
}

//树状数组

#include<cstdio>
typedef long long ll;
const int N=1e5+10;
ll sum1[N],sum2[N],sum[N];
int num,n,q,a,b,c;
char com[2];
inline int lowbit(int x){return x&(-x);}
inline void add(int x,int v){int i=x;while(i<=n)sum1[i]+=v,sum2[i]+=x*v,i+=lowbit(i);}
inline void range_add(int a,int b,int c){add(a,c);add(b+1,-c);}
inline ll ask(int x){ll ret=0;int i=x;while(i)ret+=(x+1)*sum1[i]-sum2[i],i-=lowbit(i);return ret;}
inline ll range_ask(int a,int b){return ask(b)-ask(a-1);}
int main(){
    scanf("%d%d",&n,&q);
    for(int i=1;i<=n;i++)scanf("%d",&num),sum[i]=sum[i-1]+num;
    while(q--){
        scanf("%s",com);
        if(com[0]=='C')scanf("%d%d%d",&a,&b,&c),range_add(a,b,c);
        else scanf("%d%d",&a,&b),printf("%lld\n",sum[b]-sum[a-1]+range_ask(a,b));
    }
    return 0;
}

 

 

四、主席树

 

五、莫队

codeisking
关注
专栏目录
st表、树状数组、线段树和主席树
weixin_43310882
08-10 496
//树状数组 :单点修改 （+ - * /）； // 区间修改（ + - 最值）其中+ -要利用查分数组来实现 // 单点查询 传统方法 // 区间查询 最值得话就循规蹈矩的进行查找 /* 但是如果使用区间修改过（查分数组）后再进行区间查询，则要进行计算一下 ∑ni = 1A[i] = ∑ni = 1 ∑ij = 1D[j]; 则A[1]+A[2]+...+A[n] ...
RMQ问题（ST表和线段树）
Rrrrya
07-14 359
RMQ问题：给定一个长度为n的区间（1-n），m次询问，每次询问[l,r]区间的极值。 ST表：一种利用dp求解区间最值的倍增算法。 关键：构造f[i][j] 数组代表的是以第i个元素为起点的长度为2^j的区间极值。 预处理： f[i][0] = s[i]，长度为1的区间的极值是元素本身。 状态转移：每个区间的长度都是2的次方，因此区间可以均分为两个等长区间，该区间的极值就等于均分的两个区间极值的极值。 f[i][j] 第一区间：从i开始，长度为2^(j-1)–>f[i][j-1] 第二区间：从i+2
ST表、树状数组、线段树
weixin_54385104的博客
08-23 533
适用范围 ST表： 静态区间查询，并且仅仅适用于max,min,gcd等类型的区间查询。 树状数组 可以进行单点修该，区间查询（或者利用差分进行区间修改，单点查询） 常用于求前缀和。 线段树 都可以，但是我目前只会单点修改，区间查询。 ST表 天才的记忆 从前有个人名叫 WNB，他有着天才般的记忆力，他珍藏了许多许多的宝藏。 在他离世之后留给后人一个难题（专门考验记忆力的啊！），如果谁能轻松回答出这个问题，便可以继承他的宝藏。 题目是这样的：给你一大串数字（编号为 1 到 N，大小可不一定哦！），在你看过一
线段树、树状数组、ST表
学习笔记
07-23 349
#include <iostream> using namespace std; const int n = 6; int arr[] = {1, 3, 5, 7, 9, 11}; int tree[n * 4 + 1]; /** * * @param node tree中的下标 * @param start arr中的下标 表示区间 [start, end] * @param end */ void buildTree(int node, int start...
【算法笔记】【专题】RMQ 问题：ST表/树状数组/线段树
天生我材必有用，千金散尽还复来。
01-05 2112
好久没更算法笔记专栏了，正好学了新算法来更新…… 本文将讲解用于求解 RMQ 问题的三种数据结构： - ST 表（Sparse Table） - 树状数组 - 线段树
st表、树状数组与线段树 笔记与思路整理
weixin_30352645的博客
09-21 368
已更新（2/3）：st表、树状数组 st表、树状数组与线段树是三种比较高级的数据结构，大多数操作时间复杂度为O(log n)，用来处理一些RMQ问题或类似的数列区间处理问题。 一、ST表（Sparse Table） st表预处理时间复杂度O(n log n)，查询O(1)，但不支持在线更改，否则要重新进行预处理。 使用一个二维数组：st[i][j]存储i为...
二维树状数组与数据结构优化：RMQ与线段树解析
静态数组可以用来实现线段树，通过数组 St[] 来表示二叉树结构，其中 lson(x) 和 rson(x) 分别表示节点 x 的左子节点和右子节点。构建线段树的过程包括递归地对子区间进行初始化。 线段树的主要操作包括构建（build...
RMQ问题解析：线段树与ST算法
"该资源是关于RMQ问题的讲解，主要介绍了三种解决方法：朴素搜索、线段树和ST算法。RMQ问题是指在给定数列中查询区间内最小值或最大值。朴素方法时间复杂度高，线段树在处理和查询上具有更优效率，而ST算法（实际上...
LCA与RMQ详解：从树到数组的转换与经典算法解析
在LCA与RMQ之间进行转换时，关键步骤是通过深度优先搜索遍历树，将节点及其深度存储到数组S中，并记录节点首次出现的下标R[i]和深度d(i)。这有助于后续的查询，因为可以利用这些信息来快速确定区间内的最值或公共...
数列区间最大值-RMQ/线段树/树状数组
Lineng's blog
03-01 457
题目 输入一串数字，给你 M 个询问，每次询问就给你两个数字 X,Y，要求你说出 X 到 Y 这段区间内的最大数。 输入格式 第一行两个整数 N,M 表示数字的个数和要询问的次数； 接下来一行为 N 个数； 接下来 M 行，每行都有两个整数 X,Y。 输出格式 输出共 M 行，每行输出一个数。 数据范围 1≤N≤1051≤N≤10^51≤N≤105, 1≤M≤1061≤M≤10^61≤M≤106, 1≤X≤Y≤N1≤X≤Y≤N1≤X≤Y≤N, 数列中的数字均不超过231−12^{31}−1231−1 Sol
No Time to Dry【ST表】【树状数组】
//(^ o ^)//
08-17 133
>Link luogu P7416 jzoj7232 >Description >解题思路 赛时打的暴力单调栈20分QAQ 观察一个区间[l,r][l,r][l,r]，最多的步骤为r−l+1r-l+1r−l+1，如果有cntcntcnt步可以不做，最少的步骤就为r−l+1−cntr-l+1-cntr−l+1−cnt 怎么样的步数可以合并？ 单调栈的思路，如果当前的数字aia_iai​，与上一次出现的位置之间，没有小于aia_iai​的数字，那当前的位置的步骤就可以与上一次出现的位置
数据结构：线段树及ST算法比较
aotao4494的博客
11-01 370
ST算法是一种高效的计算区间最值的方法。 他的思想是将询问区间分解成两个最长的二次幂的长度的区间并集的形式。 所以与线段树不同，这种区间分解其实存在相交的分解。 因此ST算法能维护的只是一些简单的信息，比如区间最值或者区间gcd问题 ST算法的优势： 实现简单（qwq为什么我觉得线段树更好理解） 效率比线段树高 线段树的优势： 可以更好地维护动态的信息，而ST...
ACM模板——线段树&树状数组&ST表
weixin_30415801的博客
04-05 118
int bit[maxn],n; void init() { n = maxn; memset(bit,0,sizeof(bit)); } int sum(int i) { int s = 0; while(i>0) { s += bit[i]; i -= i&-i; ...
专题 | ST表、树状数组、线段树、分块等
最新发布
hzx616599的博客
03-10 1027
分裂就是分裂一个 node，变成两个小的 node，以保证每个 node 的大小都接近。
bzoj 3521/5083 [Poi2014]Salad Bar/普及 （卡常）st表+二分+树状数组
I'm Growing
10-26 613
Description有一个长度为n的字符串，每一位只会是p或j。你需要取出一个子串S（从左到右或从右到左一个一个取出），使得 不管是从左往右还是从右往左取，都保证每时每刻已取出的p的个数不小于j的个数。你需要最大化|S|。 Input第一行一个整数n，接下来一个长度为n的只含有p,j的字符串 N<=10^6 Output输出S的最大长度 Sample Input6jpjppj Samp
线段树详解：数据结构与RMQ问题
“线段树资料，讲解线段树在数据结构和算法中的应用。” 线段树是一种高效的数据结构，常用于解决区间查询和修改的问题。它具有完全二叉树的特性，每个节点代表一个区间，从根节点表示的整体区间到叶子节点表示的...
数据结构：RMQ与LCA算法解析及线段树应用
LCA（最近公共祖先）问题是在树形结构中寻找两个指定节点的最深共同祖先。这个问题在图论和数据结构中有着广泛的应用。 线段树算法是一种高效的数据结构，它可以用来解决包括RMQ在内的多种区间问题。线段树的本质是...
线段树 + 树状数组 + ST表 模板
jesHrz的博客
08-07 239
线段树 区间修改+区间求和 logN const int N = 1e5 + 5; int a[N]; namespace SEG { struct tag { long long sum; long long lazy; } c[N &lt;&lt; 2]; void pushDown(int k, int l, int r)...
POJ3264 RMQ的ST算法
01-23 488
RMQ问题为求区间最值的问题   线段树可以在O(logN)的时间复杂度内完成询问操作。 但是ST算法可以在常数时间内完成询问操作   ST算法：基于动态规划求区间最值的算法。 分为预处理和查询两部分 预处理：定义 F[i][j] 为从 i开始到 i+2^j-1 区间内的最值 , 我们可以讲这段2^j的区间分成两部分长度都为2^(j-1)的相同区间 区间1 为  i.....i+2
实用数据结构小结（树状数组、st表、bitset）
m0_54646258的博客
07-07 310
简要小结 树状数组、st表、bitset 这些在具体场景求解高效、代码简短的小数据结构
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值