POJ-3264 RMQ 线段树与ST

题意：询问区间最大值与最小值的差。

此题解法很多，线段树是最容易理解的了，除了暴力！

#include<iostream>
#include<cstdio>
#include<cstring>
#define INF 10000000
using namespace std ;
const int N=200005;

struct node{
	int l,r;
	int Mi,Ma;
}tree[N*3];

int Max,Min,num;

void build(int l,int r,int i)
{
	tree[i].l=l;
	tree[i].r=r;
	if(l==r){
		scanf("%d",&num);
		tree[i].Mi=tree[i].Ma=num;
		return;
	}
	int mid=(l+r)>>1;
	build(l,mid,i<<1);
	build(mid+1,r,i<<1|1);
	tree[i].Mi=min(tree[i<<1].Mi,tree[i<<1|1].Mi);
	tree[i].Ma=max(tree[i<<1].Ma,tree[i<<1|1].Ma);
}

void query(int l,int r,int i)
{
	if(tree[i].Ma<=Max&&tree[i].Mi>=Min) return;
	if(tree[i].l==l&&tree[i].r==r){
		Min=min(tree[i].Mi,Min);
		Max=max(tree[i].Ma,Max);
		return;
	}
	int mid=(tree[i].l+tree[i].r)>>1;
	if(r<=mid) query(l,r,i<<1);
	else if(l>mid) query(l,r,i<<1|1);
	else{
		query(l,mid,i<<1);
		query(mid+1,r,i<<1|1);
	}
}

int main()
{
	int n,m,l,r;
    while(scanf("%d %d",&n,&m)!=EOF) {
		build(1,n,1);
		while(m--){
			Max=-INF,Min=INF;
			scanf("%d %d",&l,&r);
			query(l,r,1);
			printf("%d\n",Max-Min);
		}
    }
 	return 0 ;
}

ST算法：

也该是基于DP，倍增和分治思想。

设有：dp(i,j) 它表示[i,i+2^j-1]区间的最小值或最大，这样，i，j所表示的区间不是简单的i到j的区间。例如：dp[1][3]它所包含区间的最值在[1,8],显然减少了内存。这样的设计还有一种好处，也是该算法的关键：具有递归关系，即：

                                 dp[i][j]=max(dp[i][j-1],dp[i+(1<<(j-1))][j-1])

也就是说：把区间划分了两个区间。这样一次循环下去就可以把[1,n]的区间全部求完了。

虽然表示完所有的空间，但所求的区间(a,b)不能直接代，而a和b也不一定满足:a+2^j-1=b,在逻辑上说不通，那么就必须把区间划分成两个。先寻找最大k。k=(int)log2(b-a+1)，那么:

[a,b]=[a,a+2^k-1]+[b-2^k+1,b];这里是个数学计算可以推出这样的结果。比如:[1 13]=[1,8]+[6,13],显然有重叠的部分，但不会影响结果。

根据dp的定义:Max[a->b]=max(dp[a][k],[b-(1<<k)+1][k]);

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std ;

const int N=50005;
int Max[N][20],Min[N][20];
int a[N];

void pre(int n)
{
	for(int i=1;i<=n;i++){
		Max[i][0]=Min[i][0]=a[i];
	}
	int k=trunc(log2(n+1));
	for(int j=1;j<=k;j++)
		for(int i=1;i+(1<<j)-1<=n;i++){
			Max[i][j]=max(Max[i][j-1],Max[i+(1<<(j-1))][j-1]);
		    Min[i][j]=min(Min[i][j-1],Min[i+(1<<(j-1))][j-1]);
		}		
}

void Query(int l,int r)
{
	int k=trunc(log2(r-l+1));
	int ans1=max(Max[l][k],Max[r-(1<<k)+1][k]);
	int ans2=min(Min[l][k],Min[r-(1<<k)+1][k]);
	printf("%d\n",ans1-ans2);
}

int main()
{
	int n,m,l,r;
    while(~scanf("%d %d",&n,&m)){
		for(int i=1;i<=n;i++){
			scanf("%d",a+i);
		}
		pre(n);
		while(m--){
			scanf("%d %d",&l,&r);
			Query(l,r);
		}
    }
 	return 0 ;
}