rmq模板（poj3264验证）

搜索不收藏，线段树方法不收藏，只收藏st算法，这个算法和线段树差不多，建树和建dp数组的时间复杂度都是n*log（n）

不同的地方就是查找，线段树查找的时间复杂度是n*log（n），st算法是o（1）

所以线段树的总时间复杂度是（设m为查询次数）： m*n*log(n)

而st算法是n*log(n)，要是查询多了。。。。线段树果断悲剧。。。。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
//#define M 100010
#define MAX 50010
//#define MAXM 500

int top,n,x,y;
int f[MAX*2][30],f2[MAX*2][30];
int seq[MAX*2];

int mmin(int x,int y){
    if(x>y) return y;
    else return x;
}

int mmax(int x,int y){
    if(x>y) return x;
    else return y;
}

int rmq(){
    int i,j,k;
    int m=(int)(log(top-1+0.0)/log(2.0));
    for(i=1;i<top;++i){
        f[i][0]=seq[i];
        f2[i][0]=seq[i];
    }
    for(j=1;j<=m;++j){
        k=1<<(j-1);
        for(i=1;i+(1<<(j-1))<top;++i){
        f[i][j]=mmin(f[i][j-1],f[i+k][j-1]);
        f2[i][j]=mmax(f2[i][j-1],f2[i+k][j-1]);
        }
        //cout<<"i= "<<i<<",j= "<<j<<"  "<<f[i][j]<<endl;
    }
    return 0;
}

int query(int l,int r){
    if(l>r)swap(l,r);
    int k=(int)(log(r-l+1.0)/log(2.0));
    //cout<<mmax(f2[l][k],f2[r-(1<<k)+1][k])<<" - "<<mmin(f[l][k],f[r-(1<<k)+1][k])<<" = ";
    return mmax(f2[l][k],f2[r-(1<<k)+1][k])-mmin(f[l][k],f[r-(1<<k)+1][k]);
}
int main(){
    scanf("%d%d",&top,&n);
        top++;
        for(int i=1;i<top;i++)
            scanf("%d",&seq[i]);
        rmq();
        /*for(int i=1;i<top;i++)
            cout<<i<<": "<<f[i][0]<<endl;*/

        //scanf("%d",&n);
        while(n--){
            scanf("%d%d",&x,&y);
            cout<<query(x,y)<<endl;
        }

    return 0;
}