***HDU 3237 - Help Bubu(状压DP)

题目：

http://acm.hdu.edu.cn/showproblem.php?pid=3237

题意：

给出n本书的高度（25<=h<=32），相等高度的书为一个区间，区间的数量为这堆书的杂乱度。

现在你可以移动至多k本书，求出最小的杂乱度。

思路：

书的高度只有 0~7，所以可以使用状压来表示书的高度的状态。如果前i本书中存在高度为h的书，那么s的第h位表示为1.

dp[i][j][s][h]: 表示前i本书取走j本，书的高度状态为s，且最后一本书的高度为h的最小杂乱度。

对于第i为的书：

1.取走： dp[i&1][j+1][s][h] = min(dp[i&1][j+1][s][h], dp[(i-1)&1][j][s][h])

2.不取走：

    1）最后一本书高度等于hi：dp[i&1][j][s][h] = min(dp[i&1][j][s][h], dp[(i-1)&1][j][s][h])

    2）不等于： dp[i&1][j][s|(1<<val[i])][val[i]] = min(dp[i&1][j][s|(1<<val[i])][val[i]],  dp[(i-1)&1][j][s][h]+1)

将书架上没有书的状态表示为 高度h=8.

AC.

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;
const int inf = 0x3f3f3f3f;
const int maxn = 105;
int num[300], val[maxn];
int dp[2][maxn][300][10];

int main()
{
    //freopen("in", "r", stdin);
    int n, m, ca = 1;
    memset(num, 0, sizeof(num));
    for(int i = 0; i < (1<<8); ++i) {
        for(int j = 0; j <= 8; ++j) {
            if(i & (1<<j)) num[i]++;
        }
    }

    while(~scanf("%d%d", &n, &m)) {
        if(n == 0 && m == 0) break;
        int inil = 0;
        for(int i = 1; i <= n; ++i) {
            scanf("%d", &val[i]);
            val[i] -= 25;
            inil = inil | (1<<val[i]);
        }

        for(int j = 0; j <= m; ++j) {
            for(int s = 0; s < (1<<8); ++s) {
                for(int h = 0; h <= 8; ++h) {
                    dp[0][j][s][h] = inf;
                }
            }
        }

        dp[0][0][0][8] = 0;
        for(int i = 1; i <= n; ++i) {

            for(int j = 0; j <= m; ++j) {
                for(int s = 0; s < (1<<8); ++s) {
                    for(int h = 0; h <= 8; ++h) {
                        dp[i&1][j][s][h] = inf;
                    }
                }
            }

            for(int j = 0; j < i && j <= m; ++j) {
                for(int s = 0; s < (1<<8); ++s) {
                    for(int h = 0; h <= 8; ++h) {
                        if(dp[(i-1)&1][j][s][h] == inf) continue;
                        //取走第i本
                        if(j < m)
                            dp[i&1][j+1][s][h] = min(dp[i&1][j+1][s][h], dp[(i-1)&1][j][s][h]);
                        //不取走第i本
                        if(h == val[i])
                            dp[i&1][j][s][h] = min(dp[i&1][j][s][h], dp[(i-1)&1][j][s][h]);
                        else
                            dp[i&1][j][s|(1<<val[i])][val[i]] = min(dp[i&1][j][s|(1<<val[i])][val[i]],
                                                                    dp[(i-1)&1][j][s][h]+1);
                    }
                }
            }
        }

        int ans = inf;
        for(int j = 0; j <= m; ++j) {
            for(int s = 0; s < (1<<8); ++s) {
                for(int h = 0; h < 8; ++h) {
                    if(dp[n&1][j][s][h] == inf) continue;
                    int tmp = inil ^ s;
                    ans = min(ans, dp[n&1][j][s][h] + num[tmp]);
                }
            }
        }

        printf("Case %d: %d\n\n", ca++, ans);
    }
    return 0;
}