C++实现矩阵链乘法利用动态规划及运行实例结果

矩阵乘法满足结合律，所以哪里加括号不影响最终结果，括号加在不同的地方会对运算代价产生巨大影响，因此在进行矩阵链乘法时需要考虑那些矩阵先相乘以减少运算代价。下面使用动态规划方法实现矩阵链乘法

伪代码

c++实现代码

#include <iostream>
#define N 7

using namespace std;

void matrixChainOrder(int p[],int length,int m[][N],int s[][N])
{
    int i,l,j,k;
    int q;
    for(i=0;i<length;++i)
        m[i][i]=0;//记录矩阵所需标量乘法的最小值
    for(l=2;l<length;++l)//l为矩阵链的长度
    {
        for(i=1;i<length-l+1;++i)
        {
            j=i+l-1;
            m[i][j]=100000;
            for(k=i;k<=j-1;++k)
            {
               q=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];
               if(q<m[i][j])
               {
                   m[i][j]=q;
                   s[i][j]=k;//s[][]记录构造最优解所需信息
               }
            }
        }
    }
}

void printOptimalParents(int s[][N],int i,int j)//输出最优解函数
{
    if(i==j)
        cout<<"A"<<i;
    else
    {
        cout<<"(";
        printOptimalParents(s,i,s[i][j]);
        printOptimalParents(s,s[i][j]+1,j);
        cout<<")";
    }
}

int main()
{
    int p[N]={5,10,3,12,5,50,6};
    int m[N][N]={0},s[N][N]={0};
    matrixChainOrder(p,N,m,s);
    cout<<"m矩阵的值为"<<endl;
    for(int i=1;i<N;++i)
    {
        for(int j=1;j<N;++j)
            cout<<m[i][j]<<"\t";
        cout<<endl;
    }
    cout<<"s矩阵的值为"<<endl;
    for(int i=1;i<N;++i)
    {
        for(int j=1;j<N;++j)
            cout<<s[i][j]<<"\t";
        cout<<endl;
    }
    cout<<"矩阵链中括号放置顺序如下时运算代价最小"<<endl;
    printOptimalParents(s,1,N-1);

    return 0;
}

运行结果