状压dp入门第一题 poj3254

题目链接 http://poj.org/problem?id=3254

转自http://blog.csdn.net/harrypoirot/article/details/23163485

 #include <stdio.h>

 #include <math.h>

 #include <string.h>

 #include <stdlib.h>

 #include <iostream>

 #include <sstream>

 #include <algorithm>

 #include <string>

 #include <queue>

 #include <vector>

 using namespace std;

 const int maxn= 1e5 + ;

 const int inf = 0x3f3f3f3f;

 typedef long long ll;

 #define mod 100000000

 int M,N,top = ;

 //top表示每行最多的状态数

 int state[],num[];

 //state存放每行所有的可行状态（即没有相邻的状态

 //

 int dp[][];

 //dp[i][j]:对于前i行数据，每行有前j种可能状态时的解

 int cur[];

 //cur[i]表示的是第i行整行的情况

 inline bool ok(int x){    //判断状态x是否可行

    if(x&x<<)    return false;//若存在相邻两个格子都为1，则该状态不可行

    return true;

 }

 void init(){            //遍历所有可能的状态

    top = ;

    int total =  << N; //遍历状态的上界

    for(int i = ; i < total; ++i){

        if(ok(i))state[++top] = i;

    }

 }

 inline bool fit(int x,int k){ //判断状态x 与第k行的实际状态的逆是否有‘重合’

    if(x&cur[k])return false; //若有重合，（即x不符合要求）

    return true;  //若没有，则可行

 }

 int main(){

     while(scanf("%d%d",&M,&N)!= EOF){

        init();

        memset(dp,,sizeof(dp));

        for(int i = ; i <= M; ++i){

            cur[i] = ;

            int num;

            for(int j = ; j <= N; ++j){  //输入时就要按位来存储，cur[i]表示的是第i行整行的情况，每次改变该数字的二进制表示的一位

                 scanf("%d",&num);  //表示第i行第j列的情况（0或1）

                if(num == ) //若该格为0

                    cur[i] +=(<<(N-j)); //则将该位置为1（注意要以相反方式存储，即1表示不可放牧

            }

        }

        for(int i = ;i <= top;i++){

            if(fit(state[i],)){  //判断所有可能状态与第一行的实际状态的逆是否有重合

                 dp[][i] = ;  //若第1行的状态与第i种可行状态吻合，则dp[1][i]记为1

            }

        }

        /*

        状态转移过程中，dp[i][k] =Sigma dp[i-1][j] （j为符合条件的所有状态）

         */

        for(int i = ; i <= M; ++i){  //i索引第2行到第M行

            for(int k = ; k <= top; ++k){ //该循环针对所有可能的状态，找出一组与第i行相符的state[k]

                 if(!fit(state[k],i))continue; //判断是否符合第i行实际情况

                 for(int j = ; j <= top ;++j){ //找到state[k]后，再找一组与第i-1行符合，且与第i行（state[])不冲突的状态state[j]

                    if(!fit(state[j],i-))continue;  //判断是否符合第i-1行实际情况

                    if(state[k]&state[j])continue;  //判断是否与第i行冲突

                    dp[i][k] = (dp[i][k] +dp[i-][j])%mod;  //若以上皆可通过，则将'j'累加到‘k'上

                 }

            }

        }

        int ans = ;

        for(int i = ; i <= top; ++i){ //累加最后一行所有可能状态的值，即得最终结果！！！泥马写注释累死我了终于写完了！

            ans = (ans + dp[M][i])%mod;

        }

        printf("%d\n",ans);

    }

 }

状压dp入门第一题 poj3254的相关教程结束。