[bzoj 5143][Ynoi 2018]五彩斑斓的世界

Posted papercloud

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[bzoj 5143][Ynoi 2018]五彩斑斓的世界相关的知识,希望对你有一定的参考价值。

传送门

Descroption

给了你一个长为n的序列a,有m次操作

1.把区间[l,r]中大于x的数减去x

2.查询区间[l,r]中x的出现次数

Solution

分块

对于每个块,我们都分别维护:

  • 每种数出现的次数,需要用到并查集,即把相同的数合并在一起

  • 标记tag,表示这个块总共减了多少
  • maxx,表示这个块的最大值

对于查询,按照分块的套路直接查询就可以了

修改的时候,为了保证复杂度正确,我们考虑每次通过$O ( x ) (的复杂度,实现整个块的极差减少)x$.

这样,我们就能保证总修改的复杂度是(O(nlog n))

怎么做到呢?

假设当前数组的极差为(M),我们近似的将它看成是(maxx-tag)

  • (M< x) 显然是不存在大于x的数,直接返回即可
  • (x leq M leq 2x) 这个时候,把大于x的数直接减去x即可,fa[j]=fa[j-x],这样极差减少(x)
  • (M>2x) 直接把小于等于(x)的数加上(x),然后整个区间打上减(x)的标记,极差同样减少(x)


Code?

#include<bits/stdc++.h>
#define ll long long
using namespace std;

char B[1<<26],*S=B;
inline int read(){
    int x;char c;
    while((c=*S++)<'0'||c>'9');
    for(x=c-'0';(c=*S++)>='0'&&c<='9';)x=x*10+c-'0';
    return x;
}
const int MN=1e5,T=650;
int a[MN+5],fa[MN/T+5][MN+5],sz[MN/T+5][MN+5],minn[T+5],maxx[T+5];
int getf(int*f,int k){return f[k]?f[k]=getf(f,f[k]):k;}
int main()
{
    B[fread(B,1,1<<26,stdin)]=0;
    int n,m,i,j,o,l,r,x,lk,rk,ans;
    n=read();m=read();
    for(i=1;i<=n;++i)a[i]=read(),++sz[(i-1)/T][a[i]];
    for(i=1;i<=n;i+=T) minn[i/T]=0,maxx[i/T]=MN;
    while(m--)
    {
        o=read();l=read();r=read();x=read();
        lk=(l-1)/T;rk=(r-1)/T;
        if(o==1)
        {
            for(i=l;i<=r&&i<=lk*T+T;++i)
                if((a[i]=getf(fa[lk],a[i]))-minn[lk]>x)--sz[lk][a[i]],++sz[lk][a[i]-=x];
            for(i=lk;++i<rk;)
                if(maxx[i]-minn[i]>2*x)
                {
                    for(j=1;j<=x;++j)sz[i][fa[i][minn[i]+j]=minn[i]+j+x]+=sz[i][minn[i]+j];
                    minn[i]+=x;
                }
                else
                {
                    for(j=x+1;j<=maxx[i]-minn[i];++j)sz[i][fa[i][minn[i]+j]=minn[i]+j-x]+=sz[i][minn[i]+j];
                    maxx[i]=min(maxx[i],minn[i]+x); 
                }
            if(lk!=rk)for(i=r;i>rk*T;--i)
                if((a[i]=getf(fa[rk],a[i]))-minn[rk]>x)--sz[rk][a[i]],++sz[rk][a[i]-=x];
        }
        else
        {
            ans=0;
            for(i=l;i<=r&&i<=lk*T+T;++i)if(getf(fa[lk],a[i])-minn[lk]==x)++ans;
            for(i=lk;++i<rk;)if(x+minn[i]<=maxx[i])ans+=sz[i][x+minn[i]];
            if(lk!=rk)for(i=r;i>rk*T;--i)if(getf(fa[rk],a[i])-minn[rk]==x)++ans;
            printf("%d
",ans);
        }
    }
}



Blog来自PaperCloud,未经允许,请勿转载,TKS!

以上是关于[bzoj 5143][Ynoi 2018]五彩斑斓的世界的主要内容,如果未能解决你的问题,请参考以下文章

[Ynoi2018]五彩斑斓的世界

[Ynoi2018]五彩斑斓的世界

[Ynoi2018]未来日记

BZOJ4810: [Ynoi2017]由乃的玉米田

BZOJ4811 [Ynoi2017]由乃的OJ

BZOJ4810 [Ynoi2017]由乃的玉米田