POJ 3264 Balanced Lineup（ST）

Description
在一组数中，查询某个区间内的最大数与最小数的差
Input
第一行两个整数n和q分别表示数的个数和查询次数，之后n行每行一个整数表示该数列，最后q行每行两个整数a，b表示查询区间
Output
对于每次查询，输出该区间内最值之差
Sample Input
6 3
1
7
3
4
2
5
1 5
4 6
2 2
Sample Output
6
3
0
Solution
ST算法，大概的意思就是动态规划，以求最大值为例子，Max[i][j],代表的是从第i个开始，到第2^j个中最大的，显然Max[i][0]=a[i],即那个数本身。状态转移方程为：
Max[i][j]=max(Max[i][j-1],Max[i+(1<<(j-1))][j-1])，
意思很明显，就是前半段和后半段的最大值，[i ,2^(j-1)], [2^(j-1),j]的最大值.
Code

#include<cstdio>#include<cstring>#include<iostream>#include<cmath>using namespace std;#define maxn 200010int n,q;int l,r;int hmin[maxn][20],hmax[maxn][20];void init(){    int k=(int)(log(n*1.0)/log(2.0));    for(int j=1;j<=k;j++)        for(int i=0;i+(1<<j)<=n;i++)        {            hmin[i][j]=min(hmin[i][j-1],hmin[i+(1<<(j-1))][j-1]);            hmax[i][j]=max(hmax[i][j-1],hmax[i+(1<<(j-1))][j-1]);        }}int main(){    while(scanf("%d%d",&n,&q)!=EOF)    {        for(int i=0;i<n;i++)        {            scanf("%d",&hmin[i][0]);            hmax[i][0]=hmin[i][0];        }        init();        while(q--)        {            scanf("%d%d",&l,&r);            l--;r--;            int k=(int)(log((r-l+1)*1.0)/log(2.0));            int ans1=min(hmin[l][k],hmin[r-(1<<k)+1][k]);            int ans2=max(hmax[l][k],hmax[r-(1<<k)+1][k]);            printf("%d\n",ans2-ans1);        }    }    return 0;}