【BZOJ 2724】 蒲公英

【题目链接】

https://www.lydsy.com/JudgeOnline/problem.php?id=2724

【算法】

分块算法在线维护区间众数

分块算法的精髓就在于 ： 大段维护，局部朴素

这一题，我们可以将序列分成T段，那么每一段的长度就是(N/T)

对于每个询问，设l处于第p段，r处于第q段，那么 ：

若p = q,用朴素算法计算出区间众数即可

否则，将这个序列分为三段 ：

1.[L,R[p]] 2. [L[p+1],R[q-1]] 3. [L[q],r]

显然，区间众数只可能是 [L[p+1],R[q-1]]中的众数或[L,R[p]]中的一个数，或[L[q],r]中的一个数

不妨预处理所有以“段边界”为端点每个数出现的次数和区间众数

　　　 当T取sqrt3(N)(开三次方）时，时间复杂度是非常优秀的 : O(n^(5/3))

【代码】

#include<bits/stdc++.h>
using namespace std;
#define MAXN 40010
const int INF = 2e9;

int i,j,k,n,m,l,r,l0,r0,lastans,t,block,len;
int a[MAXN],val[MAXN],L[],R[],pos[MAXN],rk[MAXN];
int cnt[][][MAXN],mx[][],d[][];

inline int query(int l,int r)
{
        int i;
        static int sum[MAXN];
        int mx = ,ret = ;
        int p = pos[l],q = pos[r];
        if (p == q)
        {
                for (i = l; i <= r; i++)
                {
                        sum[rk[i]]++;
                         if (sum[rk[i]] > mx || (sum[rk[i]] == mx && rk[i] < ret))
                         {
                                 mx = sum[rk[i]];
                                 ret = rk[i];
                         }
                }
                for (i = l; i <= r; i++) sum[rk[i]]--;
                return ret;
        }    else
        {
                ret = d[p+][q-];
                mx = cnt[p+][q-][ret];
                for (i = l; i <= R[p]; i++)
                {
                        cnt[p+][q-][rk[i]]++;
                        if (cnt[p+][q-][rk[i]] > mx || (cnt[p+][q-][rk[i]] == mx && rk[i] < ret))
                        {
                                mx = cnt[p+][q-][rk[i]];
                                ret = rk[i];
                        }
                 }
                for (i = L[q]; i <= r; i++)
                {
                        cnt[p+][q-][rk[i]]++;
                        if (cnt[p+][q-][rk[i]] > mx || (cnt[p+][q-][rk[i]] == mx && rk[i] < ret))
                        {
                                mx = cnt[p+][q-][rk[i]];
                                ret = rk[i];
                        }
                }
                for (i = l; i <= R[p]; i++) cnt[p+][q-][rk[i]]--;
                for (i = L[q]; i <= r; i++) cnt[p+][q-][rk[i]]--;
                return ret;
        }
} 

int main()
{

        scanf("%d%d",&n,&m);
        for (i  = ; i <= n; i++)
        {
                scanf("%d",&a[i]);
                val[++len] = a[i];
        }
        sort(val+,val+len+);
        len = unique(val+,val+len+) - val - ;
        for (i = ; i <= n; i++) rk[i] = lower_bound(val+,val+len+,a[i]) - val;
        block = (int)pow(n*1.0,1.0/);
        if (block) t = n / block;
        for (i = ; i <= block; i++)
        {
                L[i] = (i - ) * t + ;
                R[i] = i * t;
        }
        if (R[block] < n)
        {
                block++;
                L[block] = R[block-] + ;
                R[block] = n;
        }
        for (i = ; i <= block; i++)
        {
                for (j = L[i]; j <= R[i]; j++)
                {
                        pos[j] = i;
                }
        }
        for (i = ; i <= block; i++)
        {
                for (j = i; j <= block; j++)
                {
                        for (k = L[i]; k <= R[j]; k++)
                                cnt[i][j][rk[k]]++;
                        for (k = ; k <= len; k++)
                        {
                                if (cnt[i][j][k] > mx[i][j] || (cnt[i][j][k] == mx[i][j] && k < d[i][j]))
                                {
                                        mx[i][j] = cnt[i][j][k];
                                        d[i][j] = k;
                                }
                        }
                }
        }
        lastans = ;
        for (i = ; i <= m; i++)
        {
                scanf("%d%d",&l0,&r0);
                l = (l0 + lastans - ) % n + ;
                r = (r0 + lastans - ) % n + ;
                if (l > r) swap(l,r);
                printf("%d\n",lastans = val[query(l,r)]);
        }

        return ;

}