HDU 4049 Tourism Planning (状压DP)
题目链接:
题目大意:
有 n(1≤n≤10) n ( 1 ≤ n ≤ 10 ) 个人要去 m(1≤m≤10) m ( 1 ≤ m ≤ 10 ) 个景点旅游, 每个景点有一个门票价格 P[i] P [ i ] , 每一个人对于一个景点有一个喜悦值 V[i][j] V [ i ] [ j ] , 对于任意两个人 i,j i , j , 如果访问同一个景点, 答案加上 B[i][j] B [ i ] [ j ] .
访问景点的顺序是固定的, 从1到m, 但是每个人都可以决定他在什么时候离开(包括一个景点都不去).
求出如何规划每个人访问要访问的景点使得最终的 V+B−P V + B − P 的和尽可能大.
解题思路:
10个人每个人有10种选法, 暴力 1010 10 10 显然是不可能的, 这道题用状压DP可以解决.
dp[i][j] d p [ i ] [ j ] 表示到第 i i 个景点, 人的状态是的答案最大值. 因为只有10个人, 那么复杂度就是 O(10∗210) O ( 10 ∗ 2 10 ) .
转移方程就是从第 i−1 i − 1 的某一个状态转移到第 i i 个景点的某一个状态.
为了做到转移可以O(1)转移, 需要提前进行预处理.
设表示第i个景点人的状态是j的获得的答案. 预处理的复杂度也是 O(10∗210) O ( 10 ∗ 2 10 )
DP转移的条件, dp[i−1][k]−>dp[i][j] d p [ i − 1 ] [ k ] − > d p [ 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;
}