[BZOJ2821]作诗(Poetize)

[BZOJ2821]作诗(Poetize)

题目大意:

一个长度为\(n(n\le10^5)\)的数列,\(m(m\le10^5)\)次询问,每次询问区间\([l,r]\)内出现次数为正偶数的数的个数。强制在线。

思路:

分块,\(f[i][j]\)表示第\(i\)块到第\(j\)块的答案,\(s[i][j]\)表示\(j\)在前\(i\)块出现次数。

时间复杂度\(\mathcal O(n\sqrt n)\)

源代码:

#include<cstdio>
#include<cctype>
#include<algorithm>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
const int N=1e5+1,B=317;
int n,c,m,a[N],bel[N],f[B][B],s[B][N],cnt[N];
inline void init() {
    for(register int i=1;i<=n;i++) {
        bel[i]=i/B+1;
        s[bel[i]][a[i]]++;
    }
    for(register int i=1;i<=bel[n];i++) {
        for(register int j=1;j<=c;j++) {
            s[i][j]+=s[i-1][j];
        }
    }
    for(register int i=1;i<=bel[n];i++) {
        for(register int j=i;j<=bel[n];j++) {
            f[i][j]=f[i][j-1];
            for(register int k=std::max(B*(j-1),1);k<B*j&&k<=n;k++) {
                cnt[a[k]]++;
                f[i][j]+=cnt[a[k]]>0&&cnt[a[k]]%2==0;
                f[i][j]-=cnt[a[k]]>1&&cnt[a[k]]%2==1;
            }
        }
        for(register int k=std::max(B*(i-1),1);k<=n;k++) {
            cnt[a[k]]--;
        }
    }
}
inline int query(const int &l,const int &r) {
    if(bel[l]==bel[r]) {
        int ret=0;
        for(register int i=l;i<=r;i++) {
            cnt[a[i]]++;
            ret+=cnt[a[i]]>0&&cnt[a[i]]%2==0;
            ret-=cnt[a[i]]>1&&cnt[a[i]]%2==1;
        }
        for(register int i=l;i<=r;i++) {
            cnt[a[i]]--;
        }
        return ret;
    }
    int ret=f[bel[l]+1][bel[r]-1];
    for(register int i=l;bel[i]==bel[l];i++) {
        cnt[a[i]]=s[bel[r]-1][a[i]]-s[bel[l]][a[i]];
    }
    for(register int i=r;bel[i]==bel[r];i--) {
        cnt[a[i]]=s[bel[r]-1][a[i]]-s[bel[l]][a[i]];
    }
    for(register int i=l;bel[i]==bel[l];i++) {
        cnt[a[i]]++;
        ret+=cnt[a[i]]>0&&cnt[a[i]]%2==0;
        ret-=cnt[a[i]]>1&&cnt[a[i]]%2==1;
    }
    for(register int i=r;bel[i]==bel[r];i--) {
        cnt[a[i]]++;
        ret+=cnt[a[i]]>0&&cnt[a[i]]%2==0;
        ret-=cnt[a[i]]>1&&cnt[a[i]]%2==1;
    }
    for(register int i=l;bel[i]==bel[l];i++) {
        cnt[a[i]]=0;
    }
    for(register int i=r;bel[i]==bel[r];i--) {
        cnt[a[i]]=0;
    }
    return ret;
}
int main() {
    n=getint(),c=getint(),m=getint();
    for(register int i=1;i<=n;i++) {
        a[i]=getint();
    }
    init();
    for(register int i=0,ans=0;i<m;i++) {
        int l=(getint()+ans)%n+1,r=(getint()+ans)%n+1;
        if(l>r) std::swap(l,r);
        printf("%d\n",ans=query(l,r));
    }
    return 0;
}
posted @ 2019-01-07 16:31 skylee03 阅读(...) 评论(...) 编辑 收藏