实现RMQ的两种常用方法

RMQ

　　RMQ(Range Maximum/Minimum Question)是指区间最值问题，在OI中较为常见，一般可以用ST表和线段树实现。

　　ST表是基于倍增思想的一种打表方法，在确定区间范围和所有的值后利用倍增预处理出$2^k$长度的区间内的最值，然后$O(1)$查询。优点是查询快且操作简便，缺点是不能进行动态操作，只支持静态查询。

　　Code：（POJ模板题）

//It is made by HolseLee on 23rd July 2018
//POJ 3264
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<iomanip>
using namespace std;
const int N=5e4+;
int n,m,a[N];
int mx[N][],mi[N][];
inline int read()
{
  char ch=getchar();int num=;bool flag=false;
  while(ch<''||ch>''){if(ch=='-')flag=true;ch=getchar();}
  while(ch>=''&&ch<=''){num=num*+ch-'';ch=getchar();}
  return flag?-num:num;
}
void ready()
{
  for(int j=;j<=;j++)
    for(int i=;i<=n;i++)
      if(i+(<<(j-))<=n){
    mx[i][j]=max(mx[i][j-],mx[i+(<<(j-))][j-]);
    mi[i][j]=min(mi[i][j-],mi[i+(<<(j-))][j-]);}
}
inline int quary(int l,int r)
{
  int maxx=-,minn=;
  int k=(int)(log((double)(r-l+))/log(2.0));
  maxx=max(mx[l][k],mx[r-(<<k)+][k]);
  minn=min(mi[l][k],mi[r-(<<k)+][k]);
  return maxx-minn;
}
int main()
{
  n=read();m=read();
  memset(mx,-,sizeof(mx));
  memset(mi,0x7f,sizeof(mi));
  for(int i=;i<=n;i++){
    a[i]=read();
    mi[i][]=mx[i][]=a[i];}
  ready();int x,y;
  for(int i=;i<=m;i++){
    x=read();y=read();
    printf("%d\n",quary(x,y));}
  return ;
}

　　线段树就不用多说了，除了可以同时维护最大和最小值意外，还能维护更多信息，操作同样也方便，而且还支持动态操作。

　　Code：

//It is made by HolseLee on 23rd July 2018
//POJ 3264
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<iomanip>
using namespace std;
const int N=5e4+;
int n,m,a[N];
struct Node{
    int mx,mi;
    Node(int xx=,int yy=)
    {mx=xx;mi=yy;}
}seg[N<<];
inline int read()
{
    char ch=getchar();int num=;bool flag=false;
    while(ch<''||ch>''){if(ch=='-')flag=true;ch=getchar();}
    while(ch>=''&&ch<=''){num=num*+ch-'';ch=getchar();}
    return flag?-num:num;
}
inline void update(int rt)
{
    seg[rt].mx=max(seg[rt<<].mx,seg[rt<<|].mx);
    seg[rt].mi=min(seg[rt<<].mi,seg[rt<<|].mi);
}
inline void build(int l,int r,int rt)
{
    if(l>r)return;
    if(l==r){
        seg[rt].mx=seg[rt].mi=a[l];return;}
    int mid=(l+r)>>;
    build(l,mid,rt<<);build(mid+,r,rt<<|);
    update(rt);
}
inline Node quary(int l,int r,int rt,int L,int R)
{
    Node ret(-,);
    if(l>R||r<L)return ret;
    if(L<=l&&r<=R){return seg[rt];}
    int mid=(l+r)>>;
    Node lc(-,);
    Node rc(-,);
    if(L<=mid)lc=quary(l,mid,rt<<,L,R);
    if(R>mid)rc=quary(mid+,r,rt<<|,L,R);
    ret.mx=max(lc.mx,rc.mx);
    ret.mi=min(lc.mi,rc.mi);
    return ret;
}
inline int get(int x,int y)
{
    int maxx=quary(,n,,x,y).mx;
    int minn=quary(,n,,x,y).mi;
    return maxx-minn;
}
int main()
{
    n=read();m=read();
    for(int i=;i<=n;i++)
    a[i]=read();
    build(,n,);int x,y;
    for(int i=;i<=m;i++){
        x=read();y=read();
        printf("%d\n",get(x,y));}
    return ;
}