今天全是奇妙的分块大法。
【BZOJ 2714】蒲公英
给一个序列,求一个区间的众数。(强制在线)
【题解】
先分块。
设f[i][j]表示第i块到第j块总的众数。
然后每次对于一个询问[l,r],x到y已经知晓,答案只可能是l~x,y~r,和块x~y的众数之一,所以接下来要做的就是判断l~x和y~r中存不存在数字出现次数比块x到y众数出现更多的数。
不过怎么判断呢?可以先用vector来保存每个数字出现的位置(因为本题的值域较大,所以要离散化一下,就是map来搞一搞),在用二分答案查找这个数字在[l,r]内出现的次数即可。
详见代码。
【代码】
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<map>
#define N 40010
#define B 200
using namespace std;
map<int,int> mp;
vector<int> pos[N];
int cnt,val[N],a[N];
int belong[N],tot[N],f[B+5][B+5];
int n;
void pre(int x)
{
int now = 0,maxnsum = 0;
memset(tot,0,sizeof(tot));
for(int i = (x-1)*B + 1;i <= n;i++) {
tot[a[i]] ++;
if(tot[a[i]] > maxnsum || (tot[a[i]] == maxnsum && val[now] > val[a[i]]) )
now = a[i],maxnsum = tot[a[i]];
f[x][belong[i]] = now;
}
}
int sigma(int l,int r,int x) {
return upper_bound(pos[x].begin(),pos[x].end(),r)-lower_bound(pos[x].begin(),pos[x].end(),l);
}
int query(int l,int r)
{
int now = f[belong[l]+1][belong[r]-1],maxnsum = sigma(l,r,now);
for(int i = l;i <= min(r,belong[l]*B);i++) {
int tmp = sigma(l,r,a[i]);
if(tmp > maxnsum || (tmp == maxnsum && val[now] > val[a[i]]) )
now = a[i],maxnsum = tmp;
}
if(belong[l] != belong[r]) {
for(int i = (belong[r]-1)*B+1;i <= r;i++) {
int tmp = sigma(l,r,a[i]);
if(tmp > maxnsum || (tmp == maxnsum && val[now] > val[a[i]]) )
now = a[i],maxnsum = tmp;
}
}
return now;
}
int main()
{
int m,l,r;
int ans = cnt = 0;
scanf("%d%d",&n,&m);
for(int i = 1;i <= n;i++) {
scanf("%d",&a[i]);
if(!mp[a[i]]) {
mp[a[i]] = ++cnt; val[cnt] = a[i];
}
a[i] = mp[a[i]];
pos[a[i]].push_back(i);
}
for(int i = 1;i <= n;i++) belong[i] = (i-1)/B + 1;
for(int i = 1;i <= belong[n];i++) pre(i);
for(int i = 1;i <= m;i++) {
scanf("%d%d",&l,&r);
l = (l+ans-1) % n + 1; r = (r+ans-1) % n + 1;
if(l > r) swap(l,r);
ans = val[query(l,r)];
printf("%d\n",ans);
}
return 0;
}
【BZOJ 2821】作诗
给一个序列,求一个区间内出现正偶次的数的个数。(强制在线)
【题解】
本题和上面一题的套路是一样的。
用f[i][j]表示第i块到第j块中符合条件的数的个数。
对于一个询问l~r,中间大块x~y已经知晓,要做的就是看l-x和y-r之间的数对答案的影响。
用vector保存每个数的位置,分别统计除了大块外的那些数(不超过2*sqrt(n)个)在x~y中出现的次数t1,以及在l~r之间出现的次数t2,根据t1和t2的奇偶性判断对答案的影响。同样,统计t1、t2也需要在vector中二分。
注:对f[i][j]的预处理可以在O(n * sqrt(n))的时间完成。
同时,块的大小取sqrt(n/log n)时最优。
【代码】
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
#define N 100010
#define B 1505
using namespace std;
vector<int> pos[N];
int a[N],belong[N],L[B],R[B];
int n,bs;
int f[B+5][B+5],tmp[N];
bool vis[N];
void pre() {
int tot;
for(int i = 1;i <= bs;i++) {
for(int j = L[i];j <= n;j++) tmp[a[j]] = 0;
tot = 0;
for(int j = L[i];j <= n;j++) {
if(!(tmp[a[j]]&1) && tmp[a[j]] != 0) tot--;
tmp[a[j]]++;
if(!(tmp[a[j]]&1) )tot++;
f[i][belong[j]] = tot;
}
}
}
int sigma(int l,int r,int x) {
return max(upper_bound(pos[x].begin(),pos[x].end(),r)-lower_bound(pos[x].begin(),pos[x].end(),l),0);
}
int query(int l,int r)
{
int x = belong[l],y = belong[r],ans = 0,t1,t2;
memset(vis,0,sizeof(vis));
if(x == y || x+1 == y) {
for(int i = l;i <= r;i++)
if(!vis[a[i]]) {
if(!(sigma(l,r,a[i])&1)) ans++;
vis[a[i]] = 1;
}
}
else {
x++; y--;
ans = f[x][y];
for(int i = l;i <= L[x]-1;i++)
if(!vis[a[i]]) {
t2 = sigma(l,L[x]-1,a[i])+sigma(R[y]+1,r,a[i]);
t1 = sigma(L[x],R[y],a[i]);
if(!(t2&1) && t2 != 0 && t1 == 0) ans++;
if((t2&1) && (t1&1)) ans++;
if((t2&1) && !(t1&1) && t1 != 0) ans--;
vis[a[i]] = 1;
}
for(int i = R[y]+1;i <= r;i++)
if(!vis[a[i]]) {
t2 = sigma(l,L[x]-1,a[i])+sigma(R[y]+1,r,a[i]);
t1 = sigma(L[x],R[y],a[i]);
if(!(t2&1) && t2 != 0 && t1 == 0) ans++;
if((t2&1) && (t1&1)) ans++;
if((t2&1) && !(t1&1) && t1 != 0) ans--;
vis[a[i]] = 1;
}
}
return ans;
}
int main()
{
int c,m,l,r;
scanf("%d%d%d",&n,&c,&m);
for(int i = 1;i <= n;i++) {
scanf("%d",&a[i]);
pos[a[i]].push_back(i);
}
int block=sqrt((double)n/log((double)n)*log(2));
if(n % block) bs = n/block+1; else bs = n/block;
for(int i = 1;i <= n;i++) belong[i] = (i-1)/block+1;
for(int i = 1;i <= bs;i++) {
L[i] = (i-1)*block+1; R[i] = i*block;
}
R[bs] = n;
pre();
int ans = 0;
while(m--) {
scanf("%d%d",&l,&r);
l = (l+ans) % n + 1;r = (r+ans) % n + 1;
if(l > r) swap(l,r);
ans = query(l,r);
printf("%d\n",ans);
}
return 0;
}