状压DP很久之前就看了一下,但是当时我连普通DP都没有做过多少,最近连番做DP题,包括多校训练和平时训练时我都能出一道DP题,于是准备先把DP的专精点满,保证有一些自己很拿手的东西。
状压DP比普通DP多了一个位运算,用二进制的01来表示取或者不取的状态,然后另一维表示处理到哪里。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
#define mod 100000000
using namespace std;
int n,m;
int begin[15];//记录每一行的二进制,判断每一位是否可取
int s[1<<13];
int top;
int dp[1<<15][15];
int jud1(int x)//用来判断相邻位置是否有1
{
if(x&(x<<1))//x<<1到下一位
{
return false;
}
return true;
}
int jud2(int x,int i)//用于判断是否不能放牧的位置被放牧了
{
if(x&begin[i])//每一行的可能取法和图求与,判断是否可取
{
return false;
}
return true;
}
void init()//初始化状态
{
memset(dp,0,sizeof(dp));
top=0;
for(int i=0;i<(1<<m);i++)
{
if(jud1(i))
{
s[top++]=i;//剪枝 只存合适的情况
}
}
}
void solve()
{
int temp;
for(int i=0;i<n;i++)
{
begin[i]=0;
for(int j=0;j<m;j++)
{
scanf("%d",&temp);
if(temp==0)//不可放牧的
{
begin[i]|=(1<<(m-j-1));
}
}
}
init();//初始化
for(int i=0;i<top;i++)
{
if(jud2(s[i],0))
{
dp[s[i]][0]=1;//最初状态
}
}
for(int i=1;i<n;i++)//每一行
{
for(int j=0;j<top;j++)
{
if(!jud2(s[j],i))
{
continue;
}
for(int k=0;k<top;k++)
{
if(jud2(s[k],i-1)&&!(s[k]&s[j]))
{
dp[s[j]][i]=(dp[s[j]][i]+dp[s[k]][i-1])%mod;
}
}
}
}
int ans=0;
for(int i=0;i<top;i++)
{
ans=(ans+dp[s[i]][n-1])%mod;
}
printf("%d\n",ans);
}
int main()
{
while(~scanf("%d%d",&n,&m))
{
solve();
}
return 0;
}