实用数据结构---RMQ问题

RMQ问题就是区间最小值问题,这是一个非常经典的题,

由他引申出来的也是不计其数最多的是给出一个区间,然后输入多组区间端点,求输入区间的最小值。

每次用循环来计算一个最小值显然不够快,怎么办呢?

实践中最常用的是Tarjan的 Sparse-Table算法,它的预处理时间是O(nlogn),但是查询只需要O(1),而且常数很小。

它的思想很简单,就是递推+二分的思想。我们先定义一个二维数组d[i][j],

代表从 i 开始,长度为2的j次幂的区间中的最小值

然后写出递推公式
d[i][j]
= min( d[i][j
– 1], d[i
+ (1 << (j - 1))][j
- 1]),

其中1
<< (j - 1) 代表2的 j – 1次幂查询操作很简单,

我们要先求出一个数k 让2的k次幂小于等于R-L+1,

那么d[L][k]与d[R-(1
<< k) + 1][k]就覆盖了整个L到R的区间,读者可自行验证。

现在就可以写出完整的代码了

#include<stdio.h>
#include<algorithm>
using namespace std;
int a[100];
int dp[100][10];
int query(int L, int R){
	int k = 0;
	while (1 << (k + 1) <= R - L + 1)
		k++;
	return min(dp[L][k], dp[R - (1 << k) + 1][k]);
}
int main(){
	int n, i, j;
	scanf("%d", &n);
	for (i = 1; i <= n; i++){
		scanf_s("%d", &a[i]);
		dp[i][0] = a[i];
	}
	for (i = 1; (1 << i) <= n; i++){
		for (j = 1; j + (1 << (i - 1)) <= n; j++){
			dp[j][i] = min(dp[j][i - 1], dp[j + (1 << (i - 1))][i - 1]);
		}
	}
	int L, R;
	scanf("%d%d", &L, &R);
	printf("%d\n", query(L, R));
	return 0;
}

附上一个相关练习[南阳oj上的一道题
很好](http://acm.nyist.net/JudgeOnline/problem.php?pid=119)

时间: 2024-11-01 15:02:11

实用数据结构---RMQ问题的相关文章

C++实用数据结构:二叉索引树

看下面这个问题(动态连续和查询): 有一个数组A(长度为n),要求进行两种操作: add(i,x):让Ai增大x: query(a,b):询问Aa+Aa+1+...+Ab的和: 若进行模拟,则每次query操作的最坏的时间复杂度为O(n),在n较大时速度较慢.用前缀和也不能提高效率(每次add操作最坏为O(n)).有一种数据结构,可以在O(n)时间里初始化,用O(logn)的速度执行add操作或查询前缀和,从而执行query操作. 首先,我们来介绍“lowbit”.对于一个数x,lowbit(x

【暑假】[实用数据结构]UVa11235 Frequent values

UVa 11235 Frequent values Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 11241   Accepted: 4110 Description You are given a sequence of n integers a1 , a2 , ... , an in non-decreasing order. In addition to that, you are given several qu

实用数据结构---并查集

题目描述 Description 若某个家族人员过于庞大,要判断两个是否是亲戚,确实还很不容易,现在给出某个亲戚关系图,求任意给出的两个人是否具有亲戚关系. 规定:x和y是亲戚,y和z是亲戚,那么x和z也是亲戚.如果x,y是亲戚,那么x的亲戚都是y的亲戚,y的亲戚也都是x的亲戚. 输入描述 Input Description 第一行:三个整数n,m,p,(n<=5000,m<=5000,p<=5000),分别表示有n个人,m个亲戚关系,询问p对亲戚关系. 以下m行:每行两个数Mi,Mj,

实用数据结构---树状数组(二叉索引树)

树状数组适用于动态连续和查询问题,就是给定一个区间, 查询某一段的和或者修改某一位置的值. 关于树状数组的结构请去百度百科,否则将看不懂下面内容 我们看这个题 士兵杀敌(二) 时间限制:1000 ms  |  内存限制:65535 KB 难度:5 描述 南将军手下有N个士兵,分别编号1到N,这些士兵的杀敌数都是已知的. 小工是南将军手下的军师,南将军经常想知道第m号到第n号士兵的总杀敌数,请你帮助小工来回答南将军吧. 南将军的某次询问之后士兵i可能又杀敌q人,之后南将军再询问的时候,需要考虑到新

实用数据结构---图的操作和算法

图算法的源代码.包含大量的注释,和最小生成树.最短路径.邻接表图深度广度优先搜索,邻接矩阵图深度广度优先搜索,欢迎借鉴 #include<stdio.h> #include<stdlib.h> #define MAXVEX 20 #define INFINITY 65535 typedef char vertexType; typedef int edgeType; typedef int Boolean; typedef int Pathmatirx[MAXVEX]; //用于存

【暑假】[实用数据结构]UVAlive

题目:   Dominating Patterns Time Limit: 3000MS   Memory Limit: Unknown   64bit IO Format: %lld & %llu Submit Status Description The archaeologists are going to decipher a very mysterious ``language". Now, they know many language patterns; each patt

【暑假】[实用数据结构]UVAlive 3026 Period

UVAlive 3026 Period 题目: Period Time Limit: 3000MS   Memory Limit: Unknown   64bit IO Format: %lld & %llu Submit Status Description For each prefix of a given string S with N characters (each character has an ASCII code between 97 and 126, inclusive),

【暑假】[实用数据结构]KMP

KMP算法 KMP算法是字符串匹配算法,可以在O(n)的时间完成,算法包含两部分,分别是:构造适配函数与两串匹配. 失配边的使用大大提高了算法效率,可以理解为已经成功匹配的字符不在重新匹配,因为我们已经知道它是什么,对应到算法中 匹配失败后应该在最大前缀之后继续匹配,因为某后缀已与最大前缀匹配成功而不用重新比较. 以下为代码实现: 1 const int maxn = 1000 + 5; 2 3 void getFail(char* P,int* f){ //构造失配边 4 int n=strl

【暑假】[实用数据结构]UVAlive 3942 Remember the Word

UVAlive 3942 Remember the Word 题目: Remember the Word Time Limit: 3000MS   Memory Limit: Unknown   64bit IO Format: %lld & %llu Submit Status Description Neal is very curious about combinatorial problems, and now here comes a problem about words. Know