Corn Fields

Corn Fields

农场主John新买新牧场，这块牧场被划分成M行N列(1 ≤ M ≤ 12; 1 ≤ N ≤ 12)，每一格都是一块正方形的土地。John打算在牧场上的某几格里种上美味的草，供他的奶牛们享用。

遗憾的是，有些土地相当贫瘠，不能用来种草。并且，奶牛们喜欢独占一块草地的感觉，于是John不会选择两块相邻的土地，也就是说，没有哪两块草地有公共边。

John想知道，如果不考虑草地的总块数，那么，一共有多少种种植方案可供他选择？（当然，把新牧场完全荒废也是一种方案）

Input

第一行：两个整数M和N，用空格隔开。

第2到第M+1行：每行包含N个用空格隔开的整数，描述了每块土地的状态。第i+1行描述了第i行的土地，所有整数均为0或1，是1的话，表示这块土地足够肥沃，0则表示这块土地不适合种草。

Output

一个整数，即牧场分配总方案数除以100,000,000的余数。

Sample Input

2 3
1 1 1
0 1 0

Sample Output

9

思路

蒙德里安那道题的一种应用, 依然是用二进制来表示该列上哪个位置种草。
总思路就是枚举所有可行状态, 然后计算转移, 最后得出结果方案数。

直接迁移思考路径, 先考虑怎样的一列是成立的：
题目要求任意两块草坪不能相邻, 相比于蒙德里安要简单一些, 只需要用个计数来判断是否存在相邻草坪即可

for(int i = 0; i < 1 << n; i++)
{
	bool cnt = false, flag = false;
	for(int j = 0; j < n; j++)
	{
		if(i >> j & 1) // 这一位是1, 种草坪
		{
			if(cnt) // 如果相邻一位也是1
			{
				flag = true;
				break;
			}
			else cnt = true;
		}
		else cnt = false;
	}
	if(!flag) ...; // 如果成立再进行下一步操作
}

接着考虑某一个状态可以从哪个状态中转移过去：
设 第i列 当前状态为 j, i - 1列 状态为 k, 初始状态可以转移

当出现相邻1时, 不能转移（不成立）：

故需要满足 k & j == 0。

题目上还有一个条件, 有的地方是不能种草坪的, 所以在判断某列是否成立时, 需要加上判定：

for(int k = 1; k <= n; k++)
for(int i = 0; i < 1 << m; i++)
{
	bool cnt = false, flag = false;
	for(int j = 0; j < m; j++)
	{
		if(i >> j & 1) // 这一位是1, 种草坪
		{// k 表示某一列, 这里无论j是否能将题目中数组对应上二进制的位置, 对结果没有影响, 因为需要在意的只是相对位置
			if(cnt || g[k][j] == 0) // 如果相邻一位也是1
			{
				flag = true;
				break;
			}
			else cnt = true;
		}
		else cnt = false;
	}
	if(!flag) ...; // 如果成立再进行下一步操作
}

最后的结果怎么求：
遍历最后一列的所有状态, 将其加起来即可。

for(int i = 0; i < 1 <<m; i++)
{
	if(check(i))
		res += f[n][i];
}

对于当前列状态是否有效还有一种更巧判断方法：
判断是否存在相邻1可以 (x&(x >> 1)) && (x&(x <<1)) 若左/右移后再&结果不为0 则说明存在相邻, 整个表达式返回0时才说明状态存在。

判断是土地是否肥沃可以先将输入处理成同样的状态表示 map[i] , 若状态 x 可以存在, 说明 x & map[i] == x。

代码

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>

using namespace std;
//------- Coding Area ---------//
const int N = 15,M = 1 << N, MOD = 100000000; // M是状态
int n, m;
int g[N][N], f[N][M], map[N];
long long res;


bool check(int k,int x) // 第k列的状态x
{
    bool flag1 = (x & (x >> 1)) && (x & (x << 1));
    bool flag2 = (x & map[k]) != x;
    return !(flag1 || flag2);
}

int main()
{
    cin >> m >> n;
    for (int i = 0; i < m; i++)
        for (int j = 0; j < n; j++)
            cin >> g[j][i];
   // 处理输入
    for (int i = 0; i < n; i++)
    {
        for(int j = 0; j < m; j++)
            map[i] = (map[i] << 1) + g[i][j];
    }
    
    // 起始状态计算
    for (int i = 0; i < 1 << m; i++)
        if (check(0,i)) f[0][i] = 1;
    
    // 开始DP
    for (int i = 1; i < n; i++)
    {
        for (int j = 0; j < 1 << m; j++)
        {
            for (int k = 0; k < 1 << m; k++)
            {
                if ((j & k) == 0 && check(i, j) && check(i - 1, k))
                    f[i][j] = (f[i][j] + f[i - 1][k]) % MOD;
            }
        }
    }

	// 找出答案
    for (int i = 0; i < 1 << m; i++)
    {
        if (check(n - 1, i)) res  = (res + f[n - 1][i]) % MOD;
    }
    cout << res << endl;
    return 0;
}