Poj 3264 Balanced Lineup (数据结构_RMQ)

题目链接:http://poj.org/problem?id=3264


题目大意给定一个序列,n个数,给出m个询问,每次询问[li,ri]区间的最大值与最小值的差。


解题思路:一直想花点时间补下数据结构方面的知识,但因为状态dp和树形dp较难,花了大把时间。想放弃又害怕痛苦!到现在也没能搞定这两个家伙,今天早晨终于归零,开始学RMQ.

    这题应该是有两种做法,建两可线段树复杂度o(nlogn),o(logn)查询,总复杂度O((n+m)logn,但这应该不会超时,没去敲。我用的是St算法,O(nlogn)初始化一个dp[MAX][MIN]数组,O(1)查询,St算法的优越之处就在于这个查询速度。

    那么它是怎么实现的呢?其实就是DP,运用倍增思想,把区间压缩。dp[i][j]表示从位置i开始(包括i自己)长度为1<<j的区间的最值,那么dp[i][0] = arr[i],这是显然的,那么dp[i][j]该怎么转移呢?dp[i][j] = max(dp[i][j-1],dp[i+(1<<(j-1))][j-1]),把长度为1<<j的区间分成两段长度为1<<(j-1)的区间,取其中的最大值。长度长的从长度短的转移过来,所以在编码中要先写j的循环再写i循环。初始完成后就坐等高效查询。当需要查询[l,r]之间最值的时候,先计算出区间长度len,然后求k = log2len(即log(r-l+1)/log(2)),最后max(dp[l][k],dp[r-(1<<k)+1][k])即是答案。为什么可以用dp[l][k]和dp[r-(1<<k)+1][k]表示一个区间?他们表示的是两个至少一个交点的区间,如果大家不是很明白,可以想下查询【1,8】或【1,10】的情况,前者8个交点,后者6个交点。咳,口齿不清,表达得不清楚,一句话概括下--那是必然的。


测试数据:

6 3
1
7
3
4
2
5
1 5
4 6
2 2

代码:
#include <stdio.h>
#include <string.h>
#include <math.h>
#define MIN 20
#define MAX 51000
#define min(a,b) (a)<(b)?(a):(b)
#define max(a,b) (a)>(b)?(a):(b)


int n,m,arr[51000];

//模板
struct RMQ {

    int n,mmax[MAX][MIN],mmin[MAX][MIN];
    void Create(int flag);
    int Query(int x,int y,int flag);
}rmq;
void RMQ::Create(int flag) {
    //flag = 0建立最大rmq,flag = 1最小
    int i,j,k;

    if (flag == 0) {

        for (i = 1; i <= n; ++i)
            mmax[i][0] = arr[i];
        for (j = 1; (1<<j) <= n; ++j)
            for (i = 1; i+(1<<j)-1 <= n; ++i)
                mmax[i][j] = max(mmax[i][j-1],mmax[i+(1<<(j-1))][j-1]);
    }
    else {

        for (i = 1; i <= n; ++i)
            mmin[i][0] = arr[i];
        for (j = 1; (1<<j) <= n; ++j)
            for (i = 1; i+(1<<j)-1 <= n; ++i)
                mmin[i][j] = min(mmin[i][j-1],mmin[i+(1<<(j-1))][j-1]);
    }
}
int RMQ::Query(int l,int r,int flag){
    //flag = 0查询最大rmq,flag = 1最小
    int k = (int)(log(r-l+1.0)/log(2.0));
    if (flag == 0)
        return max(mmax[l][k],mmax[r-(1<<k)+1][k]);
    else
        return min(mmin[l][k],mmin[r-(1<<k)+1][k]);
}


int main()
{
    int i,j,k;


    scanf("%d%d",&n,&m);
    for (i = 1; i <= n; ++i)
        scanf("%d",&arr[i]);


    rmq.n = n;
    rmq.Create(0);
    rmq.Create(1);
    for (i = 1; i <= m; ++i) {

        scanf("%d%d",&j,&k);
        printf("%d\n",rmq.Query(j,k,0) - rmq.Query(j,k,1));
    }
}





  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值