pku 3264 Balanced Lineup(RMQ)

首先依次给定一条线段上每个点的高度。然后执行查询操作，要求输出在范围[i,j]内最高点和最低点的高度之差。

 

这是一个很典型的RMQ问题，关于这个问题有个很出名的DP解答方式。

ST（动态规划） O(nlogn)-O(1) offline

通过O(nlogn)的预计算，可以达到O(1)的查询。

 

关于这个算法，在百度百科里面摘抄一点:

ST算法（Sparse Table），以求最大值为例，设d[i,j]表示[i,i+2^j-1]这个区间内的最大值，那么在询问到[a,b]区间的最大值时答案就是 max(d[a,k], d[b-2^k+1,k])，其中k是满足2^k<=b-a的最大的k,即k=[ln(b-a+1)/ln(2)]。

（这里d[a,k]和d[b-2^k+1,k])可能存在重合区段，但这不影响最后的结果）

 

这个算法虽然说是动态规划，但其实就是二分的思想，先计算出包含2^i个数的区间上的最大值（或最小值），然后再计算包含2^(i+1)个数的区间上的最大值（或最小值）。最后查询的时候就找到两个区间，都处在那个大的区间内而又完全覆盖了那个区间。

 

同时这个算法和线段树很像，优点是因为没了树的结构，每一区间的关键字可以在O(1)的时间内查询到，缺点是不能动态的插入，如果需要插入，需要重新运行DP的预计算过程。

 

#include "math.h" #include <iostream> using namespace std; const int maxn=50005; const int maxi=16; int height[maxn]; int min_height[maxn][maxi]; int max_height[maxn][maxi]; #define MinN(x,y) (x<y?x:y) #define MaxN(x,y) (x>y?x:y) void DP(int n) { for(int i=1;i<=n;i++) { min_height[i][0]=max_height[i][0]=height[i]; } int floor=log((double)n)/log(2.0)+1; for(int i=1;i<floor;i++) { for(int j=1;j<=n;j++) { int mid=j+(1<<(i-1)); if(mid<=n) { max_height[j][i]=MaxN(max_height[j][i-1],max_height[mid][i-1]); min_height[j][i]=MinN(min_height[j][i-1],min_height[mid][i-1]); } else { max_height[j][i]=max_height[j][i-1]; min_height[j][i]=min_height[j][i-1]; } } } } int Query(int l,int r) { int m=log((double)(r-l+1))/log(2.0); int front=r-(1<<m)+1; int m_max=MaxN(max_height[l][m],max_height[front][m]); int m_min=MinN(min_height[l][m],min_height[front][m]); return m_max-m_min; } int main() { int n,q; int l,r; scanf("%d%d",&n,&q); for(int i=1;i<=n;i++) { scanf("%d",&height[i]); } DP(n); for(int i=0;i<q;i++) { scanf("%d%d",&l,&r); printf("%d/n",Query(l,r)); } return 0; }