POJ3254 状压DP

转载：http://blog.csdn.net/accry/article/details/6607703

初学状压就找个比较简单的dp，应该是最好想的状态压缩了。

状态压缩dp也就是如果正常开数组，必然会导致内存过大，所有需要另外找一种办法来代替这个问题，那么就用二进制来替代，这应该就是对于状态压缩最通俗的解释了。关于一些二进制的操作比如

1 如何判断数字x第i位是否为1    1<<(i-1) & x

2 将一个数字x二进制下第i位更改成1  x= x|(1<<(i-1))

3 把一个数字二进制下最靠右的第一个1去掉  x=x&(x-1)  (BIT里面我记得用到的就是这个）

总之，二进制真的是一个十分巧妙的东西，也是一个非常好用的东西

然后是经典的TSP（旅行商问题）

代码什么的白书上都有，然后白书上一些对于位运算的操作

http://blog.csdn.net/qq_32400847/article/details/51814131这篇文章都有很详细的解释就不写了。

poj3254

题目大意是：给一个n*m的矩阵，上面1表示肥沃的土地，0表示贫瘠的土地，然后要让牛必须站在肥沃的土地上，又不能相邻（可以整个都不站牛）

题解：那就用二进制来表示整个土地的样子，状态转移是dp[i][j] = sum(dp[i-1][0-cnt]) 表示第i行的站法是由i-1行所有可能的站法并且两行不能冲突的和

有一个技巧就是如果判断一个数字是否有相邻的两个1  x&（x<<1） 如果是0则表示没有1则表示有。

#include <stdio.h>
#include <algorithm>
#include <stdlib.h>

using namespace std;
typedef long long ll;
const ll MOD = 100000000;

int sz ;
ll dp[14][400];
int cur[14];
int mark[400];
bool ok(int i)
{
    return !(i & (i << 1));
}
void add(ll &x, ll y)
{
    x += y;
    if(x >= MOD) x -= MOD;
}
bool judge(int x, int y)
{
    return x == (x | y);
}
int main()
{
    int n, m;
    scanf("%d%d",&n,&m);
    for(int i = 1;i <= n;i ++) {
        for(int j = 1;j <= m;j ++) {
            int tmp;
            scanf("%d",&tmp);
            if(tmp) cur[i] += 1<<(j-1);
        }
    }
//    for(int i = 1;i <= n;i ++) printf("%d\n",cur[i]);

    for(int i = 0;i < 1 << m; i ++) {
        if(ok(i)) {
            mark[++sz] = i;
        }
    }
    for(int i = 1;i <= sz;i ++) {
        if(judge(cur[1], mark[i])) dp[1][i] = 1;
    }
    for(int i = 2;i <= n;i ++) {
        for(int now = 1;now <= sz;now ++) {
            if(!judge(cur[i], mark[now])) continue;
            for(int last = 1;last <= sz;last ++) {
                if(!judge(cur[i-1], mark[last])) continue;
                if(!(mark[last] & mark[now])) {
                    add(dp[i][now], dp[i-1][last]);
                }
            }
        }
    }
    ll res = 0;
    for(int i = 1;i <= sz;i ++) add(res, dp[n][i]);
    printf("%lld\n",res);
    return 0;
}