矩阵连乘问题（动态规划算法+备忘录方法）

动态规划算法

#include <iostream>

using namespace std;

const int N=6;
int p[]={30,35,15,5,10,20,25};
int m[N+1][N+1];
int s[N+1][N+1];

void matrixChain(int *p,int m[][N+1],int s[][N+1])
{
    //初始化
    for(int i=1;i<=N;i++) m[i][i]=0;
    //矩阵连乘问题规模，r个连乘；
    for(int r=2;r<=N;r++)
        for(int i=1;i<=N-r+1;i++)       //遍历矩阵连乘开始矩阵
        {
            int j=i+r-1;                //确定矩阵连乘结束矩阵。
            m[i][j]=m[i][i]+m[i+1][j]+p[i-1]*p[i]*p[j];       //k此时为i；
            s[i][j]=i;
            for(int k=i+1;k<j;k++)
            {
                int t=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];
                if(t<m[i][j])
                {
                    m[i][j]=t;
                    s[i][j]=k;
                }
            }
        }
}

void traceback(int i,int j,int s[][N+1])
{
    if(i==j||i==N)
    {
        cout<<"A"<<i;
        return ;
    }
    cout<<"(";
    traceback(i,s[i][j],s);
    cout<<")(";
    traceback(s[i][j]+1,j,s);
    cout<<")";
}

int main()
{

    matrixChain(p,m,s);
    cout<<m[1][N]<<endl;
    traceback(1,N,s);
    cout<<endl;
    return 0;
}

备忘录方法

#include<iostream>

using namespace std;

int m[7][7],s[7][7];
int p[]={30,35,15,5,10,20,25};
const int N=6;

int Lookupchain(int i,int j)
{
    if(m[i][j])
        return m[i][j];
    if(i==j)
        return 0;
    int u=Lookupchain(i+1,j)+p[i-1]*p[i]*p[j];
    s[i][j]=i;
    for(int k=i+1;k<j;k++)
    {
        int t=Lookupchain(i,k)+Lookupchain(k+1,j)+p[i-1]*p[k]*p[j];
        if(t<u)
        {
            u=t;
            s[i][j]=k;
        }
    }
    m[i][j]=u;
    return u;
}

void traceback(int i,int j)
{
    if(i==j)
    {
        cout<<"A"<<i;
        return ;
    }

    cout<<"(";
    traceback(i,s[i][j]);
    cout<<")(";
    traceback(s[i][j]+1,j);
    cout<<")";
}

int main()
{

    Lookupchain(1,N);
    cout<<m[1][N]<<endl;
    return 0;
}