POJ 3254 Corn Fields (状压DP)

状压DP第一发!

http://www.tuicool.com/articles/JVzMVj讲解较好的博客。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<map>
#include<queue>
#include<cmath>
#include<algorithm>
#define MAX 1000
#define inf 0x3f3f3f3f
using namespace std;
const int mod=100000000;
int dp[20][5000],state[5000],m,n,top,cur[20];
int look(int t){
   if(t&t<<1) return false;///&的优先级小于<<
   else return true;
}

void init(){
    top=0;
    int tmp=1<<n;
    for(int i=0;i<tmp;++i)///所有的可行状态都用一个10进制数来表示
       if(look(i))
            state[++top]=i;///存放当前一行的所有状态
}

bool fi(int x,int t){
    if(x&cur[t]) return false;
    else  return true;
}

int main(){
    int i,j,k;
    while(~scanf("%d%d",&m,&n)){
       // printf("%d\n",1<<0);
        init();
        memset(dp,0,sizeof(dp));
        memset(cur,0,sizeof(cur));
        for(i=1;i<=m;i++){
            cur[i]=0;
            int x;
            for(j=1;j<=n;j++){
                scanf("%d",&x);
                if(x==0)
                    cur[i]+=(1<<(n-j));///第i行整行的情况,为了方便的描述出1(实际上是0所表示的值)如 1 1 1
            }///0的位置存1因为,如果当牧场不放牛的时候,也算是一种解法 0 1 0 那么第二行的第一个0就是5最后一个是1
        }


        for(i=1;i<=top;i++){
            if(fi(state[i],1))
                dp[1][i]=1;
        }

        for(i=2;i<=m;i++){
            for(int k=1;k<=top;k++){
                if(!fi(state[k],i))continue;
                for(j=1;j<=top;j++){
                    if(!fi(state[j],i-1))continue;
                    if(state[k]&state[j])continue;
                    dp[i][k]=(dp[i][k]+dp[i-1][j])%mod;
                }
            }
        }
        int ans=0;
        for(i=1;i<=top;i++)
            ans=(ans+dp[m][i])%mod;
        printf("%d\n",ans);
    }
    return 0;
}

你可能感兴趣的:(poj,状压dp)