在线求区间众数和离线求区间众数(uestc的题 冬马和纱天下第一和 我永远喜欢冬马和纱)

两种做法都是要用到分块

在线的题的地址是  https://acm.uestc.edu.cn/problem/dong-ma-he-sha-tian-xia-di-yi/description

题意简化就是求一个区间内的最大众数,本来是看了网上的在线区间众数求法(https://blog.csdn.net/BerryKanry/article/details/75478987?locationNum=4&fps=1),复杂度是nsqrt(n)log(n),但竟然t了,因为出题人专门把这种做法给卡了,正解的复杂度nsqrt(n)而已,接下来讲下做法(来自该题的题解)

我们用???[?][?]表示前?个块中j这个数出现的次数,注意,因为我们按照√?来分 块,所以?最多是200,而?最大是40000,因此总的空间大小是8e6,空间是完全 足够的。实际上???[?][?]相当于j这个数的前缀和,我们只需要先将每个数放在自 己所在的块中,然后暴力扫描每个数,然后做前缀和就好了,总的复杂度可以看 出来是? ?所以在时间上没有任何问题。具体实现希望大家认真思考一下。  然后我们用???[?][?]表示第?个块到第?个块的区间众数是哪个数。我们先枚举左端 点?,然后从?开始枚举到?,中途中我们只需要一边统计每个数并且更新到达每个 区间右端点的区间众数就好了,这部分可能有点复杂,但是我认为稍加思索并不 困难。(再强调一下我是直接把该题的题解复制来的,不是原创,原创的只是看了这段话后写出的代码)

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;

const int maxn = 4e4+9;

int a[maxn],b[maxn];
int n,m,num,len,kuai;

int sum[209][maxn];
int xx[209][209];
int cnt[maxn];

void init(){
    len=sqrt(n);
    kuai=len*len==n?len:len+1;
    for(int i=1;i<=kuai;i++){
        for(int j=0;j<=num;j++)cnt[j]=0;
        int ma=0,my=0;
        for(int j=(i-1)*len+1;j<=n;j++){
            cnt[a[j]]++;
            if(my<cnt[a[j]]){
                my=cnt[a[j]];
                ma=a[j];
            }
            xx[i][(j-1)/len+1]=ma;
        }
    }
    for(int i=1;i<=kuai;i++){
        for(int j=(i-1)*len+1;j<=i*len;j++){
            sum[i][a[j]]++;
        }
    }
    for(int i=2;i<=kuai;i++){
        for(int j=1;j<=num;j++){
            sum[i][j]+=sum[i-1][j];
        }
    }
    for(int i=0;i<=num;i++)cnt[i]=0;
}

int query(int l,int r){
    int nl=(l-1)/len+2,nr=(r-1)/len;
    int mx,my=0;
    if(nl>nr){
        for(int i=l;i<=r;i++){
            cnt[a[i]]++;
            my=max(cnt[a[i]],my);
        }
        for(int i=l;i<=r;i++){
            cnt[a[i]]--;
        }
        return my;
    }

    mx=xx[nl][nr];
    my=sum[nr][mx]-sum[nl-1][mx];

    for(int i=l;i<=min(r,(nl-1)*len);i++){
        cnt[a[i]]++;
//        cout<<sum[nr][a[i]]-sum[nl-1][a[i]]<<" "<<nl<<" "<<nr<<endl;
        if(my<sum[nr][a[i]]-sum[nl-1][a[i]]+cnt[a[i]]){
            my=sum[nr][a[i]]-sum[nl-1][a[i]]+cnt[a[i]];
        }
    }
    for(int i=max(l,(nr*len)+1);i<=r;i++){
        cnt[a[i]]++;
        if(my<sum[nr][a[i]]-sum[nl-1][a[i]]+cnt[a[i]]){
            my=sum[nr][a[i]]-sum[nl-1][a[i]]+cnt[a[i]];
        }
    }
    for(int i=l;i<=min(r,(nl-1)*len);i++){
        cnt[a[i]]--;
    }
    for(int i=max(l,(nr*len)+1);i<=r;i++){
        cnt[a[i]]--;
    }
    return my;
}

int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        b[i]=a[i];
    }
    sort(b+1,b+1+n);
    num=unique(b+1,b+1+n)-(b+1);
    for(int i=1;i<=n;i++){
        a[i]=lower_bound(b+1,b+1+num,a[i])-b;
    }
    init();
    int l,r,x=0;
    for(int i=0;i<m;i++){
        scanf("%d%d",&l,&r);
        l=(l+x-1)%n+1,r=(r+x-1)%n+1;
        if(l>r)swap(l,r);
        x=query(l,r);
        printf("%d\n",x);
    }
    return 0;
}

 

然后就是离线的做法了,题目衔接:https://acm.uestc.edu.cn/problem/wo-yong-yuan-xi-huan-dong-ma-he-sha/description

一个很普通的莫队,就不讲解,直接看代码,这个是一发过的,然后在线就t成鬼了

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;

const int maxn = 2e5+9;

int cnt[maxn],num[maxn];
int a[maxn],b[maxn];
int ans[maxn];
int n,m;
int kuai;
struct rt{
    int id;
    int l,r;
}arr[maxn];

bool cmp(rt a,rt b){
    if(a.l/kuai==b.l/kuai)return a.r<b.r;
    return a.l/kuai<b.l/kuai;
}

int main(){
    scanf("%d%d",&n,&m);
    kuai=sqrt(n);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        b[i]=a[i];
    }
    sort(b+1,b+1+n);
    int tot=unique(b+1,b+1+n)-(b+1);
    for(int i=1;i<=n;i++){
        a[i]=lower_bound(b+1,b+1+tot,a[i])-b;
    }
    for(int i=0;i<m;i++){
       scanf("%d%d",&arr[i].l,&arr[i].r);
       arr[i].id=i+1;
    }
    sort(arr,arr+m,cmp);
    int ma=-INF;
    int l=1,r=0;
    for(int i=0;i<m;i++){
//        cout<<arr[i].l<<" "<<arr[i].r<<endl;
        while(l<arr[i].l){
            num[cnt[a[l]]]--;
            if(ma==cnt[a[l]]&&!num[cnt[a[l]]]){
                ma=cnt[a[l]]-1;
            }
            cnt[a[l]]--;
            num[cnt[a[l]]]++;
            l++;
        }
        while(r<arr[i].r){
            ++r;
            num[cnt[a[r]]]--;
            cnt[a[r]]++;
            num[cnt[a[r]]]++;
            if(ma<cnt[a[r]])ma=cnt[a[r]];
        }
        while(l>arr[i].l){
            --l;
            num[cnt[a[l]]]--;
            cnt[a[l]]++;
            num[cnt[a[l]]]++;
            if(ma<cnt[a[l]])ma=cnt[a[l]];
        }
        while(r>arr[i].r){
            num[cnt[a[r]]]--;
            if(ma==cnt[a[r]]&&!num[cnt[a[r]]]){
                ma=cnt[a[r]]-1;
            }
            cnt[a[r]]--;
            num[cnt[a[r]]]++;
            r--;
        }
        ans[arr[i].id]=ma;
    }
    for(int i=1;i<=m;i++){
        printf("%d\n",ans[i]);
    }
    return 0;
}

 

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值