【bzoj4939】【YNOI2016】掉进兔子洞（莫队）

题目描述

您正在打galgame，然后突然发现您今天太颓了，于是想写个数据结构题练练手： 一个长为 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。

输入输出格式

输入格式：

第一行两个数表示 n , m。

第二行 n个数表示 a[i]​。

之后 m 行，每行 6 个数 l1​​ , r1​​ , l2, r2​​ , l3​​ , r3​​ 表示这三个区间。

输出格式：

对于每个询问，输出一个数表示答案。

输入输出样例

输入样例#1： 复制

5 2
1 2 2 3 3
1 2 2 3 3 4
1 5 1 5 1 5

输出样例#1： 复制

3
0

说明

n , m <= 100000 , 1 <= a[i] <= 1000000000

题解

　　大毒瘤题名不虚传……

　　做莫队的时候因为先del再add已经快RE到死了……

　　据某加藤大佬说这题一看就是莫队+bitset维护并集，然而我啥都看不出来……

　　先把原数组给离散，然后总共只有$10^5$个数，可以对每一个询问维护一个bitset，表示每一位有几个，然后只要把三个询问的bitset并起来就行了。然而bitset怎么表示数的个数呢？我们可以给每一个数很多位置，位置数为它在原数组中的个数。比方说1,1,4,3,1，那么我们就给1这个数字3个位置，做莫队的时候，维护一个bitset，设当前已经有x个1，且1在bitset中的位置为1，那么要add的时候，就让第x+1位变为1，表示加了一个1，del的话，就让第x位变为0，表示少了一个1，同时更新x

　　然后这样有可能两个数的区间重叠，那么我们在离散化的之后可以不去重，直接用在排序后的数组的位置表示新数，这样每一个数就不会和其他区间重叠了

 //minamoto
 #include<iostream>
 #include<cstdio>
 #include<bitset>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 using namespace std;
 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
 char buf[<<],*p1=buf,*p2=buf;
 inline int read(){
     #define num ch-'0'
     char ch;bool flag=;int res;
     while(!isdigit(ch=getc()))
     (ch=='-')&&(flag=true);
     for(res=num;isdigit(ch=getc());res=res*+num);
     (flag)&&(res=-res);
     #undef num
     return res;
 }
 char sr[<<],z[];int C=-,Z;
 inline void Ot(){fwrite(sr,,C+,stdout),C=-;}
 inline void print(int x){
     if(C><<)Ot();if(x<)sr[++C]=,x=-x;
     while(z[++Z]=x%+,x/=);
     while(sr[++C]=z[Z],--Z);sr[++C]='\n';
 }
 const int N=1e5+,NN=;
 int n,m,len,rt[N],ans[N];
 bitset<N> f[NN+],tmp;
 int a[N],b[N],L1[N],R1[N],L2[N],R2[N],L3[N],R3[N],cnt[N],l,r,tot;
 bool flag[NN+];
 struct node{
     int l,r,id;
     node(){}
     node(int l,int r,int id):l(l),r(r),id(id){}
 }q[N];
 inline bool operator <(node a,node b){
     return rt[a.l]==rt[b.l]?rt[a.l]&?a.r<b.r:a.r>b.r:rt[a.l]<rt[b.l];
 }
 inline void init(){
     memset(cnt,,sizeof(cnt));
     memset(flag,,sizeof(flag));
     tmp.reset();
     l=,r=,tot=;
 }
 inline void add(int x){
     tmp[x+cnt[x]]=,++cnt[x];
 }
 inline void del(int x){
     tmp[x+cnt[x]-]=,--cnt[x];
 }
 inline void solve(int lx,int rx){
     init();
     for(int i=lx;i<=rx;++i){
         q[++tot]=node(L1[i],R1[i],i),ans[i]+=R1[i]-L1[i]+;
         q[++tot]=node(L2[i],R2[i],i),ans[i]+=R2[i]-L2[i]+;
         q[++tot]=node(L3[i],R3[i],i),ans[i]+=R3[i]-L3[i]+;
     }
     sort(q+,q++tot);
     for(int i=;i<=tot;++i){
         while(l>q[i].l) add(a[--l]);
         while(r<q[i].r) add(a[++r]);
         while(l<q[i].l) del(a[l++]);
         while(r>q[i].r) del(a[r--]);
         if(!flag[q[i].id-lx+]) flag[q[i].id-lx+]=,f[q[i].id-lx+]=tmp;
         else f[q[i].id-lx+]&=tmp;
     }
     for(int i=lx;i<=rx;++i) ans[i]-=f[i-lx+].count()*;
 }
 int main(){
     n=read(),m=read(),len=sqrt(n);
     for(int i=;i<=n;++i) a[i]=b[i]=read(),rt[i]=(i-)/len+;
     sort(b+,b++n);
     for(int i=;i<=n;++i) a[i]=lower_bound(b+,b++n,a[i])-b;
     for(int i=;i<=m;++i)
     L1[i]=read(),R1[i]=read(),L2[i]=read(),R2[i]=read(),L3[i]=read(),R3[i]=read();
     for(int i=;i<=m;i+=NN) solve(i,min(m,i+NN-));
     for(int i=;i<=m;++i) print(ans[i]);
     Ot();
     return ;
 }