状压妙啊...

本题的主体思路:状压+容斥原理(或状压+数位dp)

记g[i]表示按位与后结果所有位上至少有i个1的方案数

那么根据容斥原理,ans=g[0]-g[1]+g[2]-g[3]+g[4]...

于是如果我们求出了g,就可以求出ans

可是怎么求出g呢

我们记f[i]表示a&i==i这样的a的个数,那么如果i某一位上为1,则a这一位上也为1

于是我们可以枚举所有可能的结果(0-10^6),然后观察这个结果是否是某一个可能结果的子集,如果是的话就累计个数

详细说一下,就是我首先读入所有数据,每读入一个数据x记录f[x]++作为初始值,然后不断更新

在更新的时候,我们首先枚举每一位j,然后枚举1~10^6的所有值i,如果某个值这一位上是1,则更新:

f[i^(1<<j)]+=f[i];

就是去掉j位的个数加上i

什么?怎么证明这样统计是不重不漏的?

首先,我们是按位枚举的,一开始只有初始读入的部分有值,剩下的没有值。那么,当我们枚举第一位时,我们只会更新由初值去掉第一位所能获得的所有值

以此类推,当我们更新第二位时,我们只会更新初值去掉前两位和初值只去掉第二位能获得的所有值

也就是说,我们在更新每一位时,都不会产生重复的状态,都是原来的状态+这一位,所以是不重的,而由于这样的枚举能够遍历所有数位的组合,所以也是不漏的

好,我们处理出了f,接下来?

我们可以枚举所有结果,统计他有几位上是1,那么如果有1位上是1,就会对g[1]产生贡献,等等,以此类推

然后我们再考虑,产生多少贡献?

我们会发现,如果这个结果对应的数有k个,那么答案应为2^k-1(即每个数都有选或不选两种状态,但不能全不选)

所以他产生的贡献就是2^k-1

什么?这种方法的正确性何在?

首先,根据容斥原理,答案的正确性是很显然的

那么我们只需证明g求解的正确性即可

首先回顾一下g的定义:“至少”包含i个1的取法的方案数

也就是说,我所找出的东西数位中1的个数只需>=i即可

那这个是很显然能够保证的

于是为什么不重呢?

由于每个结果互不相同,而我们最后事实上是按结果取的,所以每一种取法都是互不相同的,保证了正确性。

最后代码:

#include <cstdio>
#include <iostream>
#define ll long long
#define mode 1000000007
#define maxx 1000000
ll v[];
ll dp[];
ll f[];
int n;
int main()
{
v[]=;
scanf("%d",&n);
for(int i=;i<=maxx;i++)
{
v[i]=(v[i-]<<)%mode;
}
for(int i=;i<=n;i++)
{
int x;
scanf("%d",&x);
dp[x]++;
}
for(int j=;j<=;j++)
{
for(int i=;i<=maxx;i++)
{
if(((<<j)&i))
{
dp[i^(<<j)]+=dp[i];
}
}
}
ll ans=;
for(int i=;i<=maxx;i++)
{
int cnt=;
for(int j=;j<=;j++)
{
if((<<j)&i)
{
cnt++;
}
}
f[cnt]+=((v[dp[i]]-)%mode+mode)%mode;
}
for(int i=;i<=;i++)
{
if(i%)
{
ans-=f[i];
ans%=mode;
}else
{
ans+=f[i];
ans%=mode;
}
}
printf("%I64d\n",(ans%mode+mode)%mode);
return ;
}

最新文章

  1. Myeclipse中导入项目后java类中汉字注释出现乱码问题(已解决)
  2. Github上十大C#开源项目排行榜
  3. 使用for打印小九九
  4. Gunicorn快速入门
  5. flex4.6事件分派+组件+参数传递
  6. phpstudy虚拟主机配置
  7. LinkButton( 按钮)
  8. git push error: RPC failed; result=56, HTTP code = 0 ,the remote end hung up unexpectedly
  9. 关于js中两种定时器的设置及清除(转载)
  10. spring boot+kafka整合
  11. php 阳历转农历优化版
  12. js的以及前端框架
  13. Scala学习之路 (一)Scala的安装
  14. 【题解】洛谷P1070 道路游戏(线性DP)
  15. C++类模板的三种特化
  16. hihocoder1513
  17. div 隐藏和显示
  18. javascript-实现小抽奖程序
  19. PAT 1055 The World&#39;s Richest
  20. kettle文件以邮件附件的形式发送报告

热门文章

  1. Study 3 —— Python运算符
  2. JAVA锁的优化和膨胀过程
  3. 嵌入式iframe子页面与父页面js通信方式
  4. Python中crypto模块进行AES加密和解密
  5. python初步学习-import和datetime模块
  6. shell 终端常用插件
  7. tidb 架构 ~Tidb学习系列(5)
  8. 根据href给当前导航添加样式
  9. script标签中type为&quot;text/x-template&quot;或&quot;text/html&quot;
  10. java 多线程面试