HDU 4049 Tourism Planning (状压DP)

HDU 4049 Tourism Planning (状压DP)

题目链接:

Tourism Planning

题目大意:

有$n (1\le n \le 10)$个人要去$m (1 \le m \le 10)$个景点旅游, 每个景点有一个门票价格$P[i]$, 每一个人对于一个景点有一个喜悦值$V[i][j]$, 对于任意两个人$i, j$, 如果访问同一个景点, 答案加上$B[i][j]$.

访问景点的顺序是固定的, 从1到m, 但是每个人都可以决定他在什么时候离开(包括一个景点都不去).

求出如何规划每个人访问要访问的景点使得最终的$V + B - P​$的和尽可能大.

解题思路:

10个人每个人有10种选法, 暴力$10^{10}$ 显然是不可能的, 这道题用状压DP可以解决.

$dp[i][j]$ 表示到第$i$个景点, 人的状态是$j$的答案最大值. 因为只有10个人, 那么复杂度就是$O(10 * 2 ^{10})$.

转移方程就是从第$i-1$的某一个状态转移到第$i$个景点的某一个状态.

为了做到转移可以O(1)转移, 需要提前进行预处理.

设$value[i][j]$表示第i个景点人的状态是j的获得的答案. 预处理的复杂度也是$O(10 * 2 ^{10})$

DP转移的条件, $dp[i-1][k] - > dp[i][j]$ 的条件是保证k的状态是j状态的一个子集.

/********************************************
 *Author*        :ZZZZone
 *Created Time*  : 三  7/18 21:29:23 2018
 * Ended  Time*  : 三  7/18 22:03:27 2018
*********************************************/

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <stack>
using namespace std;
#define debug(x) std::cerr << #x << " = " << (x) << std::endl
typedef pair<int, int> PII;
typedef long long LL;
typedef unsigned long long ULL;

const int INF = 1 << 30;

int n, m;
int dp[15][1 << 11], value[15][1 << 11];
int P[15], V[15][15], B[15][15];

void Init(){
    for(int i = 1; i <= m; i++) scanf("%d", &P[i]);
    for(int i = 1; i <= n; i++){
        for(int j = 1; j <= m; j++){
            scanf("%d", &V[i][j]);
            V[i][j] -= P[j];
        }
    }
    for(int i = 1; i <= n; i++){
        for(int j = 1; j <= n; j++){
            scanf("%d", &B[i][j]);
        }
    }
    for(int i = 1; i <= m; i++){
        for(int j = 0; j < (1 << n); j++){
            value[i][j] = 0;
            for(int k = 0; k < n; k++){
                if(j & (1 << k)) value[i][j] += V[k + 1][i];
            }
            for(int k = 0; k < n; k++){
                for(int l = k + 1; l < n; l++){
                    if((j & (1 << k)) && (j & (1 << l)))
                        value[i][j] += B[k+1][l+1];
                }
            }
        }
    }
}

void Solve(){
    for(int i = 0; i <= m; i++)
        for(int j = 0; j < (1 << n); j++)
            dp[i][j] = -INF;
    dp[0][(1 << n)-1] = 0;
    for(int i = 1; i <= m; i++){
        for(int j = 0; j < (1 << n); j++){
            for(int k = 0; k < (1 << n); k++){
                if((j & k) == j) dp[i][j] = max(dp[i][j], dp[i-1][k] + value[i][j]);
            }
        }
    }
    int ans = 0;
    for(int i = 0; i < (1 << n); i++) ans = max(ans, dp[m][i]);
    if(ans <= 0) printf("STAY HOME\n");
    else printf("%d\n", ans);
}

int main()
{
    while(~scanf("%d %d", &n, &m)){
        if(n == 0 && m == 0) break;
        Init();
        Solve();
    }
    return 0;
}