poj 3368（ST）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define lowbit(x) x&-x
#define LL long long
#define Max(a,b) a>b?a:b
#define mem(a,b) memset(a,b,sizeof(a))
#define maxn 100005
using namespace std;
int f[maxn][20];
int A[maxn];
int B[maxn];
int n,m;
void init()
{
    mem(f,0);
    int k=log2(n);
    for(int i=1;i<=n;i++)//等于n
       f[i][0]=B[i];
    for(int j=1;j<=k;j++)
    {
        for(int i=1;i+(1<<j)-1<=n;i++)//i+（1<<j)-1<=n
        {
            f[i][j]=Max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
        }
    }
}
int query(int l,int r)
{
    int k=log2(r-l+1);
    return Max(f[l][k],f[r-(1<<k)+1][k]);
}
int main()
{

    while(~scanf("%d",&n) && n)
    {
        mem(B,0);//清0
        scanf("%d%d",&m,&A[1]);
        B[1]=1;
        for(int i=2;i<=n;i++)
        {
            scanf("%d",&A[i]);
            if(A[i]==A[i-1])B[i]=B[i-1]+1;
            else B[i]=1;
        }
        init();
        while(m--)
        {
            int l,r;
            int t=1,be=0,res=0;
            scanf("%d%d",&l,&r);
            for(int i=l;i<r;i++)
            {
                if(A[i]==A[i+1])t++;
                else {be=i+1;
                break;}
            }
            res=query(be,r);
            res=Max(res,t);
            if(t==r-l+1)printf("%d\n",t);
            else printf("%d\n",res);
        }
    }
}
/*画斜杠的部分容易出错，此题ST算法与线段树时间基本没有差别，cin会超时，取消同步依旧超时。

Frequent values  POJ - 3368 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define lowbit(x) x&-x
#define LL long long
#define Max(a,b) a>b?a:b
#define mem(a,b) memset(a,b,sizeof(a))
#define maxn 100005
using namespace std;
int tr[maxn*4];
int A[maxn];
int B[maxn];
void build(int root,int l,int r)
{
    if(l==r)tr[root]=A[l];
    else{
        int m=(l+r)/2;
        build(root*2,l,m);
        build(root*2+1,m+1,r);
        tr[root]=Max(tr[root*2],tr[root*2+1]);
    }
}
int query(int rt,int l,int r,int lf,int rg)
{
    int q1,q2;
    if(r<lf||l>rg)return 0;
    if(lf<=l&&r<=rg)return tr[rt];
    int m=(l+r)/2;
    q1=query(rt*2,l,m,lf,rg);
    q2=query(rt*2+1,m+1,r,lf,rg);
    return Max(q1,q2);
}
int main()
{
    int n,m;
    while(~scanf("%d",&n) && n)
    {
        mem(tr,-1);
        mem(B,0);
        mem(A,0);
        scanf("%d",&m);
        scanf("%d",&B[1]);
        A[1]=1;
        for(int i=2;i<=n;i++)
        {
            scanf("%d",&B[i]);
            if(B[i]==B[i-1])A[i]=A[i-1]+1;
            else A[i]++;
        }//预处理
        build(1,1,n);
        int l,r;
        for(int i=0;i<m;i++)
        {
            int be=0,w=0,res=0,t=1;//变量要每次更新
            scanf("%d%d",&l,&r);
            for(int i=l;i<r;i++)
            {
                if(B[i]==B[i+1])t++;
                else {be=i+1;
                //cout<<"be="<<be<<endl;
                break;}
            }
             w=query(1,1,n,be,r);
             res=Max(t,w);
            if(t==r-l+1)printf("%d\n",t);
            else printf("%d\n",res);
        }
    }
}
/*线段树：
由于变量没更新wa了好几次，把要查询的区间分成两个区间去查需，左边的区间对答案有影响，拿出来单独处理，剩下的交给线段树查询。