[2018雅礼集训1-20]B 分块

题面
考虑把序列分成 N−−√ 块，记 cnti,j 表示第 i 块颜色 j 的个数，设 si,j 为 i,j 两块之间产生的贡献，通过 cnt 可以 O(N−−√) 求出，然后对其二维前缀和一下，就可以 O(1) 求出一段块的贡献了。
对于单独的块再考虑其内部，和其对中间整块的贡献，加上去即可。
复杂度 O(NN−−√)
代码：

#include
#include
#include
#include
#define L(x) (ID(x)*B)
#define R(x) min((ll)n-1,L(x)+B-1)
#define ID(x) (x/B)
#define N 50010
#define ll long long
using namespace std;
const int B=300;
int n,m,type,a[N],cnt[200][N],cct[N];
ll s[200][200];
ll read()
{
    ll x=0;char ch=getchar();
    for(;ch<'0'||ch>'9';ch=getchar());
    for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0';
    return x;
}
ll cal(int l,int r,int L,int R)
{
    ll re=0;

    for(int i=l;i<=r;i++)
    {
        if(!cct[a[i]]) cct[a[i]]=cnt[R][a[i]]-cnt[L-1][a[i]];
        re+=cct[a[i]],cct[a[i]]++;
    }
    return re;      
}
void clr(int l,int r)
{
    for(int i=l;i<=r;i++)
        cct[a[i]]=0;    
}
int main()
{
    n=read();m=read();type=read();
    for(int i=0;iread(),cnt[ID(i)][a[i]]++;
    for(int i=0;i*Bfor(int j=0;j<=i;j++)
            if(i==j)
                for(int k=1;k<=50000;k++)
                    s[i][j]+=cnt[i][k]*(cnt[i][k]-1)/2;
            else
            {
                for(int k=i*B;k*B+B,n);k++)
                    if(!cct[a[k]]) cct[a[k]]=1,s[i][j]+=cnt[i][a[k]]*cnt[j][a[k]];  
                clr(i*B,min(i*B+B,n)-1);
            }
    for(int i=1;i*Bfor(int j=1;j<=50000;j++)
            cnt[i][j]+=cnt[i-1][j]; 
    for(int i=0;i*Bfor(int j=0;j*Bs[i][j]+=(i?s[i-1][j]:0)+(j?s[i][j-1]:0)-(i&&j?s[i-1][j-1]:0);      
    ll ans=0;
    while(m--)  
    {
        ll l=read(),r=read();
        if(type) l^=ans,r^=ans;
        l=l%n;r=r%n;
        if(l>r) swap(l,r);
        ans=0;
        if(ID(l)+1<=ID(r)-1) ans+=s[ID(r)-1][ID(r)-1]-s[ID(r)-1][ID(l)]-s[ID(l)][ID(r)-1]+s[ID(l)][ID(l)];
        if(ID(l)==ID(r)) ans+=cal(l,r,1,0),clr(l,r);

        else ans+=cal(l,R(l),ID(l)+1,ID(r)-1)+cal(L(r),r,ID(l)+1,ID(r)-1),clr(l,R(l)),clr(L(r),r);
        printf("%lld\n",ans);
    }
    return 0;
}