很有意思的一题,题意是找到有多少个子序列,使得子序列每个元素的&值为0。
我们先考虑,如果每个都是0,那么方案显然是2^n-1,但是可能算多了。分别考虑第1位到第20位一定是1的情况,因为是&运算,选的一定是该位为1。然后统计两位的,三位的。。。答案就是 ∑ -1^d(x) * ( 2^f(x )-1 ) ,x的范围的[0,1<<20),d(x)表示x改为二进制有多少个0,f(x)表示的是满足 ai&x == x的数字个数。
现在的问题关键是怎么计算f(x),因为直接暴力显然是不靠谱的,我们用dp[i][j] 表示某个状态,j是一个二进制数,我们不妨记j的低 i 位为x0,其他高位为x1,即 x1x0 == j ,x0中的1表示某个数的此位一定是1,0表示这个数的这一位可以是任意的数;而在x1中,1表示这个数的此位一定是1,0表示这个数的此位一定是0,满足这种情况的数字个数。显然dp[19][j] = f(j)。我们可以很容易得到转移方程(具体见代码)。
//#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<cmath>
#include<cctype>
#include<string>
#include<algorithm>
#include<iostream>
#include<ctime>
#include<map>
#include<set>
using namespace std;
#define MP(x,y) make_pair((x),(y))
#define PB(x) push_back(x)
typedef long long LL;
//typedef unsigned __int64 ULL;
/* ****************** */
const LL INF=1LL<<60;
const double INFF=1e100;
const double eps=1e-8;
const int mod=1000000007;
const int NN=100005;
const int MM=1000010;
/* ****************** */
int dp[20][1<<20];
int a[1<<20];
int q_pow(int x,int n,LL mod)
{
LL ans=1,xx=x;
for(;n>0;n>>=1)
{
if(n&1)ans=ans*xx%mod;
xx=xx*xx%mod;
}
return (int)ans;
}
int main()
{
int n,i,j,mask=1<<20;
int ans,t,tt;
scanf("%d",&n);
for(i=1;i<=n;i++)
scanf("%d",&a[i]);
memset(dp,0,sizeof(dp));
for(i=1;i<=n;i++)
{
if(a[i]&1)
{
dp[0][ a[i] ]++;
dp[0][ a[i]^1 ]++;
}
else
{
dp[0][ a[i] ]++;
}
}
for(i=0;i<19;i++)
for(j=0;j<mask;j++)
{
t=1<<(i+1);
if(j&t)
{
dp[i+1][j]+=dp[i][j];
dp[i+1][ j-t ]+=dp[i][j];
}
else
{
dp[i+1][j]+=dp[i][j];
}
}
ans=0;
for(j=0;j<mask;j++)
{
t=1;
for(i=0;i<20;i++)
if((1<<i)&j)
t=-t;
tt=q_pow(2,dp[19][j],mod);
tt--;
if(tt<0)tt+=mod;
ans+=t*tt;
if(ans>=mod)ans-=mod;
if(ans<0)ans+=mod;
}
printf("%d\n",ans);
return 0;
}