状压DP的核心是状态压缩,用先前求出的状态推导当前状态,用dp[i][j]表示第i层状态j的方案数
这题是上下左右不能有相邻边,因此前一个状态A 与当前状态B 应该是A&B=0 用B的状态子集A的集合之和就可求出
例如 第三层101 的子集是 000 010 因此dp[3][101]=dp[2][000]+dp[2][010];
还有一些小优化在代码注释里
AC代码如下
#include <iostream>
#include <cstdio>
#include <vector>
using namespace std;
const int N = (1<<12)-1;
vector<int> s[N];
int dp[14][N];
int pas[14];
int n,m;
bool cmp(int a)
{
for(int i=0;i<m;i++)
{
if((a>>i)&1 == 1 && (a>>(i+1))&1 == 1) return false;
}
return true;
}//不能有相邻的草地1
int main()
{
scanf("%d%d",&n,&m);
int t;
for(int i=1;i<=n;i++)
for(int j=0;j<m;j++)
{
scanf("%d",&t);
if(t)
pas[i]+=1<<(m-j-1);
}
for(int i=0;i< 1<<m;i++)
{
if(cmp(i))
for(int j=0;j< 1<<m;j++)
{
if(cmp(j)&& (i&j)==0)
s[i].push_back(j);
}
}//将每种状态合法的子集状态求出来
dp[0][0]=1;
for(int i=1;i<=n+1;i++)
{
for(int p=0;p<s[0].size();p++)
{
dp[i][0]+=dp[i-1][s[0][p]];//当前状态为零的时候
}
for(int z=pas[i];z;z=pas[i]&(z-1))//不为零的时候(这里用了求状态子集的方法)
{
if(cmp(z))
for(int p=0;p<s[z].size();p++)
dp[i][z]=(dp[i][z]+dp[i-1][s[z][p]])%100000000;//状态转移
}
}
printf("%d",dp[n+1][0]%100000000);
return 0;
}
国庆做一道简单的状压DP找点自信哈哈