[状态压缩DP] PKU 2411 Mondriaan's Dream

参考http://wenku.baidu.com/view/e262a86f1eb91a37f1115c26.html

状态设计：0表示横放，1表示竖放，可以预处理出所有合法状态，f[i, s] 表示到达第 i 层的状态（类似俄罗斯方块，前 i 层合并之后剩余的状态）。

状态转移：f[i, s] = sum{f[i-1, s']},s' 为和 s 符合的状态 j 与 s 的异或值。

# include <cstdio>
# include <cstring>

# define N 11 + 1

typedef long long int LL;

int m, n;
LL f[N][1<<N];

bool check(int x)
{
    int tmp = 0;
    for (int i = 0; i < n; ++i)
    {
        if ((x>>i)&0x1) continue;
        else
        {
            while (i<n && (((~x)>>i)&0x1))
            {
                ++i;
                ++tmp;
            }
            if (tmp & 0x1) return false;
            else tmp = 0;
        }
    }
    //printf("%04X\n", x);
    return true;
}

LL dp(int i, int s)
{
        //printf(" : %d\n", i);
    LL &ans = f[i][s];
    if (ans != -1) return ans;
    ans = 0;
    for (int j = 0; j < (1<<n); ++j)
    {
        if ((j^s)==j-s && check(j)) ans += dp(i-1, j^s);
    }
    return ans;
}

void solve(void)
{
    memset(f[0], 0, sizeof(LL)*(1<<n));
    f[0][0] = 1;
    for (int i = 1; i <= m; ++i)
        memset(f[i], -1, sizeof(LL)*(1<<n));
    LL ans = dp(m, 0);                        // 0 ~ n-1
    printf("%lld\n", ans);
}

int main()
{
    while (scanf("%d%d", &m, &n), m||n)
        solve();

    return 0;
}