【分块】【Violet】蒲公英

【描述】
亲爱的哥哥:
你在那个城市里面过得好吗?
我在家里面最近很开心呢。昨天晚上奶奶给我讲了那个叫「绝望」的大坏蛋的故事的说!它把人们的房子和田地搞坏,还有好多小朋友也被它杀掉了。我觉得把那么可怕的怪物召唤出来的那个坏蛋也很坏呢。不过奶奶说他是很难受的时候才做出这样的事的……
最近村子里长出了一大片一大片的蒲公英。一刮风,这些蒲公英就能飘到好远的地方了呢。我觉得要是它们能飘到那个城市里面,让哥哥看看就好了呢!
哥哥你要快点回来哦!
爱你的妹妹 Violet
Azure 读完这封信之后微笑了一下。
“蒲公英吗……” 在乡下的小路旁种着许多蒲公英,而我们的问题正是与这些蒲公英有关。
为了简化起见,我们把所有的蒲公英看成一个长度为n的序列 ( a1 , a2 …… an ) ,其中ai 为一个正整数,表示第i棵蒲公英的种类编号。
而每次询问一个区间 [l,r],你需要回答区间里出现次数最多的是哪种蒲公英,如果有若干种蒲公英出现次数相同,则输出种类编号最小的那个。
注意,你的算法必须是在线的
【输入】
第一行两个整数n,m ,表示有n株蒲公英,m 次询问。
接下来一行n个空格分隔的整数ai,表示蒲公英的种类
再接下来m 行每行两个整数l,r,我们令上次询问的结果为 x(如果这是第一次询问, 则 x=0)。
令 l=( l + x - 1) mod n+1 ,r =( r + x - 1 ) mod n+1,如果 l>r,则交换 l,r 。
最终的询问区间为[l,r]。
【输出】
输出m 行。每行一个整数,表示每次询问的结果。

【思路】

我们考虑分块。
1.处理 1 i 1-i 块中每个数出现了多少次。
2.处理 1 i 1-i 块中的众数。
3.查询:
答案要么是中间整块的答案,要么在两边出现过。另开一个数组维护两边每个数出现的次数,再加上中间整块的出现次数即可比较得到答案。注意,这个维护两边每个数出现次数的数组不能直接memset。不然复杂度就错了,要逐一删除以清空。其他细节见代码。
代码:(自以为写得不错 )

#include<bits/stdc++.h>
#include<tr1/unordered_map>
#define re register
using namespace std;
const int N=4e4+5;
tr1::unordered_map<int,int>lsr;
inline int red(){
    int data=0;bool w=0; char ch=getchar();
    while(ch!='-' && (ch<'0' || ch>'9')) ch=getchar();
    if(ch=='-') w=1,ch=getchar();
    while(ch>='0' && ch<='9') data=(data<<3)+(data<<1)+ch-'0',ch=getchar();
    return w?-data:data;
}
int n,m,a[N],b[N],c,tot=0,l,r,las=0,mx,id,L=1,R=0;
int cnt[205][N],col[N],f[205][205];
int num[N];
inline void check(const int&val){
	int tot=num[val]+cnt[R][val]-cnt[L-1][val];
	if(tot>mx||(mx==tot&&val<id))id=val,mx=tot;
}
inline void get(int l,int r){
	for(int re i=l;i<=r;++i)
	++num[a[i]],check(a[i]);
}
inline void Del(int l,int r){
	for(int re i=l;i<=r;++i)--num[a[i]];
}
int main(){
	n=red();m=red();int S=sqrt(n);
	for(int re i=1;i<=n;i++)b[i]=a[i]=red(),col[i]=(i-1)/S+1;
	sort(b+1,b+n+1);tot=unique(b+1,b+n+1)-b-1;
	for(int re i=1;i<=tot;i++)lsr[b[i]]=i;
	for(int re i=1;i<=n;i++)a[i]=lsr[a[i]];
	for(int re i=1;i<=col[n];i++)
		for(int re j=1;col[j]<=i&&j<=n;++j)
			++cnt[i][a[j]];
	for(int re i=1;i<=col[n];i++){
		int now=(i-1)*S;mx=-1,id=1e9;memset(num,0,sizeof(int)*(tot+1));
		for(int re j=i;j<=col[n];++j){
			while(col[++now]==j&&now<=n)++num[a[now]],check(a[now]);
			f[i][j]=id;--now;
		}
	}memset(num,0,sizeof(num)); 
	while(m--){
		l=(red()+las-1)%n+1,r=(red()+las-1)%n+1;
		if(r<l)r^=l^=r^=l;
		mx=-1,id=1e9,L=col[l]+1,R=col[r]-1;
		if(L>R)L=R+1;
		get(l,min(r,col[l]*S));
		if(col[l]^col[r]){
			get((col[r]-1)*S+1,r);
			check(f[L][R]);
			Del((col[r]-1)*S+1,r);
		}cout<<(las=b[id])<<"\n";
		Del(l,min(r,col[l]*S));
	}return 0;
}
发布了106 篇原创文章 · 获赞 22 · 访问量 5485

猜你喜欢

转载自blog.csdn.net/weixin_44111457/article/details/102655411