C. Boboniu and Bit Operations(状压)

一个超级显然的状压

可惜我昨天没打

这题有个很简单的思维解法,想看戳我哦,比状压简单很多

言归正传,反正思维的解法我没想到,用状压草过去的,状压毕竟具有普适性

$定义dp[i][j]为构造好第i个c_i后,前i个c元素相或运算是否可能得到二进制j$

$这个j的范围就是[0,1<<9]嘛,因为或和与不会改变二进制范围$

$那么对于当前的a_i,枚举一个b_q构造成此时的c_i$

$再来一层for循环枚举与运算完c_i后的二进制j$

$再来一层for循环枚举二进制d表示要从dp[i-1][d]转移而来$

于是转移方程是

$dp[i][j]=dp[i][j]$ | $dp[i-1][d]$

中间的竖线表示或运算

$意思是,只要存在一个d使得dp[i-1][d]可行,那么dp[i][j]也可行$

于是我们可以写下如下代码

dp[0][0]=1;
for(int i=1;i<=n;i++)
for(int q=1;q<=m;q++)
{
	int s = (a[i]&b[q]) ;//当前构造的c[i]
	for(int j=0;j<=(1<<9);j++)//本状态 
	{
		 if( (j&s)!=s )	continue;//本状态由于或上了c[i],二进制一定包含c[i]
		 for(int d=0;d<=j;d++)//上一状态d|s=j,所以不可能大于j
		 {
		 	if( (d|s)!=j )	continue;
		 	dp[i][j]|=dp[i-1][d];
		 }
	}
}

拿起小拇指一算复杂度,哇呜呜

$2^9\ast 2^9\ast 200\ast 200=........$

$但是仔细看,枚举d的那个循环,很多工作是不需要的$

$现在已知s和j,我们需要枚举的都是s|d==j的那些d$

d或运算s变成j

说明$s$只贡献了自己的某个子集给$j$,设贡献了$y$

然后$y$一定是$s$的子集,所以直接枚举子集就好了

这样来枚举数$v$的子集

for(int w=v;w;w=(w-1)&v )

#include <bits/stdc++.h>
using namespace std;
int n,m;
int f[209][1<<10],a[209],b[209];
int main()
{
	//定义f[i][j]表示前i个c[i]或运算能否得到二进制j
	cin >> n >> m;
	for(int i=1;i<=n;i++)	cin >> a[i];
	for(int i=1;i<=m;i++)	cin >> b[i];
	f[0][0] = 1;
	for(int i=1;i<=n;i++) 
	for(int j=1;j<=m;j++)//枚举这次选择的b[j]
	{
		int v = a[i]&b[j];//当前选的是v 
		for(int u=0;u<(1<<9);u++)//枚举当前状态f[i][u]
		{
			if( (u&v)!=v )	continue;		
			for(int w=v;w;w=(w-1)&v )//枚举v的子集贡献
				f[i][u] |= f[i-1][u^w];	
			f[i][u] |= f[i-1][u];
		}
	} 
	for(int i=0;;i++)
		if( f[n][i] ) { cout << i; return 0; }
}