poj 3368 && hdu 1806 Frequent values（数据结构：RMQ）

给出一段非递减序列，让你统计任意（L，R）区间中出现最多的数出现的次数

直接暴力毫无疑问会超时，当然这道题也可以用线段树做，但是会比较慢

正好看到了白书的RMQ，就对照着写了下， RMQ（范围最小值）对一段区间进行预处理，复杂度为O（nlogn）

而查询的复杂度为O（1）

这道题的正解是先对序列进行游程编码（Run Length Encoding）

划为形似（a,b）的二元组表示有b个连续的a

再记录每个位置所在的段号num，以及左侧left和右侧边界right

对于查询（L， R）其结果即为max(right[L]-L+1, R-left[R]+1, RMQ(count, num[L]+1, num[R]-1))

516ms如下：

#include <cstdio>
#include <cstring>
#include <iostream>
#define MAXN 101000
using namespace std;

struct NODE {
    int val, time;
}node[MAXN+10];//第i段的值以及连续出现的次数

int a[MAXN+10];//输入
int d[MAXN+10][25];//RMQ存储

struct Point {
    int pos, left, right;
}point[MAXN+10];

void RMQ_init(NODE b[], int cnt) {
    for(int i=1; i<=cnt; ++i) {
        d[i][0] = b[i].time;
        //printf("a[%d].time = %d\n", i, a[i].time);
    }
    for(int j=1; (1<<j)<=cnt; ++j) {
        for(int i=1; i+(1<<j)-1<=cnt; ++i) {
            d[i][j] = max(d[i][j-1], d[i+(1<<(j-1))][j-1]);
            //printf("d[%d][%d] = %d\n", i, j, d[i][j]);
        }
    }
}

int RMQ(int L, int R) {
    int k = 0;
    while((1<<(k+1)) <= R-L+1)
        ++k;
    return max(d[L][k], d[R-(1<<k)+1][k]);
    //这里写了一个bug，写成了max(d[L][k], d[L+(1<<k)-1][k]),可以发现这里的右侧边界很有可能超出了R，这是不对的
}

int main(void) {
    int n, q, x, cnt;
    while(scanf("%d", &n) && n) {
        scanf("%d", &q);

        memset(d, 0, sizeof(d));
        memset(node, 0, sizeof(node));
        memset(point, 0, sizeof(point));

        cnt = 0;
        for(int i=1; i<=n; ++i) {
            scanf("%d", &a[i]);

            if(a[i] != node[cnt].val)
                ++cnt;
            node[cnt].val = a[i];
            ++node[cnt].time;

            //表示该位置所在段的编号和左右端点位置
            if(i>1 && a[i] == a[i-1]) {
                point[i].pos = point[i-1].pos;
                point[i].left = point[i-1].left;
                point[i].right = i;
            } else {
                point[i].pos = cnt;
                point[i].left = point[i].right = i;
            }
        }

        for(int i=n-1; i>0; --i) {
            if(a[i] == a[i+1]) {
                point[i].right = point[i+1].right;
            }
            //printf("point[%d]: pos=%d\tleft=%d\tright=%d\n", i, point[i].pos, point[i].left, point[i].right);
        }
/*
        for(int i=1; i<=cnt; ++i) {
            printf("node[%d] : val = %d\ttime = %d\n", i, node[i].val, node[i].time);
        }
*/

        RMQ_init(node, cnt);

        int L, R, ans;
        for(int i=1; i<=q; ++i) {
            scanf("%d%d", &L, &R);
            //printf("point[%d]: pos=%d\tleft=%d\tright=%d\n", L, point[L].pos, point[L].left, point[L].right);
            //printf("point[%d]: pos=%d\tleft=%d\tright=%d\n", R, point[R].pos, point[R].left, point[R].right);
            if(a[L] == a[R]) 
                ans = R-L+1;
            else if(point[L].pos+1 == point[R].pos) {
                ans = max(R-point[R].left+1, point[L].right-L+1);
            } 
            else {
                ans = RMQ(point[L].pos+1, point[R].pos-1);
                //printf("RMQ = %d\n", ans);
                ans = max(ans, R-point[R].left+1);
                ans = max(ans, point[L].right-L+1);
            }
            printf("%d\n", ans);
        }

    }
}