-
题目描述:
-
给定一个数字序列,查询任意给定区间内数字的最小值。
-
输入:
-
输入包含多组测试用例,每组测试用例的开头为一个整数n(1<=n<=100000),代表数字序列的长度。
接下去一行给出n个数字,代表数字序列。数字在int范围内。
下一行为一个整数t(1<=t<=10000),代表查询的次数。
最后t行,每行给出一个查询,由两个整数表示l、r(1<=l<=r<=n)。
-
输出:
-
对于每个查询,输出区间[l,r]内的最小值。
-
样例输入:
-
5 3 2 1 4 3 3 1 3 2 4 4 5
-
样例输出:
-
1 1 3
RMQ (Range Minimum/Maximum Query)问题是指:对于长度为n的数列A,回答若干询问RMQ(A,i,j)(i,j<=n),返回数列A中下标在i,j里的最小(大)值,也就是说,RMQ问题是指求区间最值的问题。
主要方法及复杂度如下:
1、朴素(即搜索),O(n)-O(qn) online。
2、
线段树,O(n)-O(qlogn) online。
3、ST(实质是
动态规划),O(nlogn)-O(q) 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(q) online。
首先根据原
数列,建立
笛卡尔树,从而将问题在线性时间内规约为LCA问题。LCA问题可以在线性时间内规约为约束RMQ,也就是数列中任意两个相邻的数的差都是+1或-1的RMQ问题。约束RMQ有O(n)-O(1)的在线解法,故整个算法的时间复杂度为O(n)-O(1)。
下面是方法3的实现:
1 #include <climits> 2 #include <cmath> 3 #include <cstdio> 4 #include <iostream> 5 using namespace std; 6 7 int n, t; 8 int l, r; 9 int dp[100001][32]; 10 11 void init() { 12 for (int j = 1; (1<<j) <= n; ++j) { 13 for (int i = 1; i + (1<<j) - 1 <= n; ++i) { 14 dp[i][j] = min(dp[i][j-1], dp[i+(1<<(j-1))][j-1]); 15 } 16 } 17 } 18 19 void getRes() { 20 int k = (int) (log((double)(r-l+1))/log(2.0)); 21 int res = min(dp[l][k], dp[r-(1<<k)+1][k]); 22 printf("%d\n", res); 23 } 24 25 int main() { 26 while (scanf("%d", &n) != EOF) { 27 for (int i = 1; i <= n; ++i) { 28 scanf("%d", &dp[i][0]); 29 } 30 init(); 31 scanf("%d", &t); 32 for (int i= 0; i < t; ++i) { 33 scanf("%d %d", &l, &r); 34 getRes(); 35 } 36 } 37 return 0; 38 } 39 /************************************************************** 40 Problem: 1544 41 User: hupo250 42 Language: C++ 43 Result: Accepted 44 Time:320 ms 45 Memory:14068 kb 46 ****************************************************************/