[bzoj 4939][Ynoi 2016]掉进兔子洞
Posted papercloud
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[bzoj 4939][Ynoi 2016]掉进兔子洞相关的知识,希望对你有一定的参考价值。
传送门
Description
一个长为 n 的序列 a。
有 m 个询问,每次询问三个区间,把三个区间中同时出现的数一个一个删掉,问最后三个区间剩下的数的个数和,询问独立。
注意这里删掉指的是一个一个删,不是把等于这个值的数直接删完,
比如三个区间是 [1,2,2,3,3,3,3] , [1,2,2,3,3,3,3] 与 [1,1,2,3,3],就一起扔掉了 1 个 1,1 个 2,2 个 3。
Solution
弱化版是luoguP3674,这里放上它的题解 戳我
要找相同的元素就是集合求交集,莫队+bitset即可
可是这道题是要算上重复的元素的,很简单,离散化的时候给每个数留下相应的空位就行了
即每个数的离散值等于小于它的数个数+1.
几个例子:
我们把序列
3 3 5 5 6 7 7
,离散化成1 1 3 3 5 6 6
,而不是1 1 2 2 3 4 4
那么,我们在加入数的时候,把它加到“它的离散值+它出现的次数-1”这个位置上
那么求交就是
now.count()
还有一个问题,维护(100000*100000)的bitset!虽然看上去不太可能,但其实上我们分成3次莫队就行了
Code?
#include<bits/stdc++.h>
#define min(a,b) ((a)<(b)?(a):(b))
inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*f;
}
#define MN 100005
int n,m,a[MN],T,nums[MN];
std::bitset<MN> now;
std::bitset<MN> quer[33343];
struct ques{
int l,r,id,pl;
bool operator <(const ques&o)const{return (pl^o.pl)?(pl<o.pl):((pl&1)?r<o.r:r>o.r);}
}q[MN+33343];
int num[MN],len[MN];
inline void solve(int L,int R)
{
if(L>R) return;
register int tt=0,i;
for(i=1;i<=R-L+1;++i)
{
quer[i].set();
q[++tt].l=read(),q[tt].r=read();q[tt].pl=(q[tt].l-1)/T+1;q[tt].id=i;len[i] =q[tt].r-q[tt].l+1;
q[++tt].l=read(),q[tt].r=read();q[tt].pl=(q[tt].l-1)/T+1;q[tt].id=i;len[i]+=q[tt].r-q[tt].l+1;
q[++tt].l=read(),q[tt].r=read();q[tt].pl=(q[tt].l-1)/T+1;q[tt].id=i;len[i]+=q[tt].r-q[tt].l+1;
}
std::sort(q+1,q+tt+1);
register int l=1,r=0;now.reset();
memset(num,0,sizeof num);
for(i=1;i<=tt;++i)
{
for(;r<q[i].r;++r) num[a[r+1]]++,now[a[r+1]+num[a[r+1]]-1]=1;
for(;l>q[i].l;--l) num[a[l-1]]++,now[a[l-1]+num[a[l-1]]-1]=1;
for(;r>q[i].r;--r) now[a[r]+num[a[r]]-1]=0,--num[a[r]];
for(;l<q[i].l;++l) now[a[l]+num[a[l]]-1]=0,--num[a[l]];
quer[q[i].id]&=now;
}
for(i=1;i<=R-L+1;++i) printf("%d
",len[i]-3*quer[i].count());
}
int main()
{
n=read();m=read();
register int i;T=ceil(sqrt((double)n));
for(i=1;i<=n;++i) nums[i]=a[i]=read();
std::sort(nums+1,nums+n+1);
for(i=1;i<=n;++i) a[i]=std::lower_bound(nums+1,nums+n+1,a[i])-nums;
int N1=m/3+1,N2=2*N1;
solve(1,min(N1,m));solve(min(N1,m)+1,min(N2,m));solve(min(N2,m)+1,m);
return 0;
}
Blog来自PaperCloud,未经允许,请勿转载,TKS!
以上是关于[bzoj 4939][Ynoi 2016]掉进兔子洞的主要内容,如果未能解决你的问题,请参考以下文章