分块基本描述
一个数列进行操作是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);
}
}