题目:http://poj.org/problem?id=3264
这个题目意思就是给你一个包含你个数列,然后问你m(m<=200000)次i到j之间的最大值和最小值的差值。显然典型的Sparse_Table算法。
一年前就看过着个算法,当时只是很浅显的浏览了一遍,一直感觉这个方法很巧妙,一直想看看,结果忘记是什么知识了,终于在看RMQ时看到了,所以就好好看了了,这个方法真的很伟大,想出来这个方法的人更伟大、、、
Sparse_Table算法其实运用的是动态规划的思想(DP),把一个区间分成两个相等的区间,然后这个区间的最大(最小)值就是这两小区间的最大(最小)值的最大(最小)值。dp[i][j]表示的是区间[i,i+(1<<j)-1]这个区间的最值,显然这个区间包含的个数是偶数个,所以可以分解为两个小的区间dp[i][j-1]]和dp[i+(1<<(j-1))][j-1]这两个区间,计算后就会发现这个区间就是原来[i,i+(1<<(j-1))-1]和[i+(1<<(j-1)),i+(1<<j)-1]这个两个区间,这两个区间和起来刚好等于前面的那个区间。因此可以得到状态转移方程dp[i][j]=max(min)(dp[i][j-1],dp[i+(1<<(j-1))][j-1])。而我们在查找区间 i到j的时候,同样可以找到两个区间(可能会有重复),使得这俩个区间并起来等会i到j这个区间。这连个区间也很好找到,只需找到看,使其满足(1<<(k+1))>j-i&&(1<<k)<=j-i;然后这个两个区间就可以得到了,即:[i,i+(1<<k)-1]和[j-(1<<k),j]。也就是dp[i][k]和dp[j-(1<<k)+1][k]。
以下是这个题的代码
#include <stdio.h>
#include <iostream>
#include <string.h>
#include <math.h>
#define N 50005
using namespace std;
struct ss
{
int max,min;
}maxx[N][30];
int main()
{
int n,m,i,j;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)
{
scanf("%d",&maxx[i][0].max);
maxx[i][0].min=maxx[i][0].max;
}
for(j=1;(1<<j)<=n;j++)
{
for(i=1;i<=n&&i+(1<<j)-1<=n;i++)
{
maxx[i][j].max=max(maxx[i][j-1].max,maxx[i+(1<<(j-1))][j-1].max);
maxx[i][j].min=min(maxx[i][j-1].min,maxx[i+(1<<(j-1))][j-1].min);
}
}
int star,end;
for(i=0;i<m;i++)
{
scanf("%d%d",&star,&end);
j=log(end-star+0.0)/log(2.0);
/*if(end==star)
j=0;*/
int ans=max(maxx[star][j].max,maxx[end-(1<<j)+1][j].max)-min(maxx[star][j].min,maxx[end-(1<<j)+1][j].min);
printf("%d\n",ans);
}
}