【容斥&&状压】 CF449D Jzzhu and Numbers

题目描述:

在这里插入图片描述


首先我们考虑这么一个问题:
n n n个数中,任意 j j j,使得 a [ j ] & i = i a[j]\&i = i a[j]&i=i,有多少个数满足。
f i f_i fi表示与值为 i i i的数的个数,首先可以得出以下转移:
f i = ∑ j ∉ i f i ∣ ( 1 < < j ) f_i=\sum_{j∉i}f_{i|(1<<j)} fi=j/ifi(1<<j)

但是这样转移很显然会有重复。
一般的转移式这样的:

	for (int i = Max; i >= 0; i--)
	  for (int j = 0; j <= 19; j++)
	    if (((i>>j)&1) == 0)
	      f[i] = (f[i] + f[i|(1<<j)]*x)%P;

但是我们神奇的发现,如果把 i i i j j j的循环翻一下,就可以避免重复,而是计算子集的答案了。

	for (int j = 0; j <= 19; j++)
	  for (int i = Max-1; i >= 0; i--)
	    if (((i>>j)&1) == 0)
	      f[i] = (f[i] + f[i|(1<<j)]*x)%P;

因此就是一个结论:
在这里插入图片描述
这是从我孙女的博客上摘下来的。
所以,对于这道题,做法也十分显然了

我们先用上面的方法做一遍,求出以 i i i为子集的数的个数。
而后,从这些数中任意挑几个数,相与之后也一定以 i i i为子集。
所以方案数是:
g [ i ] = 2 f [ i ] − 1 g[i]=2^{f[i]}-1 g[i]=2f[i]1

在类比之前的问题,只需要把加好变为减号即可。


#include<bits/stdc++.h>
using namespace std;

#define int long long
const int N = 2e6+10,P=1e9+7;
int n,a[N];
int f[1<<20];
int Max;

void dp(int x){
	for (int j = 0; j <= 19; j++)
	  for (int i = Max-1; i >= 0; i--)
	    if (((i>>j)&1) == 0)
	      f[i] = (f[i] + f[i|(1<<j)]*x)%P;
}

int power(int x,int y){
	int sum = 1;
	while (y){
		if (y&1) sum = (sum * x) % P; 
		y>>=1 , x = (x*x)%P;
	}
	return sum;
}

signed main(){
	scanf("%lld",&n);
	for (int i = 1; i <= n; i++) scanf("%lld",&a[i]) , f[a[i]]++;
	Max = (1<<20) - 1;
	dp(1);
	for (int i = 0; i <= Max; i++) f[i] = power(2,f[i])-1;
	dp(-1);
	printf("%lld",(f[0]+P)%P);
}
1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 、4下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。、可私 6信博主看论文后选择购买源代码。 1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 、4下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。、可 6私信博主看论文后选择购买源代码。 1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 、4下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。、可私 6信博主看论文后选择购买源代码。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值