poj 3254 Corn Fields(动态规划:状压DP)

来源:互联网 发布:剑灵茱莉亚捏脸数据 编辑:程序博客网 时间:2024/06/03 22:52

题意是给出一块草地,分为m*n格

某个格子为1代表可以养牛,为0代表不可以养牛

相邻的草地不可以同时养牛

问有多少种放牛的方法?

看着别人的解题报告写的

因为相邻草地不可以同时放牛,所以我们保存可以放牛对应的状态

竖着相邻的草地也不可以同时放牛,所以如果同一列相邻两行&为真,不可以放牛

同时对于当前为0的草地也不可以放牛

把这三种情况剔除,就可以得到结果

状态转移方程为:dp[i][j] = sum(dp[i-1][k])//dp[i][j]代表当前第i行,按照j的二进制放牛对应结果数

这里输入的时候把0 1转换了下,看了很久才看明白

因为我们要保证当前的状态同时满足草地上对应值为1

所以反转后判断比较容易

代码如下:

#include <cstdio>#include <cstring>#include <algorithm>#define MAXN 10010#define MOD 100000000using namespace std;int cur[13];int dp[13][1<<13];int state[1<<13];bool ok(int i, int j) {    if(i & j) return 0;    return 1;}int main(void) {    int m, n, tmp;    while(scanf("%d%d", &m, &n) != EOF) {        memset(dp, 0, sizeof(dp));        int top = 0;        for(int i=0; i<(1<<n); ++i) {            if(ok(i, i<<1))                state[++top] = i;        }        for(int i=1; i<=m; ++i) {            cur[i] = 0;            for(int j=1; j<=n; ++j) {                scanf("%d", &tmp);                if(tmp == 0)                    cur[i] += 1<<(n-j);            }        }        for(int i=1; i<=top; ++i) {            if(ok(state[i], cur[1])) {                dp[1][i] = 1;                //printf("dp[1][%d] = %d\n", i, dp[1][i]);            }        }        for(int i=2; i<=m; ++i) {            for(int j=1; j<=top; ++j) {                if(!ok(state[j], cur[i]))                    continue;                for(int k=1; k<=top; ++k) {                    if(ok(state[k], cur[i-1]) && ok(state[k], state[j]))                        dp[i][j] = (dp[i][j]+dp[i-1][k])%MOD;                }            }        }        int ans = 0;        for(int i=1; i<=top; ++i)            ans = (ans+dp[m][i])%MOD;        printf("%d\n", ans);    }    return 0;}


0 0
原创粉丝点击