NYOJ-士兵杀敌（三） 线段树RMQ算法求区间的最大差值

士兵杀敌（三）

时间限制： 2000 ms  |  内存限制： 65535 KB

难度： 5

描述

南将军统率着N个士兵，士兵分别编号为1~N,南将军经常爱拿某一段编号内杀敌数最高的人与杀敌数最低的人进行比较，计算出两个人的杀敌数差值，用这种方法一方面能鼓舞杀敌数高的人，另一方面也算是批评杀敌数低的人，起到了很好的效果。

所以，南将军经常问军师小工第i号士兵到第j号士兵中，杀敌数最高的人与杀敌数最低的人之间军功差值是多少。

现在，请你写一个程序，帮小工回答南将军每次的询问吧。

注意，南将军可能询问很多次。

输入

只有一组测试数据
第一行是两个整数N,Q，其中N表示士兵的总数。Q表示南将军询问的次数。(1<N<=100000,1<Q<=1000000)
随后的一行有N个整数Vi(0<=Vi<100000000)，分别表示每个人的杀敌数。
再之后的Q行，每行有两个正正数m,n，表示南将军询问的是第m号士兵到第n号士兵。

输出

对于每次询问，输出第m号士兵到第n号士兵之间所有士兵杀敌数的最大值与最小值的差。

样例输入

5 2
1 2 6 9 3
1 2
2 4

样例输出

1
7

解题思路：

这道题一定要用RMQ算法写，不然肯定超时，我开始使用线段树求最大和最小，妥妥的超时了，然后只能使用动态规划的RMQ了，关于RMQ算法大家可以去看我的博客，里面有详细介绍点击打开链接http://blog.csdn.net/wang_heng199/article/details/75103476

下面是线段树代码：

#include <stdio.h>
#define Max 100010
struct node{
	int p, q;
}sum[Max<<2];

int max(int a, int b)
{
	return a>b?a:b;
}
int min(int a, int b) 
{
	return a<b?a:b;
}
void Build(int l, int r, int rt)//建树 
{
	int t;
	 if(l == r) {
	 	scanf("%d", &t);
	 	sum[rt].p = sum[rt].q = t;
	 	return ;
	 }
	int m = (l + r) >> 1;
	Build(l, m, rt<<1);//递归建树 
	Build(m+1, r, rt<<1|1);
    sum[rt].p = max(sum[rt<<1].p, sum[rt<<1|1].p);
    sum[rt].q = min(sum[rt<<1].q, sum[rt<<1|1].q);
}
node quest(int L, int R, int l, int r, int rt)//查询 
{
	if(L <= l && r <= R)
	{
		return sum[rt];
	}
	int m = (l + r) >> 1;
	node res;
	res.p = 0;
	res.q = Max;
	if(L <= m) 
	{
	  res.p = max(res.p, quest(L,R,l,m,rt<<1).p);
	   res.q = min(res.q, quest(L,R,l,m,rt<<1).q);
	}
	if(m < R) 
	{
	   res.p = max(res.p, quest(L,R,m+1,r,rt<<1|1).p);
	   res.q = min(res.q, quest(L,R,m+1,r,rt<<1|1).q);
	}
	return res;
}
int main()
{
	int N, M;
	int a, b;
	scanf("%d%d", &N, &M);
	Build(1,N,1);
	while(M--)
	{
		scanf("%d%d", &a, &b);
		node ss = quest(a, b, 1, N, 1);
		printf("%d\n",ss.p-ss.q);
	}
	return 0;
}

下面是动态规划RMQ算法代码：

/*RMQ*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#define N 200000+10
using namespace std;
int dmax[N][30],dmin[N][30];///2^16>N
vector<int>Q;
void RMQ_init(const vector<int>& A)
{
	int n = A.size();
	for (int i = 0; i < n; i++)
		dmax[i][0]=dmin[i][0]= A[i]; //初始化
	for (int j = 1; (1 << j) <= n; j++)
		for (int i = 0; i + (1 << j) - 1 <n; i++) { ///i从0开始保证了区间的无缝覆盖
			dmax[i][j] = max(dmax[i][j - 1], dmax[i + (1 << (j - 1))][j - 1]);
			dmin[i][j] = min(dmin[i][j - 1], dmin[i + (1 << (j - 1))][j - 1]);
		}
}
int RMQ_max(int L, int R)
{
	int k = 0;
	while ((1 << (k + 1)) <= R - L + 1) k++;///如果2^(k+1)<=R-L+1,那么看还可以加1
	return max(dmax[L][k], dmax[R - (1 << k) + 1][k]);
}
int RMQ_min(int L, int R)
{
	int k = 0;
	while ((1 << (k + 1)) <= R - L + 1) k++;
	return min(dmin[L][k], dmin[R - (1 << k) + 1][k]);
}

int main()
{
	int n,m,i,a,L,R;
	scanf("%d%d",&n,&m);
	for(i=0; i<n; i++) {
		scanf("%d",&a);
		Q.push_back(a);
	}
	RMQ_init(Q);
	while(m--) {
		scanf("%d%d",&L,&R);
		printf("%d\n",RMQ_max(L-1,R-1)-RMQ_min(L-1,R-1));///预处理是从0开始的，所以区间边界都减一
	}
}

题目链接： 点击打开链接http://acm.nyist.net/JudgeOnline/problem.php?pid=119