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

一、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;
}

 

 

四、主席树

 

五、莫队