P1879 [USACO06NOV]玉米田Corn Fields-动态规划,状压dp

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

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

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

https://www.luogu.org/problemnew/show/P1879

用F[i]来表示每一行读入时的状态,也就是每行可不可种的状态,F中存二进制数。

for(int i=1;i<=m;i++)
    for(int j=1;j<=n;j++)
        F[i]=(F[i]<<1)+field[i][j];

maxn表示一行的状态最多有2^n个,但其中很多状态都不可用,即不满足不选相邻两块种草。

int maxn=1<<n;

flag[i]表示i这个状态是否可用,即满足不选相邻两块种草。

判断方法就是把这个二进制数左移一位取&,然后右移一位取&。如果这个状态是合法的,那么都应该返回0。简单来说就是这个二进制数向左错开没有相同的,向右错开也没有相同的,即没有相邻的,这个可以自己试一试。比如x=1001,(x<<1)=10010,(x>>1)=100,x&(x<<1)=0,x&(x>>1)=0,这个x就是可行的。

flag[i]=((i&(i<<1))==0)&&((i&(i>>1))==0);

从第一行开始,每一行找合法的情况,(j & F[i]) == j说明没有草地种在贫瘠的地方,这个自己体会一下。

if(flag[j]&&((j&F[i])==j))

那么接下来开始找上一行的合法情况(上下两行之间没有相邻的草地),如果(j&k)==0,就把上一行的情况数加到f[i][j]里

if((k&j)==0)  f[i][j]=(f[i][j]+f[i-1][k])%mo;

最后把最下面一行的每一列的情况书统统加起来,就是答案啦~

for(int i=0;i<maxn;i++){
    ans+=f[m][i];
    ans%=mo;
}

完整代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int mo=1e9;
int F[13],f[13][4096],field[13][13];
bool flag[4096];
int n,m,ans;
int main()
{
    scanf("%d%d",&m,&n);
    for(int i=1;i<=m;i++){
        for(int j=1;j<=n;j++){
            scanf("%d",&field[i][j]);
        }
    }
    for(int i=1;i<=m;i++){
        for(int j=1;j<=n;j++){
            F[i]=(F[i]<<1)+field[i][j];
        }
    }
    int maxn=1<<n;
    for(int i=0;i<maxn;i++){
        flag[i]=((i&(i<<1))==0)&&((i&(i>>1))==0);
    }
    f[0][0]=1;
    for(int i=1;i<=m;i++){
        for(int j=0;j<maxn;j++){
            if(flag[j]&&((j&F[i])==j)){
                for(int k=0;k<maxn;k++){
                    if((k&j)==0){
                        f[i][j]=(f[i][j]+f[i-1][k])%mo;
                    }
                }
            }
        }
    }
    for(int i=0;i<maxn;i++){
        ans+=f[m][i];
        ans%=mo;
    }
    printf("%d",ans);
    return 0;
} 

 

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值