题目链接: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)); } }