UVa 12235 状压DP Help Bubu

题解戳这

一开始没看懂题解，后来想明白以后，d(i, j, s, x)是考虑第i本书的时候，前面已经拿走了j本书，剩下的书的种类的二进制状态为s，剩下的最后一本书的编号为x，所能得到的最小混乱度。

这里状态定义的时候，先不考虑把拿出来的书放回去。

最后统计答案的时候，把那些拿出来的书再加上。

all是所有n本书的状态，s是剩下书的种类的状态。

如果拿出来的书中有和前面高度相同的，直接插到相邻的位置就行了，不会增加混乱度。

如果拿出来的书中没有和前面高度相同的，不管放在那里混乱度都会加1，这样所增加的混乱度就是bitcount(all ^ s)

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

const int maxn = 100 + 5;
const int maxs = (1 << 8);
const int INF = 0x3f3f3f3f;

int n, k;

int a[maxn];
int d[2][maxn][maxs][8];

int bitcount(int x)
{
    int ans = 0;
    while(x) { ans += (x&1); x >>= 1; }
    return ans;
}

int main()
{
    int kase = 0;

    while(scanf("%d%d", &n, &k) == 2)
    {
        if(!n && !k) break;

        for(int i = 0; i < n; i++) scanf("%d", a + i);

        memset(d[0], 0x3f, sizeof(d[0]));
        int cur = 0;
        int all = 0;

        for(int i = 0; i < n; i++)
        {
            cur ^= 1;
            memset(d[cur], 0x3f, sizeof(d[cur]));
            int h = a[i] - 25;
            d[cur][i][1<<h][h] = 1;

            for(int j = 0; j <= min(k, i); j++)
                for(int s = all; s; s = (s-1)&all)
                    for(int x = 0; (1 << x) <= s; x++)
                        if(d[cur^1][j][s][x] != INF)
                        {
                            int t = d[cur^1][j][s][x];
                            if(x == h) d[cur][j][s][x] = min(d[cur][j][s][x], t);
                            else
                            {
                                d[cur][j+1][s][x] = min(d[cur][j+1][s][x], t);
                                d[cur][j][s|(1<<h)][h] = min(d[cur][j][s|(1<<h)][h], t + 1);
                            }
                        }
            all |= (1 << h);
        }

        int ans = n;
        for(int j = 0; j <= k; j++)
            for(int s = all; s; s = (s-1)&all)
                for(int x = 0; (1 << x) <= s; x++)
                    ans = min(ans, d[cur][j][s][x] + bitcount(all^s));
        printf("Case %d: %d

", ++kase, ans);
    }

    return 0;
}