POJ3368 离散化+ST算法求解RMQ

给定一个数组，其中的元素满足非递减顺序，要求对于一对起点和终点，回答出其中某个元素重复出现的最大次数。
比如对于-1 -1 1 1 1 1 3 10 10 10，若起点为1，终点为5，则重复出现最大的数是1，其次数为3。
由于原数组是非递减的，因此相同的数字必然是连续出现的。可以对原数组进行压缩，把相同数字构成的区间抽象为一个数，这个数就是这个区间的大小。比如上述数组可以压缩为2 4 1 3。对于原数组中的每个位置，记录它属于第几个区间，并记录每个区间的起点和终点。
对于一对起点和终点，首先计算出它们各属于哪个区间，分为三种情况：
1、属于同一个区间：答案就是它们之间的大小
2、属于相邻区间：找到它们的分隔点，取两个区间大小的较大者
3、属于两个不相邻的区间：取出头尾两个区间，计算它们的大小。对于中间的一个或多个区间，它们都是完整的，因此可以使用求解RMQ问题ST算法，计算出它们中的最大值。最后的答案即为该最大值与头尾两个区间的大小，一共三个数的最大值。

由于要使用ST算法，因此在计算出每个区间大小之后要进行初始化，这样就可以在O(1)的时间内回答出上述第三种情况的中间区间的长度的最大值。

#include <iostream>#include <cmath>#include <cstdio>using namespace std;const int N = 100005;struct Seg{int beg;int end;};int num[N];int pos[N];Seg seg[N];int dpmax[N][30];int segcount;int n, q;const double log_2 = log(2.0);inline int Max(int x, int y) {return x > y ? x : y;}void rmqinit(){for (int i = 0; i <= segcount; ++i) dpmax[i][0] = seg[i].end - seg[i].beg + 1;for (int j = 1; (1 << j) <= segcount; ++j){int t = 1 << j;for (int i = 0; i + t - 1 <= segcount; ++i){dpmax[i][j] = Max(dpmax[i][j - 1], dpmax[i + (t >> 1)][j - 1]);}}}int rmqgetmax(int beg, int end){int t = (int)(log(double(end - beg + 1)) / log_2);return Max(dpmax[beg][t], dpmax[end - (1 << t) + 1][t]);}int main(){int beg, end, ans;num[0] = 1000000;while (scanf("%d", &n) != EOF && n){scanf("%d", &q);segcount = 0;for (int i = 1; i <= n; ++i){scanf("%d", num + i);if (num[i] != num[i - 1]){segcount++;pos[i] = segcount;seg[segcount].beg = seg[segcount].end = i;}else{pos[i] = segcount;seg[segcount].end++;}}rmqinit();while (q--){scanf("%d%d", &beg, &end);if (pos[end] == pos[beg]) ans = end - beg + 1;else if (pos[end] - pos[beg] == 1) ans = Max(seg[pos[beg]].end - beg + 1, end - seg[pos[end]].beg + 1);else{ans = Max(seg[pos[beg]].end - beg + 1, end - seg[pos[end]].beg + 1);ans = Max(ans, rmqgetmax(pos[beg] + 1, pos[end] - 1));}printf("%d\n", ans);}}return 0;}