poj 3264 Balanced Lineup

RMQ模板题，用ST算法

 

//DP预处理

//dp[i][j] 表示从下标i开始，长度为2^j的最大值

//状态转移方程 dp[i][j] = max{ dp[i][j-1] , dp[i+2^(j-1)][j-1] }

//也就是一个长度为2^j的区间，二分为两个2^(j-1)的长度

//对于查询[a,b]以内的最大值,先求出区间长度LEN = b-a+1

//查询结果为 res = max{dp[a][k] , dp[b-2^k+1][k]} ， 关键是k，是什么，怎么计算

//当满足 2^k >= LEN/2  ， 且k最小时，这个k就是我们要的值

//k可以用计算公式一步算得 : k = (int) (log((double)(b-a+1))/log(2.0))

//另外一个问题，对于构建dp值的时候,dp[i][j]，这个j最大去到多少？

//道理是一样的，对于总区间，长度为L，同样是 2^j >= L/2 , 取最小的j就是我们要的值

//所以也可以用公式一步计算出来  j = (int) (log((double)L) / log(2.0))



#include <iostream>

#include <cstdio>

#include <cstring>

#include <cmath>

using namespace std;

#define N 50010

#define M 20



int n,m,a[N],dp[N][M],__dp[N][M];



void TEST()

{

    cout << "最大值" << endl;

    int K = (int) (log((double)n) / log(2.0));

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

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

            printf("dp[%d,%d] = %d\n",i,j,dp[i][j]);

    cout << "最小值" << endl;

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

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

            printf("__dp[%d,%d] = %d\n",i,j,__dp[i][j]);



}



void ST()

{

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

        dp[i][0] = __dp[i][0] = a[i];



    int K = (int) (log((double)n) / log(2.0));

    for(int j=1; j<=K; j++)

        for(int i=1; i+(1<<j)-1 <= n; i++)

        {

            int k = i+(1<<(j-1));

            dp[i][j] = max( dp[i][j-1] , dp[k][j-1] );

            __dp[i][j] = min( __dp[i][j-1] , __dp[k][j-1]);

        }

}



int RMQ(int x ,int y)

{

    int K = (int) (log((double)(y-x+1)) / log(2.0));

    int k = y-(1<<K)+1;

    int Max = max( dp[x][K] , dp[k][K]);

    int Min = min( __dp[x][K] , __dp[k][K]);

    return Max - Min;

}



int main()

{

    while(cin >> n >> m)

    {

        for(int i=1; i<=n ;i++) cin >> a[i];

        ST();

        //TEST();

        while(m--)

        {

            int x,y;

            cin >> x >> y;

            int res = RMQ(x,y);

            cout << res << endl;

        }

    }

    return 0;

}