poj 2411 状态压缩dp

思路：将每一行看做一个二进制位，那么所有的合法状态为相邻为1的个数一定要为偶数个。这样就可以先把所有的合法状态找到。由于没一层的合法状态都是一样的，那么可以用一个数组保存。由第i-1行到第i行的状态转移是dp[i][now|num[j]]+=dp[i-1][k]，其中now为(1<<m)-1-k;也就是把k中含有0的变1,1边0。k为第i-1行的所有二进制状态，转移条件是k&num[j]==num[j]。唯一注意的是，最后一行的条件是k^num[j]==0.

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define Maxn 13
#define inf 0x7fffffff
using namespace std;
__int64 dp[Maxn][<<Maxn];
int num[<<Maxn],cnt1,cnt2,graphic[Maxn],co,n,m;
void dfs(int j,int f)
{
    int i;
    if(j==m)
    {
        int sum=;
        if(f)
            graphic[j]=;
        else
            graphic[j]=;
        for(i=m;i>=;i--)
            sum+=graphic[i]*(<<(m-i));
        num[++cnt2]=sum;
        return ;
    }
    if(!f)
    {
        graphic[j]=;
        dfs(j+,);
        graphic[j]=;
        dfs(j+,);
    }
    else
    {
        graphic[j]=;
        dfs(j+,);
    }
}
int main()
{
    int i,j,k;
    while(scanf("%d%d",&n,&m)!=EOF,n||m)
    {
        if((n*m)%)
        {
            printf("0\n");
            continue;
        }
        if(n==)
        {
            printf("1\n");
            continue;
        }
        memset(dp,,sizeof(dp));
        graphic[]=;
        cnt2=;
        dfs(,);
        for(i=;i<=cnt2;i++)
            dp[][num[i]]=;
        int temp=<<m;
        temp--;
        for(i=;i<=n-;i++)
        {
            for(j=;j<=cnt2;j++)
            {
                for(k=;k<=temp;k++)
                {
                    if((k&num[j])==num[j])
                    {
                        int now=temp-k;
                        dp[i][now|num[j]]+=dp[i-][k];
                    }
                }
            }
        }
        __int64 ans=;
        for(j=;j<=cnt2;j++)
        {
            for(k=;k<=temp;k++)
            {
                if((k^num[j])==)
                    ans+=dp[i-][k];
            }
        }
        printf("%I64d\n",ans);
    }
    return ;
}