JDOJ1098: 夏娜的菠萝包——状压DP

一.题目

Description
问题描述：夏娜很喜欢吃菠萝包，她的经纪人RC每半个月就要为她安排接下来的菠萝包计划。今天是7月份，RC又要去商场进货买菠萝包了。这次RC总共买了N种菠萝包，每种一个。每个菠萝包都有一个初始美味值Ti,每过一天就会减少Di，即第2天美味值为Ti-Di,第3天为Ti-2*Di，依此类推。一旦美味值减为负数，那个包就坏掉了，不能吃了。 RC每天都要为夏娜安排当天吃菠萝包的组合，这些组合不是随意的，而是只能从夏娜喜欢的M种搭配中挑选一种。每种搭配是由Ki个菠萝包组成的，一种搭配的总美味值是这Ki个菠萝包当天的美味值之和再加上一个额外的搭配美味值Ei。不过要注意，一旦某种搭配的其中一个菠萝包坏掉了，这个搭配就不能选用了。而且，有可能存在两个搭配，里面的组合是一样的，但额外的搭配美味值却不同。 RC想让可爱的夏娜尽可能地吃得美味，因此希望能找出一种最优的方案，让小夏娜吃上若干天的菠萝包，这些天的美味值之和最大。但RC面临着两个邪恶的敌人，一个叫bug，一个叫zzy,他们也想抢夺这个经纪人之位，因此要是他们提出更优的方案，RC就可能会失去他的夏娜了。那么，你们能帮帮这个可怜的RC吗？
数据范围
$n\le 14,m\le 20$
传送门

二.Solution

一眼状压DP，对吧，这数据范围，绝了。
考虑定义状态：

$dp[i][j]$表示第$i$天吃了的菠萝包的状态为$j$的最大美味值

于是有状态转移：
$$dp[i][j]=max(dp[i][j],dp[i-1][j异或type[k]]+w[k][i])$$
w[k][i]是第i天k种搭配的美味值，预处理一下就行了。

三.Code

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

#define M 50
#define INF 0x3f3f3f3f

int n, m, t[M], d[M], k[M], e[M], w[M][M], dp[M][(1 << 14)], type[M], ans;

int main (){
    while (scanf ("%d", &n) && n){
        ans = 0;
        memset (dp, 0, sizeof dp);
        memset (w, 0, sizeof w);
        memset (type, 0, sizeof type);
        for (int i = 1; i <= n; i ++){
            scanf ("%d %d", &t[i], &d[i]);
        }
        scanf ("%d", &m);
        for (int i = 1; i <= m; i ++){
            scanf ("%d %d", &k[i], &e[i]);
            for (int j = 1; j <= k[i]; j ++){
                int id;
                scanf ("%d", &id);
                type[i] ^= (1 << id - 1);
                for (int z = 0; z < n; z ++){
                    w[i][z + 1] += t[id] - d[id] * z;
                    if (t[id] - d[id] * z < 0)
                        w[i][z + 1] = -INF;
                }
            }
        }
        for (int i = 1; i <= m; i ++)
            for (int j = 1; j <= n; j ++)
                w[i][j] += e[i];
        for (int i = 1; i <= n; i ++){
            for (int j = 0; j < (1 << n); j ++){
                for (int z = 1; z <= m; z ++){
                    if ((j & type[z]) == type[z])
                        dp[i][j] = max (dp[i][j], dp[i - 1][j ^ type[z]] + w[z][i]);
                }
                ans = max (ans, dp[i][j]);
            }
        }
        printf ("%d\n", ans);
    }
    return 0;
}

Thanks！