经典RMQ（poj3264）

题意：随机访问任意区间的极值；

ST算法：DP预处理复杂度O(nlogn)；以后每次询问复杂度为O(1);

(F[i,j]表示从i位置开始，往右2^j个元素中的极值。动规方程F[i,j]=max(F[i,j-1],F[i+2^(j-1),j-1])，当i+2^(j-1)大于序列左端时，则F[i,j]=(F[i,j-1]；

在访问（a，b）区间时，取j使得2^j刚好小于等于（b-a）；然后求得max（F[i,j],F[b-2^j+1,j])即可；

#include <iostream>

#include <stdio.h>

using namespace std;

int maxtool[50010][20];

int mintool[50010][20];

int rem[50010];

int pow(int k)

{

    if(k==0)

    return 1;

    if(k==1)

    return 2;

    if(k==2)

    return 4;

    int num=pow(k/2);

    if(k%2)

    return num*num*2;

    return num*num;

}

int finds(int k)

{

    int i=0;

   while(pow(i)<=k)

    {

       i++;

    }

    return i-1;

}

int main()

{

    int n,m;

   scanf("%d%d",&n,&m);

    for(inti=0;i<n;i++)

   scanf("%d",&rem[i]);

    for(intj=0;j<n;j++)

   maxtool[j][0]=mintool[j][0]=rem[j];

    for(int i=1;;i++)

    {

       if(pow(i)>=n)

       break;

       for(int j=0;j<n;j++)

       {

          if(j+pow(i-1)>=n)

           {

          maxtool[j][i]=maxtool[j][i-1];

          mintool[j][i]=mintool[j][i-1];

           }

           else

           {

          mintool[j][i]=min(mintool[j][i-1],mintool[j+pow(i-1)][i-1]);

          maxtool[j][i]=max(maxtool[j][i-1],maxtool[j+pow(i-1)][i-1]);

           }

       }

    }

    for(inti=0;i<m;i++)

    {

       int a,b;

      scanf("%d%d",&a,&b);

       a--;

       b--;

       int h=finds(b-a);

      cout<<max(maxtool[a][h],maxtool[b-pow(h)+1][h])-min(mintool[a][h],mintool[b-pow(h)+1][h])<<endl;

    }

    return 0;

}