Poj -3254 Corn Fields（状压DP）

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

题意：n*m的农田，有些格子可以放牧用1表示，不能放牧的用0表示。要求不能时相邻的格子都有牛。

用二进制来表示状态：假设存在某种情况 1001，1,4可以放牛，2,3不可以。该种情况化为10进制就是5.

对于一个列数为m，可能出现的方案为11.1...1（m个1）。即：++....结果为：-1；而对于这-1中情况，1于1相邻的要去掉（代码注释有）.

对于样列第一行：1 1 1：可能出现的结果有0 0 0，0 0 1，0 1 0，0 1 1.....

对于样列第二行 0 1 1：可能出现的结果为有0 0 0，0 1 0；

DP的思路在于推导：在推导之前要确保：在矩阵中为0的格子是不能放牧的，只有为1的格子才能变为0或1.所以对于样列第二行011，像111，100，110等这种情况是不存在的，所以要去掉（代码注释有）.

推导：假设第i行的情况数为x，第i+1行的情况数为y。当第i+1行出现第1种情况时：第i行x种情况中满足的情况数，当第当第i+1行出现第2种情况时:第i行x种情况中满足的情况数，...当第i+1行出现第y中情况时:第i行x种情况中满足的情况数.（上下的1也不能相邻）

即：dp[i][state(j)]=dp[i-1][state(k1)]+dp[i-1][state(k2)]+......+dp[i-1][state(kn)]。dp[i][state(j)].表示：第i+1行出现第state（j）情况。dp[i-1][state(k)]表示：第i-1行，第state（k）种情况（满足条件的）.

#include<bits/stdc++.h>
using namespace std;
int dp[20][4096];///20表示20行，这里n最大只有12，开大一些防出错。4096=1<<12（最多的情况数）。
int cur[20];///每一行的状态，需要该数组与state[]比较，排除格子为0发生变化的情况。
int state[4096],top;///满足1与1不相邻的情况。
int n,m,num;
const int mod=100000000;

void init()
{
    top=0;
    int total=(1<<m);

    for(int i=0;i<total;i++)///找出1与1不相邻的情况，只要移位就行。如若有1与相邻，肯定会有1&1，则数!=0；
        if(!(i&i<<1)) state[top++]=i;
}
bool flag(int i,int j)
{
    if(i&cur[j]) return 0;
    return 1;
}
int main()
{
    cin>>n>>m;
    init();
    memset(dp,0,sizeof(dp));



    /*for(int i=0;i<top;i++)
        cout<<state[i]<<" ";*/
    for(int i=1;i<=n;i++)
    {
        cur[i]=0;
        for(int j=1;j<=m;j++)
        {
            cin>>num;
            if(num==0) cur[i]+=(1<<(m-j));
            /**
            flag()函数是从已满足1与1不相邻的情况中找出与原情况不矛盾的列子.
            拿样列：0 1 0.按照程序计算后，cur[2]=5,化为二进制就是1 0 1.结合函数flag(),
            如果出现0发生变化的情况,比如:110=6，6&5>0.
            */
        }
    }
    //for(int i=1;i<=n;i++) cout<<cur[i]<<" ";

    for(int i=0;i<top;i++)
        if(flag(state[i],1)) dp[1][i]=1;

    for(int i=2;i<=n;i++)
    {
        ///从第2行开始
        for(int j=0;j<top;j++)///枚举第i行所有的情况
        {
            if(!flag(state[j],i)) continue;
            for(int k=0;k<top;k++)///枚举第i-1行所有的情况
            {
                if(!flag(state[k],i-1)) continue;
                if((state[j]&state[k])) continue;
                ///若满足上下1与1不相邻的情况
                dp[i][j]=(dp[i][j]+dp[i-1][k])%mod;
               /// cout<<dp[i][j]<<" ";
            }
        }
    }


    int ans=0;
    for(int i=0;i<top;i++)
        ans=(ans+dp[n][i])%mod;
    cout<<ans<<endl;
    return 0;
}