【分塊】【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.處理1i1-i塊中每個數出現了多少次。
2.處理1i1-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;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章