通过C中的动态编程解决背包问题的麻烦

时间:2015-05-09 15:05:36

标签: c dynamic-programming knapsack-problem

我需要以递归,备忘和动态编程来解决背包问题。目前我坚持使用动态编程方法。

我改编了我在互联网上其他地方找到的代码。目前输出不正确。

问题涉及利润和质量。每个项目都有相关的利润和质量,有一个MAX_N(数量)的项目和一个MAX_CAPACITY的质量。目的是尽可能在背包中获得“利润”。

以下是练习提供的示例:

  

示例:给出容量为5的背包,以及质量为[] = {2,4,3,2}的项目   和利润[] = {45,40,25,15},最佳组合是项目0(质量2和利润45)和项目2(质量3和利润25),总利润为70。没有其他质量为5或更低的组合有更大的利润。

以下是完整的代码:

#include <stdio.h>
#include <stdlib.h>

#define MAX_N 10
#define MAX_CAPACITY 165

int m[MAX_N][MAX_CAPACITY];

int max(int x, int y) {
    return x ^ ((x ^ y) & -(x < y));
}

int min(int x, int y) {
    return y ^ ((x ^ y) & -(x < y));
}

int knapsackRecursive(int capacity, int mass[], int profit[], int n) {

    if (n < 0)
        return 0;

    if (mass[n] > capacity)
        return knapsackRecursive(capacity, mass, profit, n-1);

    else
        return max(knapsackRecursive(capacity, mass, profit, n-1), knapsackRecursive(capacity - mass[n], mass, profit, n-1) + profit[n]);

}

int knapsackMemoized(int capacity, int mass[], int profit[], int n) {

}

int knapsackDynamic(int capacity, int mass[], int profit[], int n) {

    int i;
    int j;

    for (i = 0; i <= n; i++) {

        for (j = 0; j <= capacity; j++) {

            if (i == 0 || j == 0)
                m[i][j] = 0;

            else if (mass[i-1] <= j)
                m[i][j] = max(profit[i-1] + m[i-1][j-mass[i-1]], m[i-1][j]);

            else
                m[i][j] = m[i-1][j];
        }
    }

    return m[n][capacity];

}

void test() {

    // test values
    //int M1[MAX_N] = {2, 4, 3, 2};
    //int P1[MAX_N] = {45, 40, 25, 10};

    int M1[MAX_N] = {6, 3, 2, 4};
    int P1[MAX_N] = {50, 60, 40, 20};

    int M2[MAX_N] = {23, 31, 29, 44, 53, 38, 63, 85, 89, 82};
    int P2[MAX_N] = {92, 57, 49, 68, 60, 43, 67, 84, 87, 72};

    // a)
    printf("Recursion: %d\n",knapsackRecursive(MAX_CAPACITY, M1, P1, MAX_N));
    printf("Recursion: %d\n",knapsackRecursive(MAX_CAPACITY, M2, P2, MAX_N));
    printf("\n");

    // b)
    printf("Memoization: %d\n",knapsackMemoized(MAX_CAPACITY, M1, P1, MAX_N));
    printf("Memoization: %d\n",knapsackMemoized(MAX_CAPACITY, M2, P2, MAX_N));
    printf("\n");

    // c)
    printf("Dynamic Programming: %d\n",knapsackDynamic(MAX_CAPACITY, M1, P1, MAX_N));
    printf("Dynamic Programming: %d\n",knapsackDynamic(MAX_CAPACITY, M2, P2, MAX_N));

}

int main() {
    test();
}

这是我目前得到的输出。递归方法应该提供正确的结果,但动态编程当前不输出相同的结果。记忆尚未完成,因此也无法正确输出。

Recursion: 170
Recursion: 309

Memoization: 2686680
Memoization: 2686600

Dynamic Programming: 0
Dynamic Programming: 270

Process returned 25 (0x19)   execution time : 0.269 s
Press any key to continue.

1 个答案:

答案 0 :(得分:0)

事实证明,我用于编写动态编程部分的代码应该与int m[MAX_N+1][MAX_CAPACITY+1];而不是int m[MAX_N][MAX_CAPACITY];一起使用。

改变这一点让我得到了一个正常工作的代码,如果不是我想要的代码。