1725: [Usaco2006 Nov]Corn Fields牧场的安排

Description

Farmer John新买了一块长方形的牧场，这块牧场被划分成M列N行(1<=M<=12; 1<=N<=12)，每一格都是一块正方形的土地。FJ打算在牧场上的某几格土地里种上美味的草，供他的奶牛们享用。遗憾的是，有些土地相当的贫瘠，不能用来放牧。并且，奶牛们喜欢独占一块草地的感觉，于是FJ不会选择两块相邻的土地，也就是说，没有哪两块草地有公共边。当然，FJ还没有决定在哪些土地上种草。 作为一个好奇的农场主，FJ想知道，如果不考虑草地的总块数，那么，一共有多少种种植方案可供他选择。当然，把新的牧场荒废，不在任何土地上种草，也算一种方案。请你帮FJ算一下这个总方案数。

Input

* 第1行: 两个正整数M和N，用空格隔开

* 第2..M+1行: 每行包含N个用空格隔开的整数，描述了每块土地的状态。输入的第i+1行描述了第i行的土地。所有整数均为0或1，是1的话，表示这块土地足够肥沃，0则表示这块地上不适合种草

Output

* 第1行: 输出一个整数，即牧场分配总方案数除以100,000,000的余数

Sample Input

2 3
1 1 1
0 1 0

Sample Output

9

输出说明:

按下图把各块土地编号：

1 2 3
4

只开辟一块草地的话，有4种方案：选1、2、3、4中的任一块。开辟两块草地的话，有3种方案：13、14以及34。选三块草地只有一种方案：134。再加把牧场荒废的那一种，总方案数为4+3+1+1=9种。

HINT

 

Source

Gold

 

 

状压DP。。。其实就是把每行都看作一个单独的状态。。。然后枚举前一行的局面然后在转移到这一行。。。

#include<iostream>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<string>
#include<map>
#include<queue>
#include<vector>
#include<set>
#define inf 1000000000
#define maxn 10000+5
#define maxm 10000+5
#define mod 100000000
#define eps 1e-10
#define ll long long
#define for0(i,n) for(int i=0;i<=(n);i++)
#define for1(i,n) for(int i=1;i<=(n);i++)
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
#define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)
using namespace std;
int ans,ed;
int t[13],mp[13],f[13][4096];
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
    return x*f;
}
int main(){
    //freopen("input.txt","r",stdin);
    //freopen("output.txt","w",stdout);
    int m=read(),n=read();
    for1(i,m)
        for1(j,n){
            mp[i]<<=1;mp[i]+=read();
        }
    ed=(1<<n)-1;
    for0(i,ed)
        if((i&(i>>1))==0&&(i|mp[1])==mp[1])f[1][i]=1;
    for(int i=2;i<=m;i++)
        for(int j=0;j<=ed;j++){
            if(f[i-1][j])
                for(int k=0;k<=ed;k++){
                    if((j&k)==0&&(k|mp[i])==mp[i]&&(k&(k>>1))==0)
                    f[i][k]=(f[i][k]+f[i-1][j])%mod;
                }
        }
    for(int i=0;i<=ed;i++)
        ans=(ans+f[m][i])%mod;
    printf("%d",ans);
    return 0;
}