先说说地图,地图上每一行的01代表一个状态,比如输入样例中的111、010,表示第一行的三个位置都可以种稻子,第二行中间的位置可以种稻子,然后,不能种稻子的地方一定不能种稻子(废话…)
可以种稻子的地方可以选择种也可以选择不种,然后有一个前提条件,就是上下左右相邻的地方不能种稻子。
再说说怎么状态压缩,状态压缩就是把每一个状态压缩成二进制,二进制就是由01组成的,0代表不种,1代表种。二进制就要牵扯到位运算,位运算我就不想说了,百度吧。因此,一串01的二进制数就
可以代表一个状态,例如输入样例第一行是111,那么可以放入第一行的状态有,100、010、001、101、000,因为相邻位置不能放所以只有5种方法,那么第二行就只有2种方法000、010(不考虑其他行)
那么看第一行和第二行(第一行——第二行),100——000,010——000,001——000,101——000,000——000,这是5种对应方法,还可以100——010,001——010,101——010,000——010这是另外的4种对应方法(第一行5种状态对吧?第二行2种状态,按照乘法原理,应该有5*2 = 10种方法,但是111——010是不合法的,因此样例的答案是10-1 = 9)。
dp[i][j]意思是推到第i行状态为j的方案总数。
那么“100——000”即为dp[2][000]可以由dp[1][100]得到,那么dp[2][000] = dp[2][000] + dp[1][100];
那么“010——000”即为dp[2][000]可以由dp[1][010]得到,那么dp[2][000] = dp[2][000] + dp[1][010];
……
以此类推,逐行递推。
总结一下思路:先枚举第一行,把所有可能的状态和第一行的地图对比,如果成功,则在循环里继续枚举第二行,把所有可能的状态和第二行的地图对比,如果成功,再和第一行填入的状态对比,如果又匹配成功,则dp[2][000] = dp[2][000] + dp[1][100];方法数加到第二行。这就是一次循环结束了,从新枚举第二行…
把思路转换成代码
can[]代表可行的状态,稍后解释。cur[i]代表地图的第i行
for(int i=1;i<m;i++)//枚举每一行
{
for(int j=0;j<tot;j++)//对第i行枚举所有可行的状态j
{
if((can[j]&cur[i])==0)//如果状态j和第i行匹配了
{
for(int k=0;k<tot;k++)//枚举第i+1行的所有可行的状态k
{
if(((can[k]&cur[i+1])==0)&&((can[k]&can[j])==0))//状态k和第i+1行匹配且和状态j匹配
dp[i+1][can[k]] = dp[i+1][can[k]]+dp[i][can[j]];//状态数相加
}
}
}
}
这样核心代码就实现了。
有一个小方法,就是枚举可行状态的时候,假如一行是8列,不必从00000000枚举到11111111,这样很麻烦,所以要预处理。
就是在一开始把,一行的可行状态先求出来就拿“11111111”来说,这肯定是不可能的,因为有相邻的1,所以在一开始就可以舍弃掉。怎么做呢?
假如一行是8列,先从00000000枚举到11111111,对于每一个状态把它左移1位,再和他自己&运算,假如结果>0,就说明有有相邻的1,举个简单的例子:
01011要判断有没有相邻的1,if(((01011<<1) & (01011)) > 0 )则有相邻的1,(01011<<1) & (01011) 就是 010110和01011按位且运算,这两个红色地方1&1 == 1,因此结果大于0。
算法实现代码:
#include <cstring>
#include <cstdio>
#include <iostream>
#include <algorithm>
#define mod 100000000
using namespace std;
int dp[13][1<<12], cur[13];
int can[1<<12], tot, m, n;
int main() {
while (scanf("%d%d", &m, &n)) {
tot = 0;
memset(cur, 0, sizeof(cur));
memset(dp, 0, sizeof(dp));
for (int i = 0; i < (1<<n); i++) {
if ((i&(i<<1)) == 0) can[tot++] = i;
}
for (int i = 1; i <=m; i++) {
for (int j = 0; j < n; j++) {
int num;
scanf("%d", &num);
if (num == 0) cur[i] = (cur[i] | (1 << j));
}
}
for (int i=0; i < tot; i++) {
if ((cur[1]&can[i]) == 0) dp[1][can[i]] = 1;
}
for (int i = 1; i < m; i++) {
for (int j = 0; j < tot; j++) {
if ((can[j]&cur[i]) == 0) {
for (int k = 0; k < tot; k++) {
if (((can[k]&cur[i+1])==0) && ((can[k]&can[j])==0)) {
dp[i+1][can[k]] = dp[i+1][can[k]] + dp[i][can[j]];
}
}
}
}
}
int ans = 0;
for (int i = 0; i < tot; i++) {
ans += dp[m][can[i]];
ans = ans % mod;
}
cout << ans << endl;
}
return 0;
}