我的新博客: https://blog.zcmimi.top

已搬迁至新博客,此博客暂停更新

×

[Violet]蒲公英

题意:

给出一个长度为 \(n\) 序列\(a\)\(m\)次询问,每次询问区间 \(l,r\) 里的众数(出现次数最多的数)。若有多个,输出最小的

\(a_i \leq 10^9, n \leq 40000, m \leq 50000\),强制在线。

题解:

看了题解才懂的。根据https://www.cnblogs.com/acfunction/p/10051345.html

hzwer给出了更巧妙的方法http://hzwer.com/3582.html

\(a_i \le 10^9\),所以先离散化。

分块大法好!!!

分析一下:

求出现次数最多的数,那就预处理一下次数,\(s[i][j]\)表示离散化后\(i\)个块中\(j\)出现了多少次

方法:类似前缀和

然后为了提高效率,再预处理p数组,\(p[i][j]\)表示离散化后\(i\)个块到第\(j\)个块(最小的)众数

方法:枚举\(i,j\),开个桶,记录每个数出现的次数,然后去出现次数最多的(数值最小的)。

接下来:分类讨论

  1. \(pr\)(\(r\)所在的块)-\(pl\)(\(l\)所在的块)$\geq 2 $

    暴力枚举,和处理\(p\)数组时差不多。

  2. \(pr-pl>2\)

    枚举块外,块内的已经预处理好了。

    \(每个数的出现次数=块内的(已预处理)+块外的(暴力)\)

总的效率为$ O(2n \sqrt n + m \sqrt n) $

要注意细节比较多,QwQ。

代码有点啰嗦,长了点

#include<bits/stdc++.h>
#define Fur(i,x,y) for(int i=x;i<=y;i++)
#define clr(x,y) memset(x,y,sizeof(x))
int MIN(int x,int y){return x<=y?x:y;}
void SWAP(int &x,int &y){x^=y;y^=x;x^=y;}
using namespace std;
#define N 40040
int n,m,L,len,s[220][N],val[N],b[N];
bool vis[N];
struct nod{int id,val,s;}a[N];
struct node{int val,s;}p[220][220];
bool cmp1(nod x,nod y){return x.val<y.val;}
bool cmp2(nod x,nod y){return x.id<y.id;}
void pre(){
    Fur(i,1,len){
        clr(b,0);
        node tmp;tmp.val=tmp.s=0;
        Fur(j,i,len){
            Fur(k,(j-1)*L+1,MIN(j*L,n)){
                b[a[k].s]++;
                if(b[a[k].s]>tmp.s)tmp.val=a[k].s,tmp.s=b[a[k].s];
                else if(b[a[k].s]==tmp.s)tmp.val=MIN(tmp.val,a[k].s);
            }
            p[i][j]=tmp;
        }
    }
    Fur(i,1,len){
        Fur(j,1,n)s[i][a[j].s]=s[i-1][a[j].s];
        Fur(j,(i-1)*L+1,MIN(n,i*L))s[i][a[j].s]++;
    }
}
int main(){
    ios::sync_with_stdio(false);
    cin>>n>>m;L=sqrt(n);len=(n+L-1)/L;
    Fur(i,1,n)cin>>a[i].val,a[i].id=i;
    sort(a+1,a+n+1,cmp1);a[0].val=-1;
    Fur(i,1,n){
        a[i].s=a[i-1].s+(a[i-1].val!=a[i].val);
        val[a[i].s]=a[i].val;
    }
    sort(a+1,a+n+1,cmp2);
    pre();
    int l,r,la=0,ans=0;
    Fur(i,1,m){
        cin>>l>>r;
        l=(l+la-1)%n+1;r=(r+la-1)%n+1;
        if(l>r)SWAP(l,r);
        int pl=(l-1)/L+1,pr=(r-1)/L+1;
        if(pr-pl<=2){
            ans=0;
            Fur(j,l,r)b[a[j].s]=0;
            Fur(j,l,r){
                b[a[j].s]++;
                if(b[a[j].s]>b[ans])ans=a[j].s;
                else if(b[a[j].s]==b[ans])ans=MIN(ans,a[j].s);
            }
        }
        else{
            ans=p[pl+1][pr-1].val;
            b[ans]=vis[ans]=0;
            Fur(j,l,MIN(n,pl*L))b[a[j].s]=vis[a[j].s]=0;
            Fur(j,(pr-1)*L+1,r)b[a[j].s]=vis[a[j].s]=0;
            Fur(j,l,MIN(n,pl*L))b[a[j].s]++;
            Fur(j,(pr-1)*L+1,r)b[a[j].s]++;
            int mx=0,mxv;
            Fur(j,l,MIN(n,pl*L))if(!vis[a[j].s]){
                vis[a[j].s]=1;
                int tmp=b[a[j].s]+s[pr-1][a[j].s]-s[pl][a[j].s];
                if(mx<tmp)mx=tmp,mxv=a[j].s;
                else if(mx==tmp)mxv=MIN(mxv,a[j].s);
            }
            Fur(j,(pr-1)*L+1,r)if(!vis[a[j].s]){
                vis[a[j].s]=1;
                int tmp=b[a[j].s]+s[pr-1][a[j].s]-s[pl][a[j].s];
                if(mx<tmp)mx=tmp,mxv=a[j].s;
                else if(mx==tmp)mxv=MIN(mxv,a[j].s);
            }
            if(mx>b[ans]+p[pl+1][pr-1].s)ans=mxv;
            else if(mx==b[ans]+p[pl+1][pr-1].s)ans=MIN(ans,mxv);
        }
        cout<<(la=val[ans])<<"\n";
    }
}
posted @ 2019-01-25 21:38  zcmimi  阅读(259)  评论(0编辑  收藏