POJ 3254 状压DP

最新推荐文章于 2019-09-07 11:35:59 发布

原创最新推荐文章于 2019-09-07 11:35:59 发布 · 588 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#dp #状压

dp 专栏收录该内容

18 篇文章

订阅专栏

本文解析了一道经典的状压动态规划题目，该题目要求在特定的地图上计算出所有合法的放牧方案数量。文章提供了完整的C++代码实现，并详细介绍了关键的数据结构和算法思路。

题意：给你一片地，1代表着有粮食，0代表没有，现在要进行放牧，只有1的地方才能放牧（要不没食物），然后两个放牧的地方不能相邻，问有多少种放牧的方法。

思路：算是状压DP模板题吧，复习一下。

代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = ( 1 << 12 ) + 100;
const int mod = 100000000;
int n,m;
int mm[15];
bool can[maxn];
int dp[14][maxn];
void init(){
    for(int i = 0;i < (1 << 12);i ++){
        if(!(i&(i>>1)))can[i] = true;
    }
}
inline bool ok(int a,int b){
    return (!(a&b));
}
inline bool fit(int a,int b){
    return !((a|b) - b);//这地方改成这样反而没有原来一位一位判快了。。估计是哪里写挫了
}
int main(){
    init();
    while(~scanf("%d%d",&n,&m)){
        memset(mm,0,sizeof(mm));
        for(int i = 1;i <= n;i ++){
            for(int j = 0;j < m;j ++){
                int a;
                scanf("%d",&a);
                mm[i] = (mm[i] << 1) + a;
            }
        }//for(int i=1;i<=n;i++)cout<<mm[i]<<' ';cout<<endl;
        memset(dp,0,sizeof(dp));
        for(int i = 0 ;i < (1 << n);i ++)
            if(can[i])dp[0][i] = 1;
        for(int i = 1;i <= n;i ++){
            for(int j = 0;j < (1 << m);j ++)
                if(can[j]&&fit(j,mm[i-1]))
                for(int k = 0;k < (1 << m);k ++){
                    if(can[k]&&ok(j,k)&&fit(k,mm[i]))
                        dp[i][k] = ( dp[i][k] + dp[i-1][j] )%mod;
                }
        }
        int ans = 0;
        for(int i = 0;i < (1 << m);i ++)
            ans += dp[n][i];
        cout<<ans%mod<<endl;
    }
    return 0;
}