bzoj1725 [Usaco2006 Nov]Corn Fields牧场的安排

状压dp,预处理出:
line[i]:第i行状态是否合法
ok[i][j],表示第i行的状态到第j行的状态是否能转移。
situ[i][j],表示i状态转到j状态是否可行。
设f[i][j]表示第i行状态为j时的答案,那么我们枚举一个k,对于所有合法的,能够转移到j的k,加上去。最后答案就是sigma(f[n][i])。

#include
#include
#include
#include
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
int n,m,k;
const int N=1<<12;
const int mo=1e8;
bool situ[13][1<<12],line[1<<12],ok[1<<12][1<<12];
int f[14][1<<12];
int tot,state[13],mi[13],s[13];
int main()
{
    scanf("%d%d",&n,&m);
    mi[0]=1;
    tot=(1<<m)-1;
    fo(i,1,13)mi[i]=(mi[i-1]*2)%mo;
    fo(i,1,n)
    {
        int x;
        fd(j,m,1)scanf("%d",&s[j]);
        x=0;
        fo(j,1,m)x+=mi[j-1]*s[j];
        state[i]=x;
    }
    fo(i,0,tot)
    if (((i>>1)&i)==0)line[i]=1;
    fo(i,0,tot)if(line[i])
    fo(j,0,tot)if(line[j])
    if ((i&j)==0)ok[i][j]=1;
    fo(i,1,n)
    fo(j,0,tot)
    if ((state[i]|j)==state[i])situ[i][j]=1;
    fo(i,0,tot)
    if (situ[1][i]&&line[i])f[1][i]=1;
    fo(i,2,n)
    fo(j,0,tot)
    if (line[j]&&situ[i][j])
    fo(k,0,tot)
    if (line[k]&&situ[i-1][k])
    if (ok[j][k])
    f[i][j]=(f[i][j]+f[i-1][k])%mo;
    int ans=0;
    fo(i,0,tot)ans=(ans+f[n][i])%mo;
    printf("%d\n",ans);
    return 0;
}

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