UVA 11235 Frequent values 线段树/RMQ

　　vjudge 上题目链接：UVA 11235

*******************************************************大白书上解释************************************************************

　　题目大意：给出一个非降序排列的整数数组 a₁,a₂,a₃,...,a_n，你的任务是对于一系列询问 (i, j)，回答 a_i,a_i+1,...,a_j 中出现次数最多的值所出现的次数。

　　输入格式：包含多组数据。每组数据第一行为两个整数 n 和 q（1 <= n, q <= 100000）。第二行包含 n 个非降序排列的整数 a₁,a₂,a₃,...,a_n（-100000 <= a_i <= 100000）。以下 q 行每行包含两个整数 i 和 j（1 <= i <= j <= n），输入结束标志为 n = 0。

　　输出格式：对于每个查询，输出查询结果。

　　分析：应注意到整个数组是非降序的，所有相等元素都会聚集到一起。这样就可以把整个数组进行游程编码（Run Length Encoding, RLE）。比如 -1,1,1,2,2,2,4 就可以编码成 (-1, 1), (1, 2), (2, 3), (4, 1)，其中 (a, b) 表示有 b 个连续的 a。用 value[i] 和 count[i] 分别表示第 i 段的数值和出现次数，num[p], left[p], right[p] 分别表示位置 p 所在段的编号和左右端点位置，则在下图的情况，每次查询（L，R）的结果为以下 3 个部分的最大值：从 L 到 L 所在段的结束处的元素个数（即 right[L] - L + 1）、从 R 所在段的开始处到 R 处的元素个数（即 R - left[R] + 1）、中间第 num[L] + 1 段到第 num[R] - 1 段的 count 的最大值，如图 3-8 所示。

　　

*******************************************************大白书上解释结束************************************************************

　　我的理解：

　　预处理过程主要就 3 个数组：seq[] 就是上述提到的 count[] 数组，记录 seq[i] 第 i 段连续整数的出现次数；pos[i] 表示原数组的第 i 个元素在 seq[] 中处于第几段；preSum[] 则是 seq 数组的前缀和，用于快速求出第 L 段和第 R 段的元素个数。这 3 个数组准备好后，接下来就是求区间最值的问题而已，线段树或者 RMQ 都可以，二者复杂度一样，时间差异可以忽略不计，只不过我更熟悉线段树，感觉 RMQ 的边界有点不容易处理而已。

　　首先是线段树的代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <iostream>
using namespace std;
#define  For(i,s,t)  for(int i = (s); i < (t); ++i)
#define  root       int rt, int l, int r
#define  ll(rt)     ((rt) << 1)
#define  rr(rt)     (ll(rt) | 1)
#define  makemid    int mid = (l + r >> 1)
#define  lson       ll(rt), l, mid
#define  rson       rr(rt), mid + 1, r
const int N = 100005;

int c[N];
vector<int> seq, preSum;
int pos[N] = {0,};

int Max[N << 2];

inline void pushup(int rt) {   Max[rt] = max(Max[ll(rt)], Max[rr(rt)]);   }

void build(root)
{
    if(l == r) {
        Max[rt] = seq[l - 1];
        return;
    }
    makemid;
    build(lson);
    build(rson);
    pushup(rt);
}

int ql, qr;
int query(root)
{
    if(ql <= l && r <= qr) {
        return Max[rt];
    }
    makemid;
    int ret = 0;
    if(ql <= mid) {
        ret = max(ret, query(lson));
    }
    if(qr > mid) {
        ret = max(ret, query(rson));
    }
    return ret;
}

int main()
{
    int n,q;
    while(~scanf("%d",&n), n) {
        scanf("%d", &q);
        seq.clear();
        scanf("%d", c);
        int curValue = c[0], curNum = 1;
        For(i, 1, n) {
            scanf("%d", c + i);
            if(c[i] == curValue) {
                ++curNum;
                pos[i] = pos[i - 1];
            } else {
                seq.push_back(curNum);
                curValue = c[i];
                curNum = 1;
                pos[i] = pos[i - 1] + 1;
            }
        }
        seq.push_back(curNum);
        preSum.clear();
        preSum.push_back(seq[0]);
        int len = seq.size();
        For(i, 1, len) {
            preSum.push_back(preSum[i - 1] + seq[i]);
        }
        build(1, 1, len);
        int x,y;
        while(q--) {
            scanf("%d %d",&x,&y);
            --x;  --y;
            if(pos[x] == pos[y]) {
                printf("%d
", y - x + 1);
                continue;
            }
            int lmax = preSum[pos[x]] - x;
            int rmax = y + 1 - preSum[pos[y] - 1];
            int res = max(lmax, rmax);
            if(pos[y] == pos[x] + 1) {
                printf("%d
", res);
            } else {
                ql = pos[x] + 1  + 1;
                qr = pos[y] - 1  + 1;
                printf("%d
", max(res, query(1, 1, len)));
            }
        }
    }
    return 0;
}

　　然后是 RMQ 的：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <iostream>
using namespace std;
#define  For(i,s,t)  for(int i = (s); i < (t); ++i)
const int N = 100005;

int c[N];
vector<int> seq, preSum;
int pos[N] = {0,};

int d[N][18];
inline void init(int n)
{
    For(i, 0, n) {
        d[i][0] = seq[i];
    }
    for(int j = 1; (1 << j) < n; ++j) {
        for(int i = 0; i + (1 << j) - 1 < n; ++i) {
            d[i][j] = max(d[i][j - 1], d[i + (1 << (j - 1))][j - 1]);
        }
    }
}

inline int rmq(int L, int R)
{
    int k = 0, len = R - L + 1;
    while((1 << (k + 1)) < len)   ++k;
    return max(d[L][k], d[R - (1 << k) + 1][k]);
}

int main()
{
    int n,q;
    while(~scanf("%d",&n), n) {
        scanf("%d", &q);
        seq.clear();
        scanf("%d", c);
        int curValue = c[0], curNum = 1;
        For(i, 1, n) {
            scanf("%d", c + i);
            if(c[i] == curValue) {
                ++curNum;
                pos[i] = pos[i - 1];
            } else {
                seq.push_back(curNum);
                curValue = c[i];
                curNum = 1;
                pos[i] = pos[i - 1] + 1;
            }
        }
        seq.push_back(curNum);
        preSum.clear();
        preSum.push_back(seq[0]);
        int len = seq.size();
        For(i, 1, len) {
            preSum.push_back(preSum[i - 1] + seq[i]);
        }
        init(len);
        int x,y;
        while(q--) {
            scanf("%d %d",&x,&y);
            --x;  --y;
            if(pos[x] == pos[y]) {
                printf("%d
", y - x + 1);
                continue;
            }
            int lmax = preSum[pos[x]] - x;
            int rmax = y + 1 - preSum[pos[y] - 1];
            int res = max(lmax, rmax);
            if(pos[y] == pos[x] + 1) {
                printf("%d
", res);
            } else {
                int ql = pos[x] + 1;
                int qr = pos[y] - 1;
                printf("%d
", max(res, rmq(ql, qr)));
            }
        }
    }
    return 0;
}