hiho_1048_状态压缩2
Posted 农民伯伯-Coding
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了hiho_1048_状态压缩2相关的知识,希望对你有一定的参考价值。
题目大意
用1x2的单元拼接出 NxM的矩形,单元可以横放或者纵放,N < 1000, M <= 5. 求不同的拼接方案总数。
分析
计算机解决问题的基本思路:搜索状态空间。如果采用dfs进行搜索,则可以将当前搜索的位置(i, j)作为状态,而不保存棋盘的占用情况,每次扩展时进行判断是否可行。这样可以得到最后的结果,但当然会超时。
(1) 参考hiho_1048中的提示,可以知道在搜索的时候按照固定的顺序可以大幅度剪枝(通过将无序的搜索变成有序的搜索来实现)。然后可以对中间结果进行保存,即记忆化搜索。状态dp[i][j][s1][s2] 表示求第i行,第j列位置开始扩展,且第i行的占用状态为s1, 第i+1行的占用状态为s2 时,填满剩余的空格的方案总数。
(2)将记忆化搜索直接使用动态规划来实现,状态 dp2[i][k1][k2] 表示前i-1行都已经填满了,此时第i行的占用情况的二进制表示转换为十进制为 k1,第i+1行的占用情况的二进制表示转换为十进制为 k2,达到此时棋盘占用情况的摆放方案总数。
这样,可以进行状态转移时,枚举第i行的占用情况k1,和第i+1行的占用情况k2,看是否可以进行扩展,blabla...
实现
#include<iostream> #include<stdio.h> #include<string.h> using namespace std; #define BIT(i, j) (i & (1 << (j - 1))) const int mod = 1000000007; int dp[1005][5][35][35]; /*记忆化搜索 n,m为边界n行,m列。 求第i行,第j列位置开始扩展,且第i行的占用状态为s1, 第i+1行的占用状态为s2 时,填满剩余的空格的方案总数 */ int dfs(int n, int m, int i, int j, int s1, int s2){ if (dp[i][j][s1][s2] != -1) return dp[i][j][s1][s2]; int result = 0; //边界!! 无法继续扩展时,取值为1,表示所有扩展到此的合理情况的初始均为1 if (i > n || j > m || s1 >= (1 << m) || s2 >= (1 << m)) return 1; //如果第i行位置j(从右向左)被占用,且j < m,向左扩展 if (BIT(s1, j) && j < m) result = dfs(n, m, i, j + 1, s1, s2); //如果第i行位置j(从右向左)被占用,且j = m,向下一行扩展 else if (BIT(s1, j) && j == m) result = dfs(n, m, i + 1, 1, s2, 0); //如果第i行位置j没有被占用,且无法向左或者向下扩展 else if (BIT(s1, j) == 0 && (j == m || BIT(s1, j + 1)) && (i == n || BIT(s2, j))) result = 0; //如果第i行位置j没有被占用,且只能向左扩展 else if (BIT(s1, j) == 0 && j < m && BIT(s1, j + 1) == 0 && (i == n || BIT(s2, j))) result = dfs(n, m, i, j, s1 | (1 << (j - 1)) | (1 << j), s2); //如果第i行位置j没有被占用,且只能向下扩展 else if (BIT(s1, j) == 0 && (j == m || BIT(s1, j + 1)) && (i < n && BIT(s2, j) == 0)) result = dfs(n, m, i, j, s1 | (1 << (j - 1)), s2 | (1 << (j - 1))); //如果第i行位置j没有被占用,且可以向左或向下扩展 else if (BIT(s1, j) == 0 && j < m && BIT(s1, j + 1) == 0 && i < n && (BIT(s2, j) == 0)) result = dfs(n, m, i, j, s1 | (1 << (j - 1)) | (1 << j), s2) + dfs(n, m, i, j, s1 | (1 << (j - 1)), s2 | (1 << (j - 1))); //记忆化保存 return dp[i][j][s1][s2] = result % mod; } //动归数组, dp2[i][k1][k2] 表示前i-1行都已经填满了,此时第i行的占用情况的二进制表示转换为十进制为 k1, // 第i+1行的占用情况的二进制表示转换为十进制为 k2 int dp2[1005][35][35]; int main(){ int n, m; scanf("%d %d", &n, &m); /* 记忆化搜索 memset(dp, -1, sizeof(dp)); int result = dfs(n, m, 1, 1, 0, 0); printf("%d\n", result); */ //动态规划 memset(dp2, 0, sizeof(dp2)); //初始状态,第1行未被占用,只有1种可能。 dp2[1][0][0] = 1; int mx = 1 << m; for (int i = 1; i <= n; i++){ for (int j = 1; j <= m; j++){ //枚举第i行的所有情况(占用情况的二进制表示转换为十进制) for (int k1 = 0; k1 < mx; k1++){ //如果第i行的第j个位置(从右向左)为0,则可以从这里开始扩展 if ((k1 & (1 << (j - 1))) == 0){ //看是否可以横向扩展 if (j < m && ( ((1 << j) & k1) == 0)){ //状态转移 for (int k2 = 0; k2 < mx; k2++){ dp2[i][k1 | (1 << j) | (1 << (j - 1))][k2] = (dp2[i][k1 | (1 << j) | (1 << (j - 1))][k2] + dp2[i][k1][k2]) % mod; } } if (i < n) //看是否可以纵向扩展 //当前的第i+1行的占用情况 for (int k2 = 0; k2 < mx; k2++){ //如果第i+1行的第j个位置(从右向左)为0,则可以纵向扩展 if (((1 << (j-1)) & k2) == 0) //状态转移 dp2[i][k1 | (1 << (j - 1))][k2 | (1 << (j - 1))] = (dp2[i][k1 | (1 << (j - 1))][k2 | (1 << (j - 1))] + dp2[i][k1][k2]) % mod; } } } } //当第i行全部被占用满时,可以将第i+1行的状态 转移到 i+1行的当前行 for (int k = 0; k < mx; k++){ dp2[i + 1][k][0] = dp2[i][mx - 1][k]; } } //最后的结果是n行全部占满,此时第n行为 mx-1, 第n+1行为0 printf("%d\n", dp2[n][mx-1][0]); return 0; }
以上是关于hiho_1048_状态压缩2的主要内容,如果未能解决你的问题,请参考以下文章