NOIP模拟：Subset（分块）

最新推荐文章于 2018-08-09 15:49:02 发布

原创最新推荐文章于 2018-08-09 15:49:02 发布 · 569 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#分块

分块专栏收录该内容

23 篇文章

订阅专栏

本文介绍了一种集合操作与查询的算法实现，包括添加、删除元素及查询满足特定条件元素数量的功能。通过分块思想优化处理流程，适用于大量数据查询场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Description
一开始你有一个空集，集合可以出现重复元素，然后有 Q 个操作：
1. add s
  在集合中加入数字 s 。
2. del s
  在集合中删除数字 s 。保证 s 存在。如果有多个 s，只删除一个即可。
3. cnt s
  查询满足 a&s=a 条件的 a 的个数。
Input
第一行一个整数 Q 接下来 Q 行，每一行都是 3 个操作中的一个。
Output
对于每个 cnt 操作输出答案。
Sample Input
7
add 11
cnt 15
add 4
add 0
cnt 6
del 4
cnt 15
Sample Output
1
2
2
Hint
对于 30% 的数据满足： $1≤n≤1000$ ；
对于 100% 的数据满足： $1≤n≤200000；0＜s＜2^{16}$ 。
题解：分块
容易发现对于每一个cnt操作的s，满足的数a前后半段互相之间无影响。
考虑将每一个输入的a分为( $2^8*k_1+k_2(k_2<2^8)$ ),即 $k_1$ 为a的前八位， $k_2$ 为a的后八位。
维护 $f[pre][suf]$ ,表示满足 $k_1$ & $pre=k_1$ 中 $suf=k_2$ 的个数，加入a时更新满足条件的区间，查询s时，先找到s的前缀 $pre_1$ ,由以上性质可得 $a[pre_1]$ 中所有 $a$ 的前半段都满足,这时只需在此区间中循环查找suf满足的个数即可。

时间复杂度: $O(n\sqrt{n})$
- Code:

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
using namespace std;

const int Maxn=(1<<8);

inline int read()
{
    char ch=getchar();int i=0,f=1;
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){i=(i<<3)+(i<<1)+ch-'0';ch=getchar();}
    return i*f;
}

int n,a[Maxn+50][Maxn+50];
char ch[5];

int buf[30];
inline void W(int x)
{
    if(x==0)
    {
        putchar('0');
        return;
    }
    while(x)buf[++buf[0]]=x%10,x/=10;
    while(buf[0])putchar(buf[buf[0]--]+'0');
    return;
}

int main()
{
    n=read();
    while(n--)
    {
        scanf("%s",ch+1);
        int s=read();
        int t=s>>8,res=s-(t<<8);
        switch(ch[1])
        {
            case 'a':
                {
                    for(int i=0;i<Maxn;i++)
                    {
                        if((t&i)==t)a[i][res]++;
                    }
                    break;
                }
            case 'c':
                {
                    int r=0;
                    for(int i=0;i<Maxn;i++)
                    {
                        if((i&res)==i)r+=a[t][i];
                    }
                    W(r);
                    putchar('\n');
                    break;
                }
            case 'd':
                {
                    for(int i=0;i<Maxn;i++)
                    {
                        if((t&i)==t)a[i][res]--;
                    }
                    break;
                }
        }
    }
    return 0;
}