【bzoj2821】【作诗】【分块】

Description

神犇SJY虐完HEOI之后给傻×LYD出了一题:SHY是T国的公主,平时的一大爱好是作诗。由于时间紧迫,SHY作完诗
之后还要虐OI,于是SHY找来一篇长度为N的文章,阅读M次,每次只阅读其中连续的一段[l,r],从这一段中选出一
些汉字构成诗。因为SHY喜欢对偶,所以SHY规定最后选出的每个汉字都必须在[l,r]里出现了正偶数次。而且SHY认
为选出的汉字的种类数(两个一样的汉字称为同一种)越多越好(为了拿到更多的素材!)。于是SHY请LYD安排选
法。LYD这种傻×当然不会了,于是向你请教……问题简述:N个数,M组询问,每次问[l,r]中有多少个数出现正偶
数次。

Input

输入第一行三个整数n、c以及m。表示文章字数、汉字的种类数、要选择M次。第二行有n个整数,每个数Ai在[1, c
]间,代表一个编码为Ai的汉字。接下来m行每行两个整数l和r,设上一个询问的答案为ans(第一个询问时ans=0),
令L=(l+ans)mod n+1, R=(r+ans)mod n+1,若L>R,交换L和R,则本次询问为[L,R]。

Output

输出共m行,每行一个整数,第i个数表示SHY第i次能选出的汉字的最多种类数。

Sample Input

5 3 5
1 2 2 3 1
0 4
1 2
2 2
2 3
3 5

Sample Output

2
0
0
0
1

HINT

对于100%的数据,1<=n,c,m<=10^5

题解:

           对序列分块。

           用t[i][j]表示第i块到第j块出现次数为偶数的数的个数.

           对于每次询问先把整块的答案加上.

           然后暴力枚举剩下的部分,可以算出一个数在整块之内和块之外的出现次数.

           分情况加减一下即可.

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define N 100010
#define M 2010
using namespace std;
int a[N],bk,cnt,t[M][M],l[N],r[N],n,m,x,y,c,ans,bl[N],sum[N];
struct block{int l,r;}st[M];
struct use{int p,x;}q[N];
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
bool cmp(use a,use b){
  if (a.x==b.x) return a.p<b.p;
  else return a.x<b.x;
}
void getbk(){
  bk=sqrt((double)n/log((double)n)*log(2));
  if (n%bk) cnt=n/bk+1;else cnt=n/bk;
  for (int i=1;i<=n;i++) bl[i]=(i-1)/bk+1;
  for (int i=1;i<=cnt;i++)
    st[i].l=(i-1)*bk+1,st[i].r=i*bk;
}
void pre(){
  for (int i=1;i<=n;i++) q[i].x=a[i],q[i].p=i;
  sort(q+1,q+n+1,cmp);
  for (int i=1;i<=n;i++){
    if (!l[q[i].x]) l[q[i].x]=i;
    r[q[i].x]=i;
  }
  for (int i=1;i<=cnt;i++){
    int ans(0);memset(sum,0,sizeof(sum));
    for (int j=st[i].l;j<=n;j++){
      sum[a[j]]++;
      if (sum[a[j]]%2==0) ans++;
      else if (sum[a[j]]>1) ans--;      
      t[i][bl[j]]=ans;
    } 
  }
}
int up(int k,int x,int y){
   int l=x,r=y,ans=0;
   while(l<=r){
     int mid=(l+r)>>1;
     if (q[mid].p>k) r=mid-1;
     else {ans=mid;l=mid+1;}
   }
   return ans;
}
int down(int k,int x,int y){
  int l=x,r=y,ans=N;
  while (l<=r){
    int mid=(l+r)>>1;
    if (q[mid].p<k) l=mid+1;
    else {ans=mid;r=mid-1;}
  }
  return ans;
}
int cal(int x,int y,int v){
  return max(up(y,l[v],r[v])-down(x,l[v],r[v])+1,0);
}
int query(int x,int y){
  int ans(0);
  if (bl[x]==bl[y]||bl[x]+1==bl[y]){
    for (int i=x;i<=y;i++){
       if (sum[a[i]]) continue;sum[a[i]]++;
       int temp=cal(x,y,a[i]);
       if (temp%2==0) ans++;
    }
    for (int i=x;i<=y;i++) sum[a[i]]=0; 
  }
  else{
    ans+=t[bl[x]+1][bl[y]-1]; 
    for (int i=x;i<=st[bl[x]].r;i++){
      if (sum[a[i]]) continue;sum[a[i]]++;
      int t1=cal(x,y,a[i])&1;
      int t2=cal(st[bl[x]+1].l,st[bl[y]-1].r,a[i]);  
      if (t2==0&&t1==0) ans++;
      if (t2%2==1&&t1==0) ans++;
      if (t2!=0&&t2%2==0&&t1==1) ans--;
    }
    for (int i=st[bl[y]].l;i<=y;i++){
      if (sum[a[i]]) continue;sum[a[i]]++;
      int t1=cal(x,y,a[i])&1;
      int t2=cal(st[bl[x]+1].l,st[bl[y]-1].r,a[i]);  
      if (t2==0&&t1==0) ans++;
      if ((t2&1)&&t1==0) ans++;
      if (t2!=0&&t2%2==0&&t1==1) ans--;
    }
    for (int i=x;i<=st[bl[x]].r;i++) sum[a[i]]=0;
    for (int i=st[bl[y]].l;i<=y;i++) sum[a[i]]=0;
  }
  return ans;
}
int main(){
  n=read();c=read();m=read();
  for (int i=1;i<=n;i++) a[i]=read(); 
  getbk();pre();memset(sum,0,sizeof(sum));
  for (int i=1;i<=m;i++){
    x=read();y=read();
    x=(x+ans)%n+1;y=(y+ans)%n+1;
    if (x>y) swap(x,y);
    ans=query(x,y);
    printf("%d\n",ans);
  }
} 



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值