【题目链接】
http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=21819
【解题报告】
状态压缩dp入门题。
因为一层最多只有12个节点,容易想到我们用一串二进制数表示它,然后映射到对应的十进制数上。
那么我们可以一层一层的转移。
这里转移的时候需要满足两个条件。
一个状态j,需要满足!(j&(j>>2))。这说明,j表示的放奶牛的方案,不存在相邻的两个格子放奶牛。
另外,如果可以从dp[i-1][k]转移到dp[i][j],那么应该有!(j&k)表示j和k不存在相同位置上有奶牛。
最后考虑一下时间复杂度,最多需要1<<12个十进制数来保存状态,裸dp的时间复杂度为O( M*2^M*2^N) ),也就是大概4000*4000*12次运算,还是可以保证通过的。
【参考代码】
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int mod=100000000;
int N,M;
int fer[14][14];
int dp[14][10000];
bool check( int ans,int i )
{
if( ans&(ans>>1) )return false;
int k=0;
while( ans )
{
if( ans%2 && !fer[i][k] )return false;
ans/=2;k++;
}
return true;
}
void init()
{
for( int i=0; i<M; i++ )
for( int j=0; j<N; j++ )
scanf( "%d",&fer[i][j] );
memset( dp,0,sizeof dp );
for( int j=0; j<(1<<N); j++ ) //第一层赋值
{
if( check( j,0 ) )dp[0][j]=1;
}
}
int main()
{
while( ~scanf( "%d%d",&M,&N ) )
{
init();
for( int i=1; i<M; i++ )
for( int j=0; j<(1<<N); j++ ) //j当前这种状态
{
if( !check(j,i) )continue;
for( int k=0; k<(1<<N); k++ )
{
if( dp[i-1][k] && !( k&j ) ) dp[i][j]+=dp[i-1][k];
dp[i][j]%=mod;
}
}
int ans=0;
for( int i=0; i<(1<<N); i++ )
{
ans+=dp[M-1][i];
ans%=mod;
}
printf( "%d\n",ans );
}
return 0;
}