[POJ3254]Corn Fields（状压DP）

题目链接：http://poj.org/problem?id=3254

给n*m的格子图，图上有地方有草有地方没草，FJ希望把牛放到有草的地方，但是上下左右不能有相邻的牛在一起。问有多少中分配方式，结果对100000000取模。

状压DP，设dp(i,j)为第i行时状态j，由于m<=12所以只需要12个二进制位来表示当前行是如何占用的即可。

用一个一位数组存图的状态，初始化的时候枚举第一行，遇到这一行没有相邻1出现并且这个状态都是草地的情况，可以把这个状态初始化为1。

接下来枚举2~n行i，和i行的j个状态（有1<<m个）。如果不满足：1、这一行没有相邻1出现。2、这个状态都是草地，则这个状态下方案数为0。否则，枚举前一行的状态k，查看i-1行的状态k是否可以转移到i行的状态j。如果可以的话就累加取模。

/*
━━━━━┒ギリギリ♂ eye！
┓┏┓┏┓┃キリキリ♂ mind！
┛┗┛┗┛┃＼○／
┓┏┓┏┓┃ /
┛┗┛┗┛┃ノ)
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┃┃┃┃┃┃
┻┻┻┻┻┻
*/
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
using namespace std;
#define fr first
#define sc second
#define cl clear
#define BUG puts("here!!!")
#define W(a) while(a--)
#define pb(a) push_back(a)
#define Rint(a) scanf("%d", &a)
#define Rll(a) scanf("%I64d", &a)
#define Rs(a) scanf("%s", a)
#define Cin(a) cin >> a
#define FRead() freopen("in", "r", stdin)
#define FWrite() freopen("out", "w", stdout)
#define Rep(i, len) for(int i = 0; i < (len); i++)
#define For(i, a, len) for(int i = (a); i < (len); i++)
#define Cls(a) memset((a), 0, sizeof(a))
#define Clr(a, x) memset((a), (x), sizeof(a))
#define Full(a) memset((a), 0x7f7f7f, sizeof(a))
#define lrt rt << 1
#define rrt rt << 1 | 1
#define pi 3.14159265359
#define RT return
#define lowbit(x) x & (-x)
#define onecnt(x) __builtin_popcount(x)
typedef long long LL;
typedef long double LD;
typedef unsigned long long ULL;
typedef pair<int, int> pii;
typedef pair<string, int> psi;
typedef pair<LL, LL> pll;
typedef map<string, int> msi;
typedef vector<int> vi;
typedef vector<LL> vl;
typedef vector<vl> vvl;
typedef vector<bool> vb;

const int mod = 100000000;
const int maxn = 15;
const int maxm = 1 << maxn;
int G[maxn];
int dp[maxn][maxm];
int n, m, mm;

bool ok(int x, int y) {
    if((G[x] & y) != y) return 0;
    if(((y << 1) & y) != 0) return 0;
    return 1;
}

int main() {
    // FRead();
    int x;
    while(~Rint(n) && ~Rint(m)) {
        Cls(G); Cls(dp);
        mm = 1 << m;
        For(i, 1, n+1) {
            Rep(j, m) {
                Rint(x);
                G[i] <<= 1;
                G[i] |= x;
            }
        }
        Rep(j, mm) if(ok(1, j)) dp[1][j] = 1;
        For(i, 2, n+1) {
            Rep(j, mm) {
                if(!ok(i, j)) continue;
                Rep(k, mm) {
                    if(j & k) continue;
                    dp[i][j] = (dp[i][j] + dp[i-1][k]) % mod;
                }
            }
        }
        int ret = 0;
        Rep(i, mm) ret = (ret + dp[n][i]) % mod;
        printf("%d
", ret);
    }
    RT 0;
}