HDU 2639 Bone Collector II （背包 第K优解）

题意：给出N,V,K、 分别代表 物品的数量，背包的体积和满足条件的第K个大的值、  接下来两行分别有N个数字，分别代表着每个的价值，和体积、 最后求出第K个大的值为多少、

思路：根据背包九讲中的思路来进行思考、 如果对背包不是特别理解的可以看一下背包九讲、 链接：http://blog.csdn.net/liuqiyao_01/article/details/8521776 

对于求次优解、第K优解类的问题，如果相应的最优解问题能写出状态转移方程、用动态规划解决，那么求次优解往往可以相同的复杂度解决，第K优解则比求最优解的复杂度上多一个系数K。

其基本思想是将每个状态都表示成有序队列，将状态转移方程中的max/min转化成有序队列的合并。这里仍然以01背包为例讲解一下。

首先看01背包求最优解的状态转移方程：f[i][v]=max{f[i-1][v],f[i-1][v-c[i]]+w[i]}。如果要求第K优解，那么状态f[i][v]就应该是一个大小为K的数组f[i][v][1..K]。其中f[i][v][k]表示前i个物品、背包大小为v时，第k优解的值。“f[i][v]是一个大小为K的数组”这一句，熟悉C语言的同学可能比较好理解，或者也可以简单地理解为在原来的方程中加了一维。显然f[i][v][1..K]这K个数是由大到小排列的，所以我们把它认为是一个有序队列。

然后原方程就可以解释为：f[i][v]这个有序队列是由f[i-1][v]和f[i-1][v-c[i]]+w[i]这两个有序队列合并得到的。有序队列f[i-1][v]即f[i-1][v][1..K]，f[i-1][v-c[i]]+w[i]则理解为在f[i-1][v-c[i]][1..K]的每个数上加上w[i]后得到的有序队列。合并这两个有序队列并将结果（的前K项）储存到f[i][v][1..K]中的复杂度是O(K)。最后的答案是f[N][V][K]。总的复杂度是O(NVK)。

为什么这个方法正确呢？实际上，一个正确的状态转移方程的求解过程遍历了所有可用的策略，也就覆盖了问题的所有方案。只不过由于是求最优解，所以其它在任何一个策略上达不到最优的方案都被忽略了。如果把每个状态表示成一个大小为K的数组，并在这个数组中有序的保存该状态可取到的前K个最优值。那么，对于任两个状态的max运算等价于两个由大到小的有序队列的合并。

另外还要注意题目对于“第K优解”的定义，将策略不同但权值相同的两个方案是看作同一个解还是不同的解。如果是前者，则维护有序队列时要保证队列里的数没有重复的。

本人思路：  普通的01背包中，我们基本求解的都是当前背包容量为V时，的最优解，也就是满足背包容量的最大价值、  而这道题求得是第K大的值， 那么我们就不能用max来求解这个问题、 因为，max里面求得是，每装入一个物品，会把价值最大的那个物品放入背包，否则不放入、 那么回漏解、会漏掉不满足最大价值的部分解、 所以，这道题不应该使用max 、 那么首先，背包问题有两种情况，就是装入，或者不装入、也就是dp[i][j] 和dp[i][j-v[i]]+w[i];  如果想要求出他的第k个最大值，首先我们要求出他的每一个状态、也就是装入，或者不装入我们都要求出、 然后在两个序列里面合并、 合并成一个有序的数列放入DP中、 那么最后求出的答案便是dp[V][K];  这道题要注意的是去重问题、 如果dp[j][c]==dp[j][c-1]  这就代表重复、那么此时我不应该c++、   还有就是dp[i][j]代表的是在体积为i的大小时，第K个大的值是多少、

AC代码 ：

#include<cstdio>
#include<cstring>
const int maxn=1100;
const int maxm=35;
int T;
int a,b,c;
int N,V,K;
int dp[maxn][maxm];
int w[maxn],v[maxn];
int A[maxm],B[maxm];

int Init()
{
    memset(dp,0,sizeof(dp));
    scanf("%d %d %d",&N,&V,&K);
    for(int i=0;i<N;i++)
        scanf("%d",&w[i]);
    for(int i=0;i<N;i++)
        scanf("%d",&v[i]);
}

void Dp()
{
    for(int i=N-1;i>=0;i--){
        for(int j=V;j>=v[i];j--){
            for(int k=1;k<=K;k++){
                A[k]=dp[j-v[i]][k]+w[i];
                B[k]=dp[j][k];
            }
            A[K+1]=-1;B[K+1]=-1;
            a=b=c=1;
            while(c<=K&&(A[a]!=-1||B[b]!=-1)){
                if(A[a]>B[b]){
                    dp[j][c]=A[a],a++;
                }
                else{
                    dp[j][c]=B[b],b++;
                }
                if(dp[j][c]!=dp[j][c-1]) c++;
            }
        }
    }
    printf("%d\n",dp[V][K]);
}

int main()
{
    scanf("%d",&T);
    while(T--){
        Init();
        Dp();
    }
    return 0;
}