[bzoj2724]蒲公英——分块

题目大意:

给定一个序列和若干次查询,求若干次查询的区间的众数。要求强制在线。

思路:

一道分块的模板题,考虑预处理出来第每两个块之间的众数,然后发现区间的众数只有可能是整数块的众数或者是两边不完整的块的众数,旁边两个块的每一个数暴力求次数然后再和中间的众数的次数相比较即可。
考虑如何暴力求制定区间内某数出现的次数,可以对于每一种数开一个vector记下它每次出现的位置,这样的话在vector上二分即可。但是这种做法的时间复杂度为 O(n1.5logn) O ( n 1.5 log ⁡ n ) ,不是特别地优秀。
于是可以换一种求法:记下每一个数在每一个块内出现的次数,然后两边暴力求就好了(这样的话这题就用了两个分块了。。。)

/*=========================
* Author : ylsoi
* Problem : pugongyin
* Algrithm : Block
* Time : 2018.5.12
* =======================*/
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<climits>
#include<vector>
using namespace std;
void File(){
    freopen("[bzoj2724].in","r",stdin);
    freopen("[bzoj2724].out","w",stdout);
}
#define REP(i,a,b) for(register int i=a;i<=b;++i)
#define DREP(i,a,b) for(register int i=a;i>=b;--i)
#define MREP(i,x) for(register int i=beg[x];i;i=E[i].last)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define inf INT_MAX
#define PII pair<int,int>
template<typename T>bool chkmax(T &_,T __){return _<__ ? (_=__,1) : 0;}
const int maxn=4e4+10;
const int maxb=300+10;
int n,m,a[maxn],id[maxn],cnt_id,va[maxn],most[maxb][maxb];
int block,be[maxn],tot;
vector<int>pos[maxn];
PII b[maxn];
bool cmp(PII x,PII y){
    if(x.first!=y.first)return x.first<y.first;
    else return x.second<y.second;
}
int Count(int l,int r,int x){
    l=lower_bound(pos[x].begin(),pos[x].end(),l)-pos[x].begin();
    r=upper_bound(pos[x].begin(),pos[x].end(),r)-pos[x].begin()-1;
    return r-l+1;
}
void init(){
    scanf("%d%d",&n,&m);
    block=pow(n,0.5);
    tot=ceil(n*1.0/block);
    REP(i,1,n)be[i]=ceil(i*1.0/block);
    REP(i,1,n)scanf("%d",&a[i]);
    REP(i,1,n)b[i].first=a[i],b[i].second=i;
    sort(b+1,b+n+1,cmp);
    REP(i,1,n){
        if(b[i].first!=b[i-1].first)
            va[++cnt_id]=b[i].first;
        id[b[i].second]=cnt_id;
        pos[cnt_id].push_back(b[i].second);
    }
    int len=ceil(n*1.0/block);  
    REP(i,1,len){
        REP(j,1,len-i+1){
            if(len>2)most[j][j+i-1]=most[j+1][j+i-2];
            int l1=(j-1)*block+1,r1=min(n,j*block);
            int l2=(j+i-2)*block+1,r2=min(n,(j+i-1)*block);
            int Max=len>2 ? Count(l1,r2,most[j][j+i-1]) : 0;
            REP(k,l1,r1){
                int num=id[k],aa=Count(l1,r2,num);
                if(aa>Max || (aa==Max && num<most[j][j+i-1])){
                    most[j][j+i-1]=num;
                    Max=aa;
                }
            }
            REP(k,l2,r2){
                int num=id[k],aa=Count(l1,r2,num);
                if(aa>Max || (aa==Max && num<most[j][j+i-1])){
                    most[j][j+i-1]=num;
                    Max=aa;
                }
            }
            //cout<<j<<" "<<j+i-1<<" "<<va[most[j][j+i-1]]<<endl;
        }
    }
}
int query(int l,int r){ 
    int L=be[l],R=be[r];
    int Max=L+1<R ? Count(l,r,most[L+1][R-1]) : 0;
    int ret=L+1<R ? most[L+1][R-1] : 0;
    REP(k,l,min(L*block,n)){
        int num=id[k],aa=Count(l,r,num);
        if(aa>Max || (aa==Max && num<ret)){
            ret=num;
            Max=aa;
        }
    }
    REP(k,(R-1)*block+1,r){
        int num=id[k],aa=Count(l,r,num);
        if(aa>Max || (aa==Max && num<ret)){
            ret=num;
            Max=aa;
        }
    }
    return ret;
}
int main(){
    //File();
    init();
    int last_ans=0;
    REP(i,1,m){
        int l,r;
        scanf("%d%d",&l,&r);
        l=(l+last_ans-1)%n+1;
        r=(r+last_ans-1)%n+1;
        if(l>r)swap(l,r);
        last_ans=va[query(l,r)];
        printf("%d\n",last_ans);
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值