【算法基础课】状压DP

前置知识：

• 获得 $x$ 的二进制第 $i$ 位： x >> & 1

把二进制看为集合（1表示存在于集合）的话：

• 获得两个集合的交集，并集：x & y , x | y
• 判断其中一个是否是另一个的子集：（x & y） == x

蒙德里安的梦想

• 题解：[AcWing]， 蒙德里安的梦想

• 题意：就是将一个n * m的二维矩阵，分成若干个1 * 2的方格，有多少种分配方式（完全分配）

• 分析：只考虑横着放的。定义 $dp[i,j]$ 为放到第 $i$ 列（放置在 $i,i+1$ 两列的横向方块），且第 $i$ 列的状压方式为 $j$ 时的所有放置数量。那么 $dp[m][0]$ 就是在第i列上不放置 $m,m+1$ 的方块时的答案

• 一列的放置有两个限制：
  1.和上一列的放置方式不能有重叠
  2.留下的间隔必须为偶数，为了可以竖着放置

#include<iostream>

using namespace std;
typedef long long ll;

const int N=1<<12;

ll dp[13][N];
bool st[N];


int main()
{
	int n,m;
	while(cin>>n>>m)
	{
		if(n==0 && m==0) break;
		
		for(int i=0;i<13;i++) for(int j=0;j<N;j++) dp[i][j]=0;
		
		for(int i=0;i< 1<<n ;i++)//预处理数字是否满足第二个限制
		{
			st[i]=1;
			int cnt=0;
			for(int j=0;j<n;j++)
				if(i>>j & 1)
				{
					if(cnt & 1) st[i]=0;
					cnt=0;
				} else cnt++;
			if(cnt & 1) st[i]=0;
		}
		
		dp[0][0]=1;
		for(int i=1;i<=m;i++)
			for(int j=0;j< 1<<n;j++)
				for(int k=0;k< 1<<n;k++)
					if((j & k)==0 && st[j | k])
						dp[i][j]+=dp[i-1][k];
		
		cout<<dp[m][0]<<endl;
		
	}
	return 0;
} 

最短哈密顿距离

题解：【AcWing 91】 最短Hamilton路径 状压DP + 位运算 详解

分析：

• 定义 $dp[i,j]$ 为起点以方式 $j$ 走到点 $i$ 的最短距离。

• 由松弛定理可知，$dis[i,j]$ 可以由 $dis[i,k]$ + $dis[k,j]$ 转移而来，那么可以枚举一个路径的每一种可能，更新 $dis[i,j]$

• 可能会枚举到无效状态，即 $j=k$ 或 i & 1 == 1，但其值都为无穷，可忽略

#include<iostream>
using namespace std;
typedef long long LL;

const int N=22,M=1<<20;
int a[N][N];
int dp[M][N];

int main()
{
	int n; cin>>n;
	
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++) cin>>a[i][j];
	
	for(int i=0;i<M;i++)
		for(int j=0;j<N;j++) dp[i][j]=2e9;
	dp[1][0]=0;
	
	
	for(int i=2;i<1<<n;i++)
		for(int j=0;j<n;j++)
			if(i>>j & 1)
				for(int k=0;k<n;k++)
					if(i>>k & 1)
						dp[i][j]=min(dp[i][j],dp[i-(1<<j)][k]+a[k][j]);
	cout<<dp[(1 << n)-1][n-1]<<endl;
	return 0;
}