RMQ（st算法离线算法）求区间最大值/最小值

定义dp[i][j]为起始位置为i，终点为i+2^j-1的区间的最大值/最小值；
转移方程 dp[i][j]=max(dp[i][j-1],dp[i+(1<<j-1)][j-1])
具体代码：

void st(int n){
    for(int i=1;i<=n;++i) dp[i][0]=p[i].w;
    int len=log(n)/log(2)+1;//区间长度为n的二进制最大数，满足2^len>n
    for(int i=1;i<len;++i){
        int s=1<<(i-1);
        for(int j=1;j<=n-(1<<i)+1;++j){
            dp[j][i]=max(dp[j][i-1],dp[j+s][i-1]);
        }
    }
}
int query(int l,int r){
    int len=log(r-l+1)/log(2);//满足2^len<r-l+1<=2^(len+1)
    return max(dp[l][len],dp[r-(1<<len)+1][len]);l+
    //答案一定在[l,l+2^len-1]和[r-2^len+1,r]，这两个区间是相邻或者相交的，因为长度相加为2^(len+1)>=r-l+1>2^len。
}

练习：POJ3368
思路：
因为序列是单调的，所以我们可以离散化，对于连续相同数值的区间缩为一个点，点包含l，r和区间长度；
对于一个区间如果满足a[l]==a[r] ，就输出r-l+1；（区间的数值都相同，一个区间）
否则就在左右两边的区间跟中间的区间找最大值。（两个或多个不同区间的情况）
代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn=1e5+11;
const int mod=1e9+7;
const double eps=0.000001;
typedef long long ll;
typedef long double lb;
typedef double db;
struct D{
    int l,r,w;
}p[maxn];
int dp[maxn][25],pos[maxn],a[maxn];
void st(int n){
    for(int i=1;i<=n;++i) dp[i][0]=p[i].w;
    int len=log(n)/log(2)+1;
    for(int i=1;i<len;++i){
        int s=1<<(i-1);
        for(int j=1;j<=n-(1<<i)+1;++j){
            dp[j][i]=max(dp[j][i-1],dp[j+s][i-1]);
        }
    }
}
int query(int l,int r){
    if(r<l) return 0;
    int len=log(r-l+1)/log(2);
    return max(dp[l][len],dp[r-(1<<len)+1][len]);
}
int main(){
    int n,q,id,l,r;
    while(scanf("%d",&n)==1&&n){
        scanf("%d",&q);
        id=1;
        p[id].l=pos[1]=1;
        for(int i=1;i<=n;++i) scanf("%d",&a[i]);
        for(int i=2;i<=n;++i){
            if(a[i]!=a[i-1]){
                p[id].r=i-1;
                p[id].w=i-p[id].l;
                p[++id].l=i;
            }
            pos[i]=id;
        }
        p[id].r=n,p[id].w=n-p[id].l+1;
        st(id);
        while(q--){
            scanf("%d%d",&l,&r);
            if(pos[l]==pos[r]) printf("%d\n",r-l+1);
            else{
                int a1=p[pos[l]].r-l+1;
                int a2=r-p[pos[r]].l+1;
                int a3=query(pos[l]+1,pos[r]-1);
                int ans=max(max(a1,a2),a3);
                printf("%d\n",ans);
            }
        }
    }
    system("pause");
}