bzoj 4939: [Ynoi2016]掉进兔子洞

有 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。

 

Input

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

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

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

 

Output

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

 

Sample Input

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

Sample Output

3
0

HINT

 

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

 

首先，很显然这题是要用莫队来处理的。我们先把输入的数字另外排一下序，然后记录一下pi表示每一个数字对应在排好序的数列里面是排第几个。询问的时候要把一个询问拆成3个询问，然后再并起来。然后在莫队的时候记录cnti表示当前数字i出现的次数，再开一个bitset，假设当前需要加入bitset的数字是x，我们就令bitset中的第px+cntx位为1，然后++cntx。这样就会使得把三个bitset并起来之后，恰好剩下的就是在三个区间都出现了的数字，也就是需要删除的数字，我们就只需要看并起来之后还剩多少个1就行了

还有一个问题，就是我们必须要对每一个询问(拆之前)都要开一个bitset来存答案，但这会开不下。所以我们就要平衡一下空间复杂度和时间复杂度，把询问分为几块来处理

所以其实程序还可以更快一些，那就是在给所有的块排好序之后再分块。但这样写起来有点麻烦，所以我就懒得卡了(很有道理，但是我没写)   

 

ps: 这里我发现了bitset还是调用他的函数比较快。 

比如now[x]=1 我改成了now.flip(x)  就快了很多。 

 

 

#include<bits/stdc++.h>
using namespace std;  
int const N=100000+10;   
int const B=350; 
int const T=25000;    
struct query{
    int l,r,id;  
    bool operator < (const query &rhs)const{
        if(l/B!=rhs.l/B) return  l/B<rhs.l/B;  
        return r<rhs.r;  
    } 
}q[N];     
int a[N],b[N],n,m,sum[N],ans[N],tot,cnt[N],p[N];  
bitset<N> s[T+1],now;      
void solve(){
    for(int i=1;i<=tot/3;i++) s[i].set();  
    for(int i=1;i<=100000;i++) cnt[i]=0;        
    now.reset();  
    int l=1,r=0;  
    for(int i=1;i<=tot;i++){ 
        while (l>q[i].l){
            --l; now.flip(p[a[l]]+cnt[a[l]]);  cnt[a[l]]++;  
        }  
        while (r<q[i].r){
            ++r; now.flip(p[a[r]]+cnt[a[r]]); cnt[a[r]]++;  
        }    
        while (l<q[i].l){
            cnt[a[l]]--; now.flip(p[a[l]]+cnt[a[l]]);l++;  
        }  
        while (r>q[i].r){
            cnt[a[r]]--; now.flip(p[a[r]]+cnt[a[r]]);r--;  
        }
        s[q[i].id]&=now;    
    } 
    for(int i=1;i<=tot/3;i++)   {
        printf("%d\n",sum[i]-s[i].count()*3);  
    }
} 
int main(){
    scanf("%d%d",&n,&m);  
    for(int i=1;i<=n;i++)  
        scanf("%d",&a[i]),b[i]=a[i];  
    sort(b+1,b+n+1);  
    int k=unique(b+1,b+n+1)-b-1;  
    for(int i=1;i<=n;i++) 
        a[i]=lower_bound(b+1,b+k+1,a[i])-b; 
    memcpy(b,a,sizeof(a));   
    sort(b+1,b+n+1);    
    for(int i=n;i>=1;i--)  
        p[b[i]]=i;  
    int num=(m+T-1)/T; 
    for(int i=1;i<=num;i++){
        tot=0;  
        for(int j=(i-1)*T+1;j<=min(i*T,m);j++){
            int x1,y1,x2,y2,x3,y3;  
            scanf("%d%d%d%d%d%d",&x1,&y1,&x2,&y2,&x3,&y3);  
            q[++tot].l=x1; q[tot].r=y1; q[tot].id=j-(i-1)*T;  
            q[++tot].l=x2; q[tot].r=y2; q[tot].id=j-(i-1)*T;  
            q[++tot].l=x3; q[tot].r=y3; q[tot].id=j-(i-1)*T; 
            sum[j-(i-1)*T]=y1-x1+y2-x2+y3-x3+3; 
        } 
        sort(q+1,q+tot+1);  
        solve();  
    }  
    return 0; 
}