题目大意: 给出一个序列 a a a,一个子序列是优秀的当且仅当子序列中所有元素 按位与 后的结果为 0 0 0,一个子序列的价值为 φ ( \varphi( φ( 所有元素之和 + 1 ) +1) +1),求出所有优秀子序列的价值之和。
题解
大佬们都说是子集卷积裸题,蒟蒻只能默默懵逼……
每一个数可以被分解为若干个 2 2 2 的幂,如 10 = 2 3 + 2 1 10=2^3+2^1 10=23+21,于是可以将每个数看作一个集合,表示包含哪几个 2 2 2 的幂,那么优秀子序列的定义就变成了:子序列中的所有集合的交集为空。
考虑暴力点的
d
p
dp
dp,设
f
i
f_i
fi 表示选若干个没有交集的集合,使他们的和为
i
i
i 的方案数,对于第
i
i
i 个数(集合),枚举由不在该集合的元素组成的集合
,即枚举那些和他没有交集的数
j
j
j,然后让
f
j
+
a
i
f_{j+a_i}
fj+ai 加上
f
j
f_j
fj。
这样的复杂度是 O ( n max { a i } ) O(n \max\{a_i\}) O(nmax{ai}) 的,卡卡常大概能拿 40 40 40 分吧。
但是这个 d p dp dp 是可以优化的,发现所有相同的 a i a_i ai 是可以一起转移的,设 t i t_i ti 表示数字 i i i 的出现次数,那么转移的时候,就直接让 i i i 枚举数字 1 1 1 ~ 2 × 1 0 5 2\times 10^5 2×105,然后让 f j + i f_{j+i} fj+i 加上 f j f_j fj 乘上 t i t_i ti 即可。
但是这样会出现重复统计,比如说 i = 1 , j = 2 i=1,j=2 i=1,j=2 时, f 3 f_3 f3 会被更新一次, i = 2 , j = 1 i=2,j=1 i=2,j=1 时,又会被更新一次,所以还要限制一下 j < i j<i j<i,这样就不会重复了。
但是要注意,这个限制不能改成 j > i j>i j>i,因为只有小于 i i i 的 j j j 对应的 f j f_j fj 是统计完了的,大于 i i i 的 j j j 还可能被更新。然而为了实现的方便,代码中先枚举的 f i f_i fi,再枚举 t j t_j tj,这样限制就变成了 j > i j>i 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);
}