POJ 3254 状态压缩DP

【题目链接】
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;
}