蒲公英——分块

本文介绍了一种通过将数列分块来优化操作复杂度的方法,适用于处理区间查询和更新问题。通过将数列分为sqrt(n)个块,可以将操作复杂度从O(mn)降低到O(sqrt(n))。文章详细解释了如何实现分块,并提供了一个具体实例,展示了如何使用分块技巧求解区间众数问题。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

分块基本描述
一个数列进行操作是O(mn)O(mn)的,那么每n−−√n个分成一段,共分成n−−√n段,然后我们操作时,对于在一段内的直接单点暴力修改,跨段部分按段转移,这样每次操作复杂度就降成了O(n−−√)O(n)。
非常简单的思想,糟心的是如何实现。

建块
首先我们先求出块的数目bct = (int)sqrt(n)
然后用block数组记录某一个元素在第几块,实现:
block[i] = (i - 1) / bct + 1
然后就完成了

解析
题意就是强制在线求区间众数

首先你要知道一个性质:维护一个集合的众数,这个众数要么是当前集合内出现次数最多的那个数,要么存在于新加入的那些数中

先离散化,然后分块

预处理f[i][j]表示从第i块到第j块的众数

之后对于每个离散化后的数开vector,记录下出现的位置

查询的时候,对于整块,直接可以二分查找

两侧的零散块,对于出现的每一个数字,都需要二分判断是不是更新的众数

代码:
1.时间超限 不知道为什么

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+10;
const int maxx=1e3+10;

int n,m,a[maxn],b[maxn],f[maxx][maxx];
vector<int>c,v[maxn];
int s[maxn];

int q=200;


int cnt[maxn];
int maxm,ans;

int pre(int x)
{
    memset(cnt,0,sizeof(cnt));
    maxm=0,ans=0;
    for(int i=(x-1)*q+1;i<=n;i++)
    {
        int tmp=b[i];
        cnt[a[i]]++;
        if(maxm<cnt[a[i]]||(maxm==cnt[a[i]]&&a[i]<ans))
        {
            maxm=cnt[a[i]],ans=a[i];
        }
        f[x][tmp]=ans;
    }
}
int num(int l,int r,int x)
{
    int x1=lower_bound(v[x].begin(),v[x].end(),l)-v[x].begin();
    int x2=upper_bound(v[x].begin(),v[x].end(),r)-v[x].begin();
    return x2-x1;
}
int query(int l,int r)
{
   int tmp=min(b[l]*q,r);
   int ans=f[b[l]+1][b[r]-1];
   int maxm=num(l,r,ans);
   for(int i=l;i<=tmp;i++)
   {
       int d=num(l,r,a[i]);
       if(maxm<d||(maxm==d&&a[i]<ans))
            maxm=d,ans=a[i];

   }
   for(int i=(b[r]-1)*q+1;i<=r;i++)
   {
       int d = num(l, r, a[i]);
        if(maxm < d || (maxm == d && a[i] < ans) )
            maxm = d, ans = a[i];
   }
   return s[ans];

}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        c.push_back(a[i]);
        b[i]=(i-1)/q+1;
    }

    sort(c.begin(),c.end());
    unique(c.begin(),c.end());
    for(int i=1;i<=n;i++)
    {
        int tmp=a[i];
        a[i]=lower_bound(c.begin(),c.end(),a[i])-c.begin()+1;
        s[a[i]]=tmp;
    }

    for(int i=1;i<=n;i++)
        v[a[i]].push_back(i);

    for(int i=1;i<=b[n];i++)
    {
        pre(i);
    }
    int x=0;
    for(int i=1;i<=m;i++)
    {
        int l0,r0;
        scanf("%d%d",&l0,&r0);
        int l=(l0+x-1)%n+1,r=(r0+x-1)%n+1;
        if(l>r)
            swap(l,r);
        x=query(l,r);
        printf("%d\n",x);
    }
}

2.能过得代码,思路差不多

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
int maxn,m,n,tot,T,mode[40][40],last[40010],
a[40010],ord[40010],L[40],R[40],num[40][40][40010],cnt[40010];
int main()
{
    int i,j,k,p,q,x,y,z,K,ans,l,r,ll,rr;
    scanf("%d%d",&n,&m);
    for (i=1;i<=n;i++)
      scanf("%d",&a[i]),ord[i]=a[i];
    sort(ord+1,ord+n+1);
    maxn=unique(ord+1,ord+n+1)-ord-1;
    for (i=1;i<=n;i++)
      a[i]=lower_bound(ord+1,ord+maxn+1,a[i])-ord;
    T=pow(n+0.1,2.0/3);
    tot=pow(n,1.0/3);
    for (i=1;i<=tot;i++)
    {
        L[i]=R[i-1]+1;
        if (i==tot) R[i]=n;
        else R[i]=L[i]+T-1;
    }
    for (i=1;i<=tot;i++)
    {
        for (j=L[i];j<=R[i];j++)
          num[i][i][a[j]]++;
        for (j=1;j<=maxn;j++)
          if (num[i][i][j]>num[i][i][mode[i][i]])
            mode[i][i]=j;
    }
    for (i=tot;i;i--)
      for (j=i+1;j<=tot;j++)
        for (k=1;k<=maxn;k++)
        {
            num[i][j][k]=num[i][i][k]+num[i+1][j][k];
            if (num[i][j][k]>num[i][j][mode[i][j]])
              mode[i][j]=k;
        }
    ans=0;
    for (K=1;K<=m;K++)
    {
        scanf("%d%d",&l,&r);
        l=(l+ans-1)%n+1;
        r=(r+ans-1)%n+1;
        if (l>r) swap(l,r);
        ll=1;
        while (ll<=tot&&L[ll]<l) ll++;
        rr=tot;
        while (rr&&R[rr]>r) rr--;
        if (ll>rr)
        {
            ans=0;
            for (i=l;i<=r;i++)
            {
                if (last[a[i]]<K)
                {
                    last[a[i]]=K;
                    cnt[a[i]]=1;
                }
                else
                  cnt[a[i]]++;
                if (cnt[a[i]]>cnt[ans]||(cnt[a[i]]==cnt[ans]&&a[i]<ans)) ans=a[i];
            }
            ans=ord[ans];
            printf("%d\n",ans);
            continue;
        }
        ans=mode[ll][rr];
        last[ans]=K;
        cnt[ans]=num[ll][rr][ans];
        for (i=l;i<L[ll];i++)
        {
            if (last[a[i]]<K)
            {
                last[a[i]]=K;
                cnt[a[i]]=num[ll][rr][a[i]]+1;
            }
            else
              cnt[a[i]]++;
            if (cnt[a[i]]>cnt[ans]||(cnt[a[i]]==cnt[ans]&&a[i]<ans)) ans=a[i];
        }
        for (i=R[rr]+1;i<=r;i++)
        {
            if (last[a[i]]<K)
            {
                last[a[i]]=K;
                cnt[a[i]]=num[ll][rr][a[i]]+1;
            }
            else
              cnt[a[i]]++;
            if (cnt[a[i]]>cnt[ans]||(cnt[a[i]]==cnt[ans]&&a[i]<ans)) ans=a[i];
        }
        ans=ord[ans];
        printf("%d\n",ans);
    }
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值