2724: [Violet 6]蒲公英

题目大意：不带修改的区间众数。

%%%clj大大：区间众数

分块，首先有个重要结论：

摘自clj论文。

分块做。

那么我们就预处理f[i][j]，第i块~j块间的众数。

每次枚举最左边（右边）的数，与剩下的块的众数（已求出）比较一下谁更多，更新数组。（类似区间dp）

至于求区间内x出现次数，可以用数据结构，也可以二分下标求解。

int check(int x,int y,int c)
{
	int l=0,r=v[c].size()-1,ans;
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(v[c][mid]>=x) ans=mid,r=mid-1;
		else l=mid+1;
	}
	l=0,r=v[c].size()-1;int ans1;
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(v[c][mid]<=y) ans1=mid,l=mid+1;
		else r=mid-1;
	}
	return ans1-ans+1;
}

那么最后求答案就是分块模板了。

code：

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
#include<vector>
#include<cmath>
#include<algorithm>
using namespace std;
vector<int> v[40010];
int n,m,e[40010],f[210][210],L[210],R[210],pre[40010];
struct node1{
	int a,b,num;
}a[40010];
struct node{
	int a,num;
};
bool cmp(node1 a,node1 b){return a.b<b.b;}
int check(int x,int y,int c)
{
	int l=0,r=v[c].size()-1,ans;
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(v[c][mid]>=x) ans=mid,r=mid-1;
		else l=mid+1;
	}
	l=0,r=v[c].size()-1;int ans1;
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(v[c][mid]<=y) ans1=mid,l=mid+1;
		else r=mid-1;
	}
	return ans1-ans+1;
}
node findmax(int x,int y,int X,int Y)
{
	node ans;ans.num=0;
	for(int i=x;i<=y;i++)
	{
		int tmp=check(X,Y,a[i].a);
		if(ans.num<tmp||ans.num==tmp&&ans.a>a[i].a) ans.a=a[i].a,ans.num=tmp;
	}
	return ans;
}
int main()
{
	scanf("%d %d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&a[i].b),a[i].num=i;
	sort(a+1,a+n+1,cmp);
	int tmp=0;
	for(int i=1;i<=n;i++)
	{
		if(i==1||a[i].b!=a[i-1].b) tmp++;
		a[a[i].num].a=tmp;
		e[tmp]=a[i].b;
	}
	for(int i=1;i<=tmp;i++) v[i].push_back(0);
	for(int i=1;i<=n;i++) v[a[i].a].push_back(i);
	tmp=sqrt(n);
	for(int i=1;i<=n;i++)
	{
		pre[i]=i/tmp+1;
		if(pre[i]!=pre[i-1]) L[pre[i]]=i,R[pre[i-1]]=i-1;
	}
	R[pre[n]]=n;
	for(int l=1;l<=pre[n];l++)
		for(int i=pre[1];i<=pre[n];i++)
		{
			int j=i+l-1;if(j>pre[n]) break;
			node MAX=findmax(L[i],R[i],L[i],R[j]);
			if(i==j) f[i][j]=MAX.a;
			else
			{
				tmp=check(L[i],R[j],f[i+1][j]);
				if(MAX.num>tmp||MAX.num==tmp&&MAX.a<f[i+1][j]) f[i][j]=MAX.a;
				else f[i][j]=f[i+1][j];
			}
		}
	node ans;ans.a=0;
	for(int i=1;i<=m;i++)
	{
		int l,r;scanf("%d %d",&l,&r);
		l=(l+e[ans.a]-1)%n+1;r=(r+e[ans.a]-1)%n+1;
		if(l>r) swap(l,r);
		if(pre[l]==pre[r]) ans=findmax(l,r,l,r);
		else
		{
			int x=R[pre[l]],y=L[pre[r]];
			node tmp1=findmax(l,x,l,r),tmp2=findmax(y,r,l,r);
			if(tmp1.num>tmp2.num||tmp1.num==tmp2.num&&tmp1.a<tmp2.a) ans=tmp1;
			else ans=tmp2;
			if(pre[l]+1<=pre[r]-1)
			{
				tmp1.a=f[pre[l]+1][pre[r]-1];tmp1.num=check(l,r,tmp1.a);
				if(tmp1.num>ans.num||tmp1.num==ans.num&&tmp1.a<ans.a) ans=tmp1;
			}
		}
		printf("%d\n",e[ans.a]);
	}
}