NOI Online3 提高组 优秀子序列 题解

题目传送门

题目大意： 给出一个序列 $a$，一个子序列是优秀的当且仅当子序列中所有元素 按位与 后的结果为 $0$，一个子序列的价值为 $\phi ($ 所有元素之和 $+1)$，求出所有优秀子序列的价值之和。

题解

大佬们都说是子集卷积裸题，蒟蒻只能默默懵逼……

每一个数可以被分解为若干个 $2$ 的幂，如 $10=2^3+2^1$，于是可以将每个数看作一个集合，表示包含哪几个 $2$ 的幂，那么优秀子序列的定义就变成了：子序列中的所有集合的交集为空。

考虑暴力点的 $dp$，设 $f_i$ 表示选若干个没有交集的集合，使他们的和为 $i$ 的方案数，对于第 $i$ 个数（集合），枚举由不在该集合的元素组成的集合，即枚举那些和他没有交集的数 $j$，然后让 $f_{j+a_i}$ 加上 $f_j$。

这样的复杂度是 O ( n max ⁡ { a i } ) O(n \max\{a_i\}) O(nmax{ai​}) 的，卡卡常大概能拿 $40$ 分吧。

但是这个 $dp$ 是可以优化的，发现所有相同的 $a_i$ 是可以一起转移的，设 $t_i$ 表示数字 $i$ 的出现次数，那么转移的时候，就直接让 $i$ 枚举数字 $1$ ~ $2\times 10^5$，然后让 $f_{j+i}$ 加上 $f_j$ 乘上 $t_i$ 即可。

但是这样会出现重复统计，比如说 $i=1,j=2$ 时，$f_3$ 会被更新一次，$i=2,j=1$ 时，又会被更新一次，所以还要限制一下 $j<i$，这样就不会重复了。

但是要注意，这个限制不能改成 $j>i$，因为只有小于 $i$ 的 $j$ 对应的 $f_j$ 是统计完了的，大于 $i$ 的 $j$ 还可能被更新。然而为了实现的方便，代码中先枚举的 $f_i$，再枚举 $t_j$，这样限制就变成了 $j>i$。

代码如下：

#include <cstdio>
#include <algorithm>
using namespace std;
#define maxn 300010
#define mod 1000000007
#define ll long long

int n;
ll t[maxn],f[maxn];
int prime[maxn],tt=0,phi[maxn];
bool v[maxn];
void work()
{
	phi[1]=1;
	for(int i=2;i<=maxn-10;i++)
	{
		if(!v[i])prime[++tt]=i,phi[i]=i-1;
		for(int j=1;j<=tt&&i*prime[j]<=maxn-10;j++)
		{
			v[i*prime[j]]=true;
			if(i%prime[j]==0)
			{
				phi[i*prime[j]]=phi[i]*prime[j];
				break;
			}
			phi[i*prime[j]]=phi[i]*(prime[j]-1);
		}
	}
}
inline char cn()
{
	static char buf[1000010],*p1=buf,*p2=buf;
	return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
void read(int &x)
{
	x=0;int f1=1;char ch=cn();
	while(ch<'0'||ch>'9'){if(ch=='-')f1=-1;ch=cn();}
	while(ch>='0'&&ch<='9')x=x*10+(ch-'0'),ch=cn();
	x*=f1;
}

int main()
{
	read(n); work(); int ma=0;
	for(int i=1,x;i<=n;i++)read(x),t[x]++,ma=max(x,ma);
	int limit=1;while(limit<=ma)limit<<=1;f[0]=1;
	for(int i=0;i<limit;i++)if(f[i])
	{
		int p=(limit-1)^i;//求出集合i的补集，下面的j是在枚举p的子集，这里的j定义和上面是一样的
		for(int j=p;j>i;j=(j-1)&p)(f[i|j]+=f[i]*t[j]%mod)%mod;
	}
	ll ans=1;for(int i=1;i<limit;i++)(ans+=f[i]*phi[i+1]%mod)%mod;
	while(t[0]--)ans=2*ans%mod; printf("%lld",ans);
}