POJ 2452 Sticks Problem（ST+二分查找）

Description
给出一长度为n的序列，序列没有重复的数值，问j-i的最大值，其中j,i满足a[i]< a[k]< a[j],i< k< j，如果不存在这样的i,j则输出-1
Input
多组用例，每组用例第一行为一整数n表示序列长度，第二行n个不同的整数表示该序列，以文件尾结束输入
Output
对于每组用例，输出满足条件的j-i的最大值，如果不存在满足条件的i和j则输出-1
Sample Input
4
5 4 3 6
4
6 5 4 3
Sample Output
1
-1
Solution
对于每一个i，从区间[i,n]中找到使得i为区间最小值位置的最大区间[i,t]，然后再在区间[i,t]中找到最大值的位置j，用j-i更新答案即可，区间最值可以用ST维护，查找需要二分，最后答案如果为0说明没有满足条件的i,j
Code

#include<cstdio>#include<iostream>#include<algorithm>#include<cstring>#include<cmath>using namespace std;#define maxn 55555#define INF 0x3f3f3f3fint n,a[maxn],sta[maxn],p,smax[maxn][22],smin[maxn][22];void ST(){    for(int i=0;i<n;i++)smax[i][0]=smin[i][0]=i;    int k=(int)(log(1.0*n)/log(2.0));    for(int j=1;j<=k;j++)           for(int i=0;i+(1<<j)<=n;i++)           {            if(a[smax[i][j-1]]>a[smax[i+(1<<(j-1))][j-1]])smax[i][j]=smax[i][j-1];            else smax[i][j]=smax[i+(1<<(j-1))][j-1];            if(a[smin[i][j-1]]<a[smin[i+(1<<(j-1))][j-1]])smin[i][j]=smin[i][j-1];            else smin[i][j]=smin[i+(1<<(j-1))][j-1];         }}int query_max(int l,int r)  {      int k=(int)(log(1.0*(r-l+1))/log(2.0));      if(a[smax[l][k]]<a[smax[r-(1<<k)+1][k]]) return smax[r-(1<<k)+1][k];    return smax[l][k];}int query_min(int l,int r)  {      int k=(int)(log(1.0*(r-l+1))/log(2.0));      if(a[smin[l][k]]>a[smin[r-(1<<k)+1][k]]) return smin[r-(1<<k)+1][k];    return smin[l][k];}int binary_search(int x,int l,int r){    while(l<r)    {        int mid=(l+r+1)/2;        if(query_min(x,mid)==x)l=mid;        else r=mid-1;    }    return l;}int main(){    while(~scanf("%d",&n))    {        for(int i=0;i<n;i++)            scanf("%d",&a[i]);        ST();        int ans=0;        for(int i=0;i<n;i++)        {            int j=query_max(i,binary_search(i,i,n-1));            ans=max(ans,j-i);        }        if(ans)printf("%d\n",ans);        else printf("-1\n");    }    return 0;}