poj 3264 RMQ与线段树

用线段树和RMQ都可以做的题，给出一个数组，多次查询一个区间内之间最大值与最小值的差

RMQ算法，使用一个二维数组预处理，使得查询时实现O（1）的查询效率.

用一个二维数组f(i,j)记录区间[i,i+2^j-1](持续2^j个)区间中的最小值。其中f[i,0] = a[i]; 
所以，对于任意的一组(i,j)，f(i,j) = min{f(i,j-1),f(i+2^(j-1),j-1)}来使用动态规划计算出来。 
比如求区间[m,n]中a的最大值，找到一个数k使得2^k<n-m+1. 
这样，可以把这个区间分成两个部分：[m,m+2^k-1]和[n-2^k+1,n].我们发现，这两个区间是已经初始化好的. 
前面的区间是f(m,k)，后面的区间是f(n-2^k+1,k). 这样求这两个区间内最大值中大的那个可以求出整个区间内的最大值

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define maxn 100005
using namespace std;
int n,q,a[maxn];
int amax[maxn][18],amin[maxn][18];
void init(){
	int m=floor(log((double)n)/log(2.0));
	for(int i=1;i<=n;i++) amax[i][0]=amin[i][0]=a[i];
	for(int i=1;i<=m;i++)
		for(int j=n;j;j--){
			amax[j][i]=amax[j][i-1];
			amin[j][i]=amin[j][i-1];
			if(j+(1<<(i-1))<=n){
				amax[j][i]=max(amax[j][i],amax[j+(1<<(i-1))][i-1]);
				amin[j][i]=min(amin[j][i],amin[j+(1<<(i-1))][i-1]);
			}
		}
}
int query(int l,int r){
	int m=floor(log((double)(r-l+1))/log(2.0));
	int Max=max(amax[l][m],amax[r-(1<<m)+1][m]);
	int Min=min(amin[l][m],amin[r-(1<<m)+1][m]);
	return Max-Min;
}
int main(){
	while(scanf("%d%d",&n,&q)!=EOF){
		for(int i=1;i<=n;i++) scanf("%d",&a[i]);
		init();
		while(q--){
			int l,r;
			scanf("%d%d",&l,&r);
			printf("%d\n",query(l,r));
		}
	}
	return 0;
}

线段树做法

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
using namespace std;
const int maxv=50005;
const int maxn=300000;
int max(int a,int b)
{
	if(a>b)
	return a;
	return b;
}
int min(int a,int b)
{
	if(a<b)
	return a;
	return b;
}
int a[maxv],tl[maxn],tr[maxn],tmax[maxn],tmin[maxn];
void xdshu(int v,int f,int t)
{
	tl[v]=f;
	tr[v]=t;
	if(f==t)
	{
		tmax[v]=a[f];
		tmin[v]=a[f];
		return ;
	}
	int mid=(f+t)>>1;
	xdshu(v<<1,f,mid);
	xdshu((v<<1)|1,mid+1,t);
	tmax[v]=max(tmax[v<<1],tmax[(v<<1)|1]);
	tmin[v]=min(tmin[v<<1],tmin[(v<<1)|1]);
}
int query(int v,int l,int r,int flag)
{
	if(l==tl[v]&&r==tr[v])
	{
		return (flag?tmin[v]:tmax[v]);
	}
	int mid=(tl[v]+tr[v])>>1;
	if(r<=mid)
	{
		return query(v<<1,l,r,flag);
	}
	else if(l>mid)
	{
		return query((v<<1)|1,l,r,flag);
	}
	else
	{
		if(flag) return min(query(v<<1,l,mid,flag),query((v<<1)|1,mid+1,r,flag));
		else return max(query(v<<1,l,mid,flag),query((v<<1)|1,mid+1,r,flag));
	}
}
int main()
{
	int n,q,l,r,i;
	while(scanf("%d%d",&n,&q)!=EOF)
	{
		for(i=1;i<=n;i++)
		{
			scanf("%d",&a[i]);
		}
		xdshu(1,1,n);
		while(q--)
		{
			scanf("%d%d",&l,&r);
			printf("%d\n",query(1,l,r,0)-query(1,l,r,1));
		}
	}
	return 0;
}