【bzoj 2821】 作诗 题意&题解&代码（C++）

题目链接：

http://www.lydsy.com/JudgeOnline/problem.php?id=2821

权限题，附题意：

题解：

分块题，他要求在线处理区间问题，很明显是为了卡莫队这种神奇算法，首先将序列分块，按 sqrt(n) （即根号n）分块，然后预处理出一个数组 f[i][j] 表示从第 i 块到第 j 块内，满足题意的答案，还要处理出另一个数组 cnt[i][j] 表示前 i 块中数字 j 出现了多少次，（不要担心数组太大开不下，三千两百万它还真就是存下了，o.o！），这样的话当询问区间 （i，j）时，我们利用 f[i][j] 通过各种判断首先将包含在这个区间里的块的答案统计出来，再利用 cnt[i][j] 将两头或是一头没有达到一整快的部分暴力的统计，一边统计一边修改已经得到的答案即可，时限还是够的。

注意题中一个特别坑的地方，

他特别心机的的在mod与n+1 之间加了一个空格，让大家以为是mod（n+1）实际上却是（mod n）+1。。2333

---

代码

#include<iostream>#include<algorithm>#include<stdio.h>#include<string.h>#include<math.h>#define maxn (100005)#define maxs (320)using namespace std;int vis[maxn],ans,L,R,len,n,num,c,m,sq[maxn],f[maxs][maxs],cnt[maxs][maxn],a[maxn],st[maxs],tmp[maxn];void solve(int l,int r)//计算 l,r{    if (sq[l]==sq[r] || sq[l]+1==sq[r])    {        for (int i=l;i<=r;i++)        {            vis[a[i]]++;            if (vis[a[i]]%2==0) ans++;            else if (vis[a[i]]!=1) ans--;        }        for (int i=l;i<=r;i++)        vis[a[i]]=0;        return ;    }    if (l==st[sq[l]] && r==st[sq[r]+1]-1)    {        ans=f[sq[l]][sq[r]];        return ;        }    else if (l==st[sq[l]])    {        ans=f[sq[l]][sq[r]-1];        for (int i=st[sq[r]];i<=r;i++)        vis[a[i]]=cnt[sq[r]-1][a[i]]-cnt[sq[l]-1][a[i]];        for (int i=st[sq[r]];i<=r;i++)        {            vis[a[i]]++;            if (vis[a[i]]%2==0) ans++;            else if (vis[a[i]]!=1) ans--;               }                for (int i=st[sq[r]];i<=r;i++)                vis[a[i]]=0;        return ;                    }    else if (r==st[sq[r]+1]-1)    {        ans=f[sq[l]+1][sq[r]];        for (int i=l;i<st[sq[l]+1];i++)        vis[a[i]]=cnt[sq[r]][a[i]]-cnt[sq[l]][a[i]];        for (int i=l;i<st[sq[l]+1];i++)        {            vis[a[i]]++;            if (vis[a[i]]%2==0) ans++;            else if (vis[a[i]]!=1) ans--;           }                for (int i=l;i<st[sq[l]+1];i++)                vis[a[i]]=0;        return ;    }    else    {        int x=sq[l]+1;        int y=sq[r]-1;        ans=f[x][y];        for (int i=l;i<st[x];i++)        vis[a[i]]=cnt[y][a[i]]-cnt[x-1][a[i]];        for (int i=st[y+1];i<=r;i++)        vis[a[i]]=cnt[y][a[i]]-cnt[x-1][a[i]];        for (int i=l;i<st[x];i++)        {            vis[a[i]]++;            if (vis[a[i]]%2==0) ans++;            else if (vis[a[i]]!=1) ans--;           }        for (int i=st[y+1];i<=r;i++)        {            vis[a[i]]++;            if (vis[a[i]]%2==0) ans++;            else if (vis[a[i]]!=1) ans--;               }           for (int i=l;i<st[x];i++) vis[a[i]]=0;        for (int i=st[y+1];i<=r;i++) vis[a[i]]=0;           return ;    }}int main(){    scanf("%d%d%d",&n,&c,&m);    len=sqrt(n);    for (int i=1;i<=n;i++)    {        scanf("%d",&a[i]);        sq[i]=(i-1)/len+1;        cnt[sq[i]][a[i]]++;        if (sq[i]!=sq[i-1])        st[sq[i]]=i;    }    num=sq[n];    for (int i=1;i<=num;i++)    for (int j=1;j<=c;j++)    cnt[i][j]+=cnt[i-1][j];    //预处理 cnt数组    for (int i=1;i<=num;i++)    {        int tans=0;        for (int j=1;j<=c;j++)        tmp[j]=0;        for (int j=st[i];j<=n;j++)        {            tmp[a[j]]++;            if (tmp[a[j]]%2==0) tans++;            else if (tmp[a[j]]!=1) tans--;            f[i][sq[j]]=tans;        }    }//预处理 f 数组    ans=0;    for (int i=1;i<=m;i++)    {        scanf("%d%d",&L,&R);        L=(L+ans)%n+1;        R=(R+ans)%n+1;        ans=0;        if (L>R) swap(L,R);        solve(L,R);        printf("%d\n",ans);    }}

不得不说，csdn的博客编辑最近好像十分的恶心。。。