ST算法模板题 POJ 3264

POJ 3264

题意：给定一串数组，求某个区间范围内的最大值和最小值之差。 

ST算法：能在 O(n*log n) 时间的预处理后，以 O(1) 的时间复杂度在线回答”数列A中下标在 l 到 r 之间的数的最大值或最小值是多少。

#include<iostream>
#include<iomanip>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
#include<queue>
#include<stack>
#include<vector>
#include<map>
#include<set>
#include<sstream>
#include<cstdlib>
#include<time.h>
#define rep(i,a,n) for(int i=a;i<=n;i++)
#define repp(i,n,a) for(int i=n;i>=a;i--)
#define ll long long
#define ull unsigned long long
#define mem(a,x) memset((a),(x),sizeof ((a)))//x只能是0或-1或false或true
#define debug(x) cout<<"X: "<<(x)<<endl
#define de cout<<"************"<<endl
#define lowbit(x) ((x)&(-x))
#define lson rt<<1
#define rson rt<<1|1
#define gcd(a,b) __gcd(a,b)
#define lcm(a,b) a*b/(__gcd(a,b))
#define inf 0x3f3f3f3f//1e9+6e7
#define Eps 1e-8
#define Mod 1e9+7
#define maxn 100010
const double pi=acos(-1.0);
using namespace std;
const int SIZE=200010;
int a[50010];
int f[50010][110];
int ff[50010][110];
int n;
void ST_prework()
{
    for(int i=1;i<=n;i++)
        f[i][0]=a[i];
    int t=log(n)/log(2)+1;
    for(int j=1;j<t;j++)
        for(int i=1;i<=n-(1<<j)+1;i++)
            f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
}
void ST_preworkk()
{
    for(int i=1;i<=n;i++)
        ff[i][0]=a[i];
    int t=log(n)/log(2)+1;
    for(int j=1;j<t;j++)
        for(int i=1;i<=n-(1<<j)+1;i++)
            ff[i][j]=min(ff[i][j-1],ff[i+(1<<(j-1))][j-1]);
}
int ST_query(int l,int r)   //返回a[l]到a[r]中的最大值
{
    int k=log(r-l+1)/log(2);
    return max(f[l][k],f[r-(1<<k)+1][k]);
}
int ST_queryy(int l,int r)  //返回a[l]到a[r]中的最小值
{
    int k=log(r-l+1)/log(2);
    return min(ff[l][k],ff[r-(1<<k)+1][k]);
}
int main()
{

    int q;
    cin>>n>>q;
    rep(i,1,n)
        scanf("%d",&a[i]);
    ST_prework();
    ST_preworkk();
    int l,r;
    while(q--)
    {
        scanf("%d%d",&l,&r);
        cout<<ST_query(l,r)-ST_queryy(l,r)<<endl;
    }
    return 0;
}

或者：

#include<iostream>
#include<iomanip>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
#include<queue>
#include<stack>
#include<vector>
#include<map>
#include<set>
#include<sstream>
#include<cstdlib>
#include<time.h>
#define rep(i,a,n) for(int i=a;i<=n;i++)
#define repp(i,n,a) for(int i=n;i>=a;i--)
#define ll long long
#define ull unsigned long long
#define mem(a,x) memset((a),(x),sizeof ((a)))//x只能是0或-1或false或true
#define debug(x) cout<<"X: "<<(x)<<endl
#define de cout<<"************"<<endl
#define lowbit(x) ((x)&(-x))
#define lson rt<<1
#define rson rt<<1|1
#define gcd(a,b) __gcd(a,b)
#define lcm(a,b) a*b/(__gcd(a,b))
#define inf 0x3f3f3f3f//1e9+6e7
#define Eps 1e-8
#define Mod 1e9+7
#define maxn 100010
const double pi=acos(-1.0);
using namespace std;
const int SIZE=200010;
int a[50010];
int f[50010][110];
int ff[50010][110];
int n;
void ST_prework()
{
    for(int i=0;i<=n-1;i++)
        f[i][0]=a[i];
    int t=log(n)/log(2)+1;
    for(int j=1;j<t;j++)
        for(int i=0;i<=n-(1<<j);i++)
            f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
}
void ST_preworkk()
{
    for(int i=0;i<=n-1;i++)
        ff[i][0]=a[i];
    int t=log(n)/log(2)+1;
    for(int j=1;j<t;j++)
        for(int i=0;i<=n-(1<<j);i++)
            ff[i][j]=min(ff[i][j-1],ff[i+(1<<(j-1))][j-1]);
}
int ST_query(int l,int r)   //返回a[l]到a[r]中的最大值
{
    int k=log(r-l+1)/log(2);
    return max(f[l][k],f[r-(1<<k)+1][k]);
}
int ST_queryy(int l,int r)  //返回a[l]到a[r]中的最小值
{
    int k=log(r-l+1)/log(2);
    return min(ff[l][k],ff[r-(1<<k)+1][k]);
}
int main()
{

    int q;
    cin>>n>>q;
    //rep(i,1,n)
    for(int i=0;i<n;i++)
        scanf("%d",&a[i]);
    ST_prework();
    ST_preworkk();
    int l,r;
    while(q--)
    {
        scanf("%d%d",&l,&r);
        l--;r--;
        cout<<ST_query(l,r)-ST_queryy(l,r)<<endl;
    }
    return 0;
}