RMQ你让我懂得了dp的强大

 最近碰到一个线段树的题目，发现discuss中说可以用RMQ做，以前也了解过，但是没有深入地去理解，今天查了一些资料，好好地整理一下。

表面上RMQ，看着还挺高级，实际上就是一个简单的dp，只不过这个dp有点高端，特别是状态方程，形式上和原来没什么两样，但是内容上却差别很大。好了废话少说，看看例子。

 

1、什么是RMQ呢？

RMQ (Range Minimum/Maximum Query)问题是指：对于长度为n的数列A，回答若干询问RMQ(A,i,j)(i,j<=n)，返回数列A中下标在i,j里的最小(大）值，也就是说，RMQ问题是指求区间最值的问题。所以以后要是求区间的最大值，并且会频繁地访问的时候就可以用这种算法。

其实很多时候有些题目表面上看不能用RMQ做，但是可以做一下转化。比如上一篇说的那一个poj3368，表面上看也是区间的最大值，但是直接算却不好算，不能直接用RMQ，我们将相同的数合为一个数，这个数代表相同数字的数目，这不就成了区间最大值了吗？

2、复杂度对比

       1、朴素（即搜索），O(n)-O(qn) online。

　　2、线段树，O(n)-O(qlogn) online。

　　3、ST（实质是动态规划），O(nlogn)-O(1) online。

　　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+1(即长度)的最大的k,即k=[ln(b-a+1)/ln(2)]。

　　d的求法可以用动态规划，d[i, j]=max(d[i, j-1],d[i+2^(j-1), j-1])。

　　4、RMQ标准算法：先规约成LCA（Lowest Common Ancestor），再规约成约束RMQ，O(n)-O(1) online。

　　首先根据原数列，建立笛卡尔树，从而将问题在线性时间内规约为LCA问题。LCA问题可以在线性时间内规约为约束RMQ，也就是数列中任意两个相邻的数的差都是+1或-1的RMQ问题。约束RMQ有O(n)-O(1)的在线解法，故整个算法的时间复杂度为O(n)-O(1)。

引用百度的说明。

3、举例子，揭开神秘面纱

设a是要求区间最值的数列，f[i,j]表示从第i个数起连续2^j个数中的最大值。例如数列3 2 4 5 6 8 1 2 9 7 ,f[1,0]表示第1个数起，长度为2^0=1的最大值，其实就是3这个数。f[1,2]=5，f[1,3]=8，f[2,0]=2，f[2,1]=4……从这里可以看出f[i,0]其实就等于a。这样，DP的状态、初值都已经有了，剩下的就是状态转移方程。我们把f[i,j]（j≥1）平均分成两段（因为j≥1时，f[i,j]一定是偶数个数字），从i到i+2^(j-1)-1为一段，i+2^(j-1)到i+2^j-1为一段（长度都为2^（j-1））。用上例说明，当i=1，j=3时就是3,2,4,5 和6,8,1,2这两段。f就是这两段的最大值中的最大值。于是我们得到了动规方程F[i,j]=max(F[i,j-1],F[i+2^(j-1),j-1])。

　　接下来是得出最值，也许你想不到计算出f有什么用处，一般要想想计算max还是要O(logn)，甚至O(n)。但有一个很好的办法，做到了O(1)。还是分开来。如在上例中我们要求区间[2,8]的最大值，就要把它分成[2,5]和[5,8]两个区间，因为这两个区间的最大值我们可以直接由f[2,2]和f[5,2]得到。扩展到一般情况，就是把区间[l,r]分成两个长度为2^n的区间（保证有f对应）。直接给出表达式：

　　k:=trunc(ln(r-l+1)/ln(2));

　　ans:=max(F[l,k],F[r-2^k+1,k]);

　　这样就计算了从l开始，长度为2^k的区间和从r-2^k+1开始长度为2^k的区间的最大值（表达式比较烦琐，细节问题如加1减1需要仔细考虑），二者中的较大者就是整个区间[l,r]上的最大值。

4、总结

RMQ和LCA在算法上可以相互转换，这一块，以后碰见再总结吧。