hihoCoder #1068 : RMQ-ST算法

最新推荐文章于 2018-08-08 22:10:03 发布

fushenliao

最新推荐文章于 2018-08-08 22:10:03 发布

阅读量390

点赞数

分类专栏：动态规划

本文链接：https://blog.csdn.net/Night_13/article/details/55517176

版权

动态规划专栏收录该内容

15 篇文章 0 订阅

订阅专栏

#1068 : RMQ-ST算法

时间限制: 10000ms

单点时限: 1000ms

内存限制: 256MB

描述

小Hi和小Ho在美国旅行了相当长的一段时间之后，终于准备要回国啦！而在回国之前，他们准备去超市采购一些当地特产——比如汉堡（大雾）之类的回国。

但等到了超市之后，小Hi和小Ho发现者超市拥有的商品种类实在太多了——他们实在看不过来了！于是小Hi决定向小Ho委派一个任务：假设整个货架上从左到右拜访了N种商品，并且依次标号为1到N，每次小Hi都给出一段区间[L, R]，小Ho要做的是选出标号在这个区间内的所有商品重量最轻的一种，并且告诉小Hi这个商品的重量，于是他们就可以毫不费劲的买上一大堆东西了——多么可悲的选择困难症患者。

（虽然说每次给出的区间仍然要小Hi来进行决定——但是小Hi最终机智的选择了使用随机数生成这些区间！但是为什么小Hi不直接使用随机数生成购物清单呢？——问那么多做什么！）

提示一：二分法是宇宙至强之法！（真的么？）

提示二：线段树不也是二分法么？

输入

每个测试点（输入文件）有且仅有一组测试数据。

每组测试数据的第1行为一个整数N，意义如前文所述。

每组测试数据的第2行为N个整数，分别描述每种商品的重量，其中第i个整数表示标号为i的商品的重量weight_i。

每组测试数据的第3行为一个整数Q，表示小Hi总共询问的次数。

每组测试数据的第N+4~N+Q+3行，每行分别描述一个询问，其中第N+i+3行为两个整数Li, Ri，表示小Hi询问的一个区间[Li, Ri]。

对于100%的数据，满足N<=10^6，Q<=10^6, 1<=Li<=Ri<=N，0<weight_i<=10^4。

输出

对于每组测试数据，对于每个小Hi的询问，按照在输入中出现的顺序，各输出一行，表示查询的结果：标号在区间[Li, Ri]中的所有商品中重量最轻的商品的重量。

虽说小Ho身经百战，但是面对这样一个他从来没有接触过的问题也一时挠头——数据范围的N和Q都在10^6这个级别，这就意味着至少需要一个时间复杂度为O(Nlog(N))才能解决这个问题。

“如果我对于每个询问都扫描对应的区间，找到最小的值，那么最坏情况和平均情况下都是O(NQ)的时间复杂度，无论数据是不是随即生成的，都不是一个很好的选择。”小Ho如是想道：“这就意味着我需要事先进行一些预处理，使得一些重复计算的东西不再重复计算，才能够将复杂度降低下来。”

小Ho暗自点头，随即想道：“那么如果我对于这整个区间维护一个线段树，维护每棵子树中的最小值，是不是就可以完美解决这个问题了呢？”，于是对小Hi说道：“我觉得线段树可以很好的解决这个问题！。”

“没有错呢~但是相比于线段树来说，我有一个更为容易事先的算法，你想不想听？”小Hi点了点头道。

“更容易实现？就是代码量更小咯？”小Ho作为一个懒病患者自然是很乐意：“快教我！”

“其实你利用线段树也就是这样来减少复杂度的——先预先计算一些区间的最小值，然后把每个询问都拆成若干个计算了最小值的区间，并且统计这些区间的最小值的最小值，从而得出答案的。”小Hi先总结道：“那么其实我可以将统计的区间这样规定——统计所有长度为2的非负整数次幂的区间。”

“整数次幂？也就是长度为1、2、4、8之类的区间咯？”小Ho答道：“我想想……如果我统计的是这些长度的区间的话，用pre_calc[L, Len]表示左边界为L，长度为Len的区间中的最小值——那么对于一个询问[Li, Ri]，我只要找到小于这个区间长度的最大的2的非负整数次幂——T，那么这个区间中的最小值就是min{pre_calc[Li, T], pre_calc[Ri-T+1, T]}咯？”

