时间限制: 1 Sec 内存限制: 128 MB
题目描述
经过几个月的艰苦学习,Iserlohn终于赢得了全额奖学金。作为一个运动鞋狂热爱好者,他决定用所有的钱在运动鞋商店进行消费。
有一些球鞋Iserlohn要收集,如Air Jordan 和 Nike Pro。而每个品牌已发布各种产品。由于,Iserlohn绝对是一个运动鞋狂热,他意欲购买每个品牌至少有一个产品。
虽然每个产品的都有相应的定价,但是在Iserlohn心中也有一个价值倾向。由于他的钱有限,他希望最大限度地提高鞋的总价值。显然,作为一个收藏家,他不会购买相同产品的两次。
现在,Iserlohn需要你帮他找到他的答案,即在至少购买每个品牌的一个产品的情况下,使得所有运动鞋在他心中的总价值最高。
输入
第1行包含三个整数N,M,K(1<=N<=100,1 <= M<= 10000,1<=K<=10)分别表示共有N种产品,M的钱和K种品牌。
接下来N行,每行三个整数a,b,c(1<=a<=k, 0<=b,c<100000)分别表示该商品属于哪种品牌,商品标价以及他心目中的价格。
输出
只有一个正整数,输出他心中的最大价值,如果不能满足他的要求,输出“Impossible”
样例输入
5 10000 3
1 4 6
2 5 7
3 4 99
1 55 77
2 44 66
样例输出
255
Solution
考虑到问题的描述:
K个品牌,花J元买的最大价值,且需要保证每种品牌至少一件产品。
转化成最优子结构就是:
F[i][j][k]表示前i个品牌,买到第i个品牌的第j件物品,花k元能买到的最大收藏价值,并且每种品牌至少购买一个产品。
显然,f[i][j][k]的最优解必定包括f[i-1][num[i-1]][k-w[i][j]]的最优解
或者f[i][j-1][k-w[i][j]]的最优解
或者f[i][j-1][k]的最优解
这也是能做出的3中决策——第j件物品为第i个品牌的第一件,不是第一件,或者不买第j件物品。(这些决策实行的条件为前面的转移决策合法)
决策合法当且仅当每种品牌至少购买一个产品
所以就有了状态转移方程式:
F[i][j][k]=max{F[i-1][num[i-1]][k-w[i][j]]+c[i][j] ,
F[i][j-1][k-w[i][j]]+c[i][j], F[i][j-1][k]}
此处因为每个品牌都要有所以不能从f[i-1][num[i-1]][k]转移
具体实现时我们可以把f数组先全部置为-1,表示不合法
初始时的合法状态为f[0][0][k]=0 num[0]=0
考虑到此处可以像01背包一样进行空间上的优化
考虑把第二维省掉(因为后面两个式子去掉第一维跟01背包一样)
注意点是第三维度要倒序循环
在考虑一下能否把第一维度省掉
不能,因为第二维度不能匹配,且转移的是i-1的num[i-1],可能被覆盖掉
不过可以考虑用滚动数组
时间复杂度:() 也就是 O(KNM)
空间复杂度:(KM) 或者 (M)
离奇之处:
if (dp[i][k-w[i][j]]!=-1)
dp[i][k]=max(dp[i][k],dp[i][k-w[i][j]]+c[i][j]);
if (dp[i-1][k-w[i][j]]!=-1)
dp[i][k]=max(dp[i][k],dp[i-1][k-w[i][j]]+c[i][j]);
这两句话不能反过来,即必须先做f[i][k-w[i][j]]
此处的缘由我也不知道,WA了很多次
有人知道后望告知我。。。
知道了。。原因是价格可能为0,所以dp[i][k]=max(dp[i][k],dp[i][k-w[i][j]]+c[i][j]);
放在后面的话假如价格为零就先当本品牌第一个产品买了一次,然后又买了一次,造成双次购买,结果偏大(出题人好坑。)(来自hk同志的想法)
CODE
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxk=11,maxm=10010,maxn=105;
int dp[maxk][maxm],w[maxk][maxn],c[maxk][maxn],num[maxk];
int n,m,k;
void init()
{
int a,b,cc;
cin>>n>>m>>k;
for (int i=1;i<=n;i++)
{
cin>>a>>b>>cc;
num[a]++;
w[a][num[a]]=b;
c[a][num[a]]=cc;
}
memset(dp,-1,sizeof(dp));
for (int i=0;i<=m;i++)
dp[0][i]=0;
}
void work()
{
for (int i=1;i<=k;i++)
for (int j=1;j<=num[i];j++)
for (int k=m;k>=w[i][j];k--)
{
if (dp[i][k-w[i][j]]!=-1)
dp[i][k]=max(dp[i][k],dp[i][k-w[i][j]]+c[i][j]);
if (dp[i-1][k-w[i][j]]!=-1)
dp[i][k]=max(dp[i][k],dp[i-1][k-w[i][j]]+c[i][j]);
}
}
void outit()
{
if (dp[k][m]==-1) puts("Impossible");
else cout<<dp[k][m];;
}
int main()
{
init();
work();
outit();
}