BZOJ 1725: [Usaco2006 Nov]Corn Fields牧场的安排 状压动归
Code:
#include<bits/stdc++.h>
#define setIO(s) freopen(s".in","r",stdin)
#define mod 100000000
#define maxn 20
#define ll long long
using namespace std;
ll F[15][20000],ans;
int G[20][20],w[20];
vector<int>sta[maxn],h;
int main(){
// setIO("input");
int n,m;
scanf("%d%d",&n,&m);
for(int i=1;i<15;++i) w[i]=(1<<(i-1));
//预处理出所有合法状态
for(int k=0;k<w[m+1];++k)
{
int l=k,pre=0,flag=1;
while(l)
{
if((l&1) && pre) { flag=0; break; }
if(l&1)
pre = 1;
else
pre = 0;
l>>=1;
}
if(flag) h.push_back(k);
}
for(int i=1;i<=n;++i)
{
int cur=0;
for(int j=1;j<=m;++j)
{
scanf("%d",&G[i][j]);
if(G[i][j]==0) cur|=w[j]; //不合法
}
for(int j=0,sz=h.size();j<sz;++j)
{
if(h[j]&cur) continue;
sta[i].push_back(h[j]);
}
}
F[1][0]=1;
for(int i=0,sz=sta[1].size();i<sz;++i) F[1][sta[1][i]]=1;
for(int i=2;i<=n;++i)
{
for(int j=0,sz=sta[i].size();j<sz;++j)
{
for(int k=0,sz2=sta[i-1].size();k<sz2;++k){
if(sta[i][j]&sta[i-1][k]) continue;
F[i][sta[i][j]]=(F[i][sta[i][j]]+F[i-1][sta[i-1][k]])%mod;
}
}
}
for(int i=0,sz=sta[n].size();i<sz;++i) ans=(ans+F[n][sta[n][i]])%mod;
printf("%lld",ans);
return 0;
}