洛谷P4135作诗（静态分块_区间贡献计算洛谷-CSDN博客

本文链接：https://blog.csdn.net/weixin_45539557/article/details/116465673

该博客探讨了一种在线算法问题，即查询区间内出现偶数次的数的个数。文章介绍了如何利用静态分块策略处理此问题，详细解释了如何预处理关键点信息以加速查询，并提供了相应的代码实现。静态分块方法在不支持修改操作的情况下，作为莫队算法的变种，用于解决强制在线场景下的区间查询问题。博客还讨论了离线和在线情况下的算法选择及其优势。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

洛谷P4135作诗（静态分块

题意：

查询区间内出现偶数次的数的个数，强制在线，值域 $O (n)$

思路：

假如不强制在线的莫队板子题，但是强制在线怎么办呢，此时我们考虑静态分块，即不带修，放一些关键点，预处理关键点到关键点的信息来加速查询。显然块之间的贡献不是独立的，对于这种问题，我们一般直接预处理所有对应块的答案，然后考虑散块、整块之间的贡献。我们考虑 $f [l] [r]$ 表示从块l到块r的答案个数。可以直接通过增量法计算（枚举最左端的块l，再一直扫到n），复杂度 $O(n^2/B)$

然后是散块的贡献，我们用 $c n t [i] [j]$ 表示前 $i$ 个块中 $j$ 的出现次数，可以 $O (n B)$ 预处理出来，计算散块的时候，我们暴力枚举散块部分对整块的贡献考虑是否更新答案即可

这里我们可以发现，静态分块，即不支持修改操作，放一些关键点，预处理关键点到关键点的信息来加速查询，假如可以离线，功能莫队的子集，所以强制在线的时候分块才有优越性，具体问题如区间逆序对，区间众数

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+5;
int cnt[400][maxn],ans[400][400],a[maxn],big,num,L[400],R[400],pos[maxn],sum[maxn],n,c,m;
void init(){
    big=sqrt(n);
    num=(n-1)/big+1;
    for(int i=1;i<=num;++i){
        L[i]=(i-1)*big+1;
        R[i]=i*big;
    }
    R[num]=n;
    for(int i=1;i<=num;++i){
        for(int j=L[i];j<=R[i];++j)
            pos[j]=i;
    }
    for(int i=1;i<=n;++i){
        cnt[pos[i]][a[i]]++;
    }
    for(int i=2;i<=num;++i){
        for(int j=1;j<=c;++j)
            cnt[i][j]+=cnt[i-1][j];
    }
    for(int i=1;i<=num;++i){
        int s=0;
        for(int j=L[i];j<=n;++j){
            sum[a[j]]++;
            if((sum[a[j]]&1)==0)s++;//维护到当前为止的个数！
            else if(sum[a[j]]>2)s--;//&的优先级比==低
            ans[i][pos[j]]=s;
        }
        for(int j=L[i];j<=n;++j)sum[a[j]]=0;       
    }
}
int query(int l,int r){
    int pl=pos[l],pr=pos[r],res=0;
    if(pl==pr){
        for(int i=l;i<=r;++i){
            sum[a[i]]++;
            if((sum[a[i]]&1)==0)res++;
            else if(sum[a[i]]>2)res--;
        }
        for(int i=l;i<=r;++i)sum[a[i]]--;
    }else{
        res=ans[pl+1][pr-1];
        for(int i=l;i<=R[pl];++i)sum[a[i]]=cnt[pr-1][a[i]]-cnt[pl][a[i]];
        for(int i=l;i<=R[pl];++i){
            sum[a[i]]++;
            if((sum[a[i]]&1)==0)res++;
            else if(sum[a[i]]>2)res--;
        }
        for(int i=L[pr];i<=r;++i)
            sum[a[i]]=sum[a[i]]?sum[a[i]]:cnt[pr-1][a[i]]-cnt[pl][a[i]];
        for(int i=L[pr];i<=r;++i){
            sum[a[i]]++;
            if((sum[a[i]]&1)==0)res++;
            else if(sum[a[i]]>2)res--;
        }
        for(int i=l;i<=R[pl];++i)sum[a[i]]=0;
        for(int i=L[pr];i<=r;++i)sum[a[i]]=0;
    }
    return res;
}
int main(){
    scanf("%d%d%d",&n,&c,&m);
    for(int i=1;i<=n;++i)scanf("%d",&a[i]);
    init();
    int l,r,anss=0;
    for(int i=1;i<=m;++i){
        scanf("%d%d",&l,&r);
        l=(l+anss)%n+1;r=(r+anss)%n+1;
        if(l>r)swap(l,r);
        cout<<(anss=query(l,r))<<"\n";
    }
    return 0;
}