15.蒙德里安的梦想 状态压缩DP

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 typedef long long ll;
 4 const int N = 13, M = 1 << N;
 5 ll dp[N][M];
 6 bool st[M];
 7 int main() {
 8     int n, m;
 9     while (cin >> n >> m && n != 0) {
10         //预处理一下，所有的状态是否不存在连续奇数个0
11         for (int i = 0; i < (1 << n); i++) { //遍历所有的状态
12             st[i] = true; //假设是成立的
13             int cnt = 0; //cnt是当前这一段连续0的个数
14             for (int j = 0; j < n; j++) { //n是这个二进制数的位数
15                 if (i >> j & 1) { //如果当前这一位是1，说明上一段已经截止了
16                     if (cnt & 1) { //判断上一段连续的0是否是奇数个，若是
17                         st[i] = false; //说明第i个状态时不合法的，存在了连续奇数个0
18                     }
19                     cnt = 0; //遇到1以后，连续奇数个0已经结束了，新开始一段
20                 } else {
21                     cnt++;
22                 }
23             }
24             if (cnt & 1) { //然后判断最后一段0的个数
25                 st[i] = false;
26             }
27         }
28         //然后是dp的过程
29         memset(dp, 0, sizeof dp); //把dp数组置为0
30         dp[0][0] = 1; //边界情况
31         for (int i = 1; i <= m; i++) { //枚举所有的列
32             for (int j = 0; j < (1 << n); j++) { //枚举第i列的所有状态
33                 for (int k = 0; k < (1 << n); k++) { //再枚举第i - 1列的所有状态
34                     if ((j & k) == 0 && st[j | k]) {
35                         dp[i][j] += dp[i - 1][k];
36                     }
37                 }
38             }
39         }
40         cout << dp[m][0] << endl; //答案
41     }
42     return 0;
43 }