POJ 3368 最频繁值 （裸ST）

#include<iostream>
#include<cstdlib>
#include<sstream>
#include<cstdio>
#include<stack>
#include<cstdio>
#include<map>
#include<set>
#include<queue>
#include<cstring>
#include<cmath>
#include<vector>
#include<math.h>
#include<algorithm>
#define REP(i,a,b) for(int i=a;i<b;i++)
#define REPP(i,a,b) for(int i=a;i<=b;i++)
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=100005;
int lb[maxn];
int F[maxn][100];
void Initlog()//求解所有log值，将其存储到数组lb[]中
{
    lb[0]=-1;
    REPP(i,1,maxn)
    {
        lb[i]=(i&(i-1))?lb[i-1]:lb[i-1]+1;
    }
}
int a[maxn];
void _Cin(int n,int q)
{
    memset(a,0,sizeof(a));
    memset(F,0,sizeof(F));
    REPP(i,1,n)
    {
        scanf("%d",&a[i]);
        if(i==1) F[i][0]=1;
        else
        {
            if(a[i]==a[i-1]) F[i][0]=F[i-1][0]+1;
            else F[i][0]++;
        }
    }
}
void ST_create(int n)//每个测试用例的元素个数都不同，所以用n作为参数
{
    REPP(j,1,lb[n])
    {
        REPP(i,1,(n-(1<<j)+1))
        {
            F[i][j]=max(F[i][j-1],F[i+(1<<(j-1))][j-1]);
        }
    }
}
int RMQ(int l,int r)
{
    if(l>r) return 0;
    int k=lb[r-l+1];
    return max(F[l][k],F[r-(1<<k)+1][k]);
}
int l,r,n,q;
int main()
{
    while(scanf("%d",&n)&&n)
    {
        Initlog();
        scanf("%d",&q);
        _Cin(n,q);
        ST_create(n);
        REPP(j,1,q)
        {
            scanf("%d%d",&l,&r);
            int t=l;
            while(t<=r && a[t]==a[t-1])
                t++;
            printf("%d\n",max(t-l,RMQ(t,r)));
        }
    }
    return 0;
}