P1879 [USACO06NOV]玉米田Corn Fields（状压dp）

P1879 [USACO06NOV]玉米田Corn Fields

状压dp水题

看到$n,m<=12$，肯定是状压鸭

先筛去所有不合法状态，蓝后用可行的状态跑一次dp就ok了

#include<iostream>
#include<cstdio>
#include<cstring>
#define re register
using namespace std;
const int p=1e9;
int n,m,a[13],f[13][4098];
int t[13][4098],len[13];
int mod(int a){return a<p?a:a-p;}
void prepare(){
    for(re int i=1;i<=n;++i)
        for(re int j=0;j<(1<<m);++j){
            if((j&(j<<1))||(j&(j>>1))||(j&a[i])) continue;
            t[i][++len[i]]=j;//合法的存起来
        }
}
int main(){
    scanf("%d%d",&n,&m); int q;
    for(re int i=1;i<=n;++i)
        for(re int j=0;j<m;++j)
            scanf("%d",&q),a[i]=(a[i]<<1)+(q^1);//转存二进制并取反（有利于后面的筛）
    prepare();
    for(re int i=1;i<=len[1];++i) f[1][t[1][i]]=1;
    for(re int i=2;i<=n;++i)
        for(re int j=1;j<=len[i-1];++j)
            for(re int k=1;k<=len[i];++k){
                if(t[i-1][j]&t[i][k]) continue;//不合法
                f[i][t[i][k]]=mod(f[i][t[i][k]]+f[i-1][t[i-1][j]]);
            }
    int ans=0;
    for(re int i=1;i<=len[n];++i) ans=mod(ans+f[n][t[n][i]]);
    printf("%d",ans);
    return 0;
}