RMQ求区间最大频率

暂时只做了递增序列的求解方法，就是将原数列处理成个数的数组，类似于前缀和的东西，对这个频率数组求区间最大值，

但是仅仅维护区间max还是不够的，，好比这种情况:

原数列：     1 1 1 2 3 5 5 7 7 

频率数组：  1 2 3 1 1 1 2 1 2 

下标：         1 2 3 4 5 6 7 8 9

如果查询区间3~5，频率数组分别为3 1 1这样，max就是3明显不对，这时加一种处理，去掉开始的重复部分，就是保证查询区间的第一个数的频率数组权值为1。详情见代码：

AC代码：

#include <map>
#include <set>
#include <stack>
#include <cmath>
#include <queue>
#include <bitset>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <fstream>
#include <cstdlib>
#include <sstream>
#include <cstring>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;
//#define   maxd          1010
#define lson l , m , rt << 1
#define rson m + 1 , r , rt << 1 | 1
#define   mc(x, y)     memcpy(x, y, sizeof(x))
#define   ms(x,y)      memset(x,y,sizeof(x))
#define   rep(i,n)      for(int i=0;i<(n);i++)
#define   repf(i,a,b)   for(int i=(a);i<=(b);i++)
#define   PI           pair<int,int>
//#define   mapp            make_pair
#define   FI            first
#define   SE            second
#define   IT            iterator
#define   PB            push_back
#define   Times         10
typedef   long long     ll;
typedef   unsigned long long ull;
typedef   long double   ld;
typedef   pair<int,int > pce;
//#define N 100
const double eps = 1e-10;
const double  pi = acos(-1.0);
const  ll    mod = 1e9+7;
const  int   inf = 0x3f3f3f3f;
//const  ll    INF = (ll)1e18+300;
const int   maxd = 100500 + 10;

int n;
int Q;
int ac[maxd];
int max_[maxd][20];
int min_[maxd][20];
int pre[maxd];

void RMQ() {
    for (int i = 1; i <= n; i++) {
        max_[i][0] = min_[i][0] = pre[i];
    }
    for (int j = 1; j <= 20; j++) {
        /*外层循环要从j开始，因为ST表我们是通过j-1来更新j的，如果是i为外层循环，第一次我们计算出了：
            1, 2, 4, 8,...这些位置，而下一次更新的时候，我们要的是max([1,4], [5, 8])但是我们没有计算出这两个区间的最大值，因为一次层的时候
            没有计算3, 5, 6, 7,...的值，因此必须j为第一层循环。

        */
        for (int i = 1; i <= n; i++) {
            if(i + (1<<j) - 1 <= n) {
                max_[i][j] = max(max_[i][j-1], max_[i + (1<<(j-1))][j - 1]);
                //min_[i][j] = min(min_[i][j-1], min_[i + (1<<(j-1))][j - 1]);
            }
        }
    }
}

int main() {
    while(~scanf("%d", &n) && n) {
        scanf("%d", &Q);
        for (int i = 1; i <= n; i++) {
            scanf("%d", &ac[i]);
        }
        //sort(ac+1, ac+1+n);
        for (int i = 1; i <= n; i ++) {
            if(i == 1) {
                pre[i] = 1;
                continue;
            }
            if(ac[i] == ac[i - 1]) {
                pre[i] = pre[i - 1] + 1;
            }
            else {
                pre[i] = 1;
            }
        }
        RMQ();
        for (int i = 1; i <= Q; i++) {
            int a, b;
            scanf("%d%d", &a, &b);
            int temp = a;
            while(temp <= b && ac[temp] == ac[temp - 1]) {
                temp ++;
            }
            int res = temp - a;
            a = temp;
            if(a <= b) {
                int k = (int)(log(b - a + 1)/log(2.0));
                /*查询操作就像图片中所示，取超过（r - l）/2的长度，分区间比较，而我们就可以取log2并向上取整.
                */
                int max_ans = max(max_[a][k], max_[b - (1<<k) + 1][k]);
                //cout << max_ans << " " << res << endl;
                //int min_ans = min(min_[a][k], min_[b - (1<<k) + 1][k]);
                cout << max(res, max_ans) << endl;
            }
            else {
                cout << res << endl;
            }
        }
    }
}