洛谷P4135作诗(静态分块

洛谷P4135作诗(静态分块

题意:

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

思路:

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

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

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

#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;
}
  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值