“反应的很快嘛，你看这样对于每一个询问，是不是就可以只用O(1)的时间复杂度就可以作出回答了，而且就只用一条很简单的语句而已。”小Hi笑道。

“但是我还不知道怎么求这些pre_calc啊？”小Ho不禁问道。

“我先问你，对于所有的i满足1<=i<=N, pre_calc[i, 1]是不是很好求？”小Hi反问道。

“唔……是的，pre_calc[i, 1]就是标号为i的物品的重量weight_i。”小Ho答道。

“那么对于，所有的i, j满足1<=i<=N, 1<2^j<=N，pre_calc[i, 2^j]=min{pre_calc[i, 2^(j-1)], pre_calc[i+2^(j-1), 2^(j-1)]}是不是成立的？”小Hi继续问道。

“额……将一个区间，分解成两个一半长度的区间，利用已经求出的值进行计算，非常的合理！”小Ho赞叹道：“甚至于只需要一个二维循环就能够计算出来了。”

“怎么样！是不是比线段树简单许多呢？”小Hi笑道。

#pragma comment(linker, "/STACK:102400000,102400000")
#include <iostream>
#include <vector>
#include <string>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <ctime>
#include <cstdlib>
#include <set>
#include <map>
#include <cctype>
#include <list>
#include <cmath>
#include <bitset>
#include <queue>
#include <stack>
#include <sstream>
#include <functional>
#include <cassert>
using namespace std;

int RL(char _str[], int _maxlen = 1e8) {
	if (fgets(_str, _maxlen, stdin) == NULL) return -1;
	int i = 0;
	while (_str[i]) ++i;
	if (_str[i - 1] == '\n') _str[--i] = 0;
	return i;
}

inline int RI() {
	int temp;
	scanf("%d", &temp);
	return temp;
}

inline long long RLL() {
	long long temp;
	scanf("%I64d", &temp);
	return temp;
}

inline double RD() {
	double temp;
	scanf("%lf", &temp);
	return temp;
}

template <class T1, class T2>
inline void Max(T1 &a, T2 b) { if (b > a) a = b; }

template <class T1, class T2>
inline void Min(T1 &a, T2 b) { if (b < a) a = b; }

typedef pair<int, int> pii;
typedef pair<double, double> pdd;
typedef long long LL;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const double eps = 1e-8;
const int mod = 100007;
const int maxn = 1e6 + 9;

int N, M, a[maxn];
int dp[20][maxn];

void InitRMQ() {
	for (int i = 0; i < N; ++i) dp[0][i] = a[i];
	for (int i = 1; 1 << i <= N; ++i) {
		for (int j = 0; j + (1 << i) <= N; ++j) {
			dp[i][j] = min(dp[i - 1][j], dp[i - 1][j + (1 << (i - 1))]);
		}
	}
}

int RMQ(int l, int r) {
	int k = log(r - l + 1) / log(2);
	return min(dp[k][l], dp[k][r - (1 << k) + 1]);
}

int main() {
	//std::ios::sync_with_stdio(0);
	//std::cin.tie(0);
#ifdef NIGHT_13
	freopen("in.txt", "r", stdin);
	//freopen("myout.txt", "w", stdout);
#endif
	int CAS = 0;
	scanf("%d", &N);
	for (int i = 0; i < N; ++i) {
		scanf("%d", &a[i]);
	}
	InitRMQ();
	scanf("%d", &M);
	while (M--) {
		int a, b;
		scanf("%d%d", &a, &b);
		printf("%d\n", RMQ(--a, --b));
	}
	return 0;
}

fushenliao

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
hihoCoder #1068 : RMQ-ST算法

#1068 : RMQ-ST算法时间限制:10000ms单点时限:1000ms内存限制:256MB描述小Hi和小Ho在美国旅行了相当长的一段时间之后，终于准备要回国啦！而在回国之前，他们准备去超市采购一些当地特产——比如汉堡（大雾）之类的回国。但等到了超市之后，小Hi和小Ho发现者超市拥有的商品种类实在太多了——他们实在看不过来了！于
复制链接

扫一扫

专栏目录