在C ++程序崩溃中使用动态编程的矩阵链乘法?

时间:2014-09-30 12:22:53

标签: c++ algorithm dynamic-programming

我编写了以下C ++程序来实现使用动态编程实现MCM。但是下面的程序崩溃了。我的代码出了什么问题?

#include<iostream>
#include<cstdlib>
#define SZ 10
using namespace std;

int table[SZ][SZ];
int P[] = {2,3,3,5};

int MCM(int i, int j)
{
    if(i==j) return 0;

    else
    {
        int min = INT_MAX;

        for(int k=i;k<=j;k++)
        {
            if(table[i][k]==0)
                table[i][k] = MCM(i,k);
            if(table[k+1][j]==0)
                table[k+1][j] = MCM(k+1,j);
            int sum = table[i][k] + table[k+1][j] + P[i-1]*P[j]*P[k];
            if(sum<min)
                min = sum;  
        }
        return min;
    }
}

int main()
{
    int size = sizeof(P)/sizeof(P[0]);
    printf("Minimum number of mutiplications is %d",MCM(0,size-1));
    return 0;
}

3 个答案:

答案 0 :(得分:1)

你的代码将无限循环。除此之外,你犯了一些错误:

  1. 您从未在表格中指定最佳值(当您找到最小金额时,您不会存储它)。因此,每次检查表[i] [j] == 0时,它都是真的
  2. 循环中的
  3. k可以等于j而你使用k + 1,这是一个错误
  4. 无论如何,我认为你的代码的正确版本应该是这样的:

    #include<iostream>
    #include<cstdlib>
    #define SZ 10
    using namespace std;
    
    int table[SZ][SZ];
    int P[] = {1,2,3,4};
    
    int MCM(int i, int j)
    {
        if(i==j) return 0;
    
        else
        {
            int min = INT_MAX;
    
            for(int k=i;k<j;k++)
            {
                if(table[i][k]==0)
                    table[i][k] = MCM(i,k);
                if(table[k+1][j]==0)
                    table[k+1][j] = MCM(k+1,j);
                int sum = table[i][k] + table[k+1][j] + P[i]*P[j]*P[k];
                if(sum<min)
                    min = sum;
           }
           table[i][j] = min;
           return min;
       }
    

    }

    int main()
    {
        int size = sizeof(P)/sizeof(P[0]);
        printf("Minimum number of mutiplications is %d",MCM(0,size-1));
        return 0;
    }
    

答案 1 :(得分:0)

第一次调用MCM(0,size-1)时,参数i = 0,但是你减去1并使用结果-1来访问P [-1](第23行)。这可能会导致崩溃。

答案 2 :(得分:0)

哦我弄错了这只是一个小错误

第17行的循环应该以条件

结束
k<j 

而不是

k<=j

程序编译并成功运行

#include<iostream>
#include<cstdlib>
#define SZ 10
using namespace std;

int table[SZ][SZ];
int P[] = {1, 2, 3, 4, 3};

int MCM(int i, int j)
{
    if(i==j) return 0;

    else
    {
        int min = INT_MAX;

        for(int k=i;k<j;k++) // bug was here: for(int k=i;k<=j;k++)
        {
            if(table[i][k]==0)
                table[i][k] = MCM(i,k);
            if(table[k+1][j]==0)
                table[k+1][j] = MCM(k+1,j);
            int sum = table[i][k] + table[k+1][j] + P[i-1]*P[j]*P[k];
            if(sum<min)
                min = sum;  
        }
        return min;
    }
}

int main()
{
    int size = sizeof(P)/sizeof(P[0]);
    printf("Minimum number of mutiplications is %d",MCM(1,size-1));
    return 0;
}