POJ - 3254 Corn Fields 状压DP

题意：

把奶牛放在1的位置上，且两个奶牛之间不能相邻

求出一行当中满足两两不相邻的状态  ->  当前状态 & 当前状态左移一位， 如果结果为零则满足条件，否则不满足 存入st数组

例： 101， 1010  |   11， 110

将每一行的1，0转成二进制存到一个数里，这里反着处理一下，010 表示 第一三草坪可以放牛 存入a数组

这样处理可以保证，当st & a = 0 的时候，说明st这个状态是可以放在这一行的

然后 dp 处理一下

st[j] & st[k] == 0 保证竖直方向不相邻

dp[i][j] += dp[i-1][k]  dp[i][j] 表示第i行放第j个满足两两不相邻的状态时候能达到的方案数 

 

 

AC代码

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<queue>
#include<string>
#include<map>
#include<cmath>
#define INF 0x3f3f3f3f
#define rep(i,a,n) for(int i=a;i<n;i++)
#define per(i,a,n) for(int i=n-1;i>=a;i--)
#define fori(x) for(int i=0;i<x;i++)
#define forj(x) for(int j=0;j<x;j++)
#define memset(x,y) memset(x,y,sizeof(x))
#define memcpy(x,y) memcpy(x,y,sizeof(y))
#define sca(x) scanf("%d", &x)
#define scas(x) scanf("%s",x)
#define sca2(x,y) scanf("%d%d",&x,&y)
#define sca3(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define scl(x) scanf("%lld",&x)
#define scl2(x,y) scanf("%lld%lld",&x,&y)
#define scl3(x,y,z) scanf("%lld%lld%lld",&x,&y,&z)
#define pri(x) printf("%d\n",x)
#define pri2(x,y) printf("%d %d\n",x,y)
#define pris(x) printf("%s\n",x)
#define prl(x) printf("%lld\n",x)
//#include <bits/stdc++.h>

typedef long long ll;
const int maxn=1e6+7;
const int mod=100000000;
const double eps=1e-8;
// const double pi = acos(-1);

using namespace std;

int dp[13][1<<13];
int st[1<<13];
int a[13];

int main()
{
    int n,m;
    sca2(n,m);
    memset(a,0);
    memset(dp,0);
    rep(i,0,n)
    {
      rep(j,0,m)
      {
        int temp;
        sca(temp);
        if(temp == 0)
        {
          a[i] += ( 1<<(m-1-j) );
        }
      }
    }
    int cnt = 0;
    rep(i,0,(1<<m))
    {
      if((i&(i<<1)) == 0)
      {
        st[cnt++] = i;
      }
    }
    rep(i,0,cnt)
    {
      if((st[i] & a[0]) == 0)
      {
        dp[0][i] = 1;
      }
    }
    rep(i,1,n)
    {
      rep(j,0,cnt)
      {
          if(st[j] & a[i])  continue;
          rep(k,0,cnt)
          {
            if(st[k] & a[i-1])  continue;
            if(!(st[k] & st[j]))
            {
              dp[i][j] += dp[i-1][k];
            }
          }
      }
    }
    int ans = 0;
    rep(i,0,cnt)
    {
      if(!(a[n-1] & st[i]))
        ans = (ans + dp[n-1][i])%mod;
    }
    pri(ans);
    return 0;
}