子集dp + 容斥 好题 CF449D

link

题意

在给定的 n 个数中选任意个,使选出的数作 与 运算后为 0 的可选方案数

分析

首先若没有 与 的要求,方案数为 2^n - 1 

把每个数的二进制展开排列整齐,考虑每一列是否合法

考虑容斥 ,g[i]表示选出序列中有 i 列不合法

则 ans=g[0]-g[1]+g[2]-g[3] ......

怎样算出 g ?

考虑状态 x ,为选出序列二进制 与操作后的状态

求出每一个 x ,不用子集dp,用“超集” dp

解决

举个例子

x=(01010)_2

f[x][4]=f[x][3]+f[(11010)_2][3]

而 f[(11010)_2][3]=f[(11010)_2][2]

以此类推进行转移

代码

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int mod=1e9+7,N=1e6+100,M=21;
int n,a[N],f[N*4],ans;
int ksm(int a,int b){
	int ans=1;
	while(b){
		if(b&1) ans*=a;
		ans%=mod;
		a*=a;
		a%=mod;
		b>>=1;
	}
	return ans;
}
signed main(){
	cin>>n;
	for(int i=1;i<=n;i++) {
		cin>>a[i];
		f[a[i]]++;
	}	
	for(int i=0;i<M;i++){
		for(int j=(1<<20)-1;j>=0;j--){
			if(!(j&(1<<i))) f[j]=(f[j]+f[(j|(1<<i))])%mod;
		}
	}
	for(int i=0;i<(1<<20);i++){
		int cnt=0;
		int x=i;
		while(x){
			if(x&1) cnt++;
			x>>=1;
		}
		if(cnt&1) {
			ans=(ans-ksm(2,f[i])+1+mod)%mod;
		}
		else ans=(ans+ksm(2,f[i])-1+mod)%mod;
	}
	ans%=mod;
	cout<<ans;
	
	return 0;
}

posted @ 2024-03-16 12:30  LC_Nocl  阅读(8)  评论(0)    收藏  举报  来源