UVA 11235 Frequent Values ---RMQ

大白书上的例题，具体讲解见大白书，最好用用一个Log数组直接求k，这样就是纯O(1)了

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define N 100003

int a[N],num[N],le[N],ri[N],cnt[N];
int d[N][],n,type;

void RMQ_init()
{
    int i,j;
    for(i=;i<=type;i++)
        d[i][] = cnt[i];
    for(j=;(<<j)<=type;j++)
    {
        for(i=;i+(<<j)-<=type;i++)
            d[i][j] = max(d[i][j-],d[i+(<<(j-))][j-]);
    }
}

int LOG[];
void getLog(int n)
{
    for(int i=;i<=n;i++)
        LOG[i] = (int)(log((double)i)/log(2.0));
}

int RMQ(int l,int r)
{
    int k = LOG[R-L+];
    return max(d[l][k],d[r-(<<k)+][k]);
}

int main()
{
    int q,i,j,pos;
    int l,r;
    while(scanf("%d",&n)!=EOF && n)
    {
        scanf("%d",&q);
        for(i=;i<=n;i++)
            scanf("%d",&a[i]);
        a[] = -;
        type = ;
        for(pos=;pos<=n;pos++)
        {
            if(a[pos] != a[pos-])
            {
                if(pos !=  )
                    cnt[type] = pos-le[pos-];
                num[pos] = ++type;
                for(j=le[pos-];j<=pos-;j++)
                    ri[j] = pos-;
                le[pos] = pos;
                if(pos == n)
                    cnt[type] = ,ri[pos] = pos;
            }
            else
            {
                le[pos] = le[pos-];
                num[pos] = num[pos-];
                if(pos == n)
                {
                     cnt[type] = pos - le[pos] + ;
                     for(j=le[pos];j<=n;j++)
                        ri[j] = n;
                }
            }
        }
        RMQ_init();
        while(q--)
        {
            scanf("%d%d",&l,&r);
            if(num[l] == num[r])
            {
                printf("%d\n",r-l+);
                continue;
            }
            int lmax = ri[l]-l+;
            int rmax = r-le[r]+;
            int mmax;
            if(num[l]+ > num[r]-)
                mmax = ;
            else
                mmax = RMQ(num[l]+,num[r]-);
            printf("%d\n",max(mmax,max(lmax,rmax)));
        }
    }
    return ;
}