动态规划 -- 矩阵链乘法

参考资料

1 算法导论

2 xiaowei_cqu

3 ujs_abc

矩阵链乘法简介

对于给定的由n个要相乘的矩阵构成的序列 A1, A2, ..., An，要计算乘积 A1 A2 ... An 。根据矩阵的性质，可以在矩阵之间添加括号，最后的结果都是相同的。但是矩阵链中添加括号的顺序对求积运算的代价有很大的影响（可参考算法导论动态规划章节）。矩阵链乘法所要解决的问题是通过怎样的矩阵组合求出最小的乘法次数。

动态规划算法

通过寻找最优加全部括号的结构是动态规划的第一步，即将问题分割为两个子问题（最优加全部括号Ai, Ai+1,...Ak 和 Ak,Ak+1 ....Aj），寻找子问题的实例的最优解，然后合并得到整个问题的最优解。

该问题的递归定义为：   

    

代码实现

void matrixchain(int *p,int m[][N],int s[][N])
{
    int i ,j , len , q;
    for(i= 0; i < N ; j++)
        m[i][i] = 0;   //m 用来记录矩阵连乘的最小代价的
    for(len = 2; len <= N; len++)
    {
        for(i = 0; i < N-len + 1; i++)
        {
            j = i + len -1;
            m[i][j] = inf;
            //计算Ai -Aj最小代价
            for(k = i; k < j; k++)
            {
                q = m[i][k] + m[k+1][j] + p[i]*p[k+1]*p[j+1];
                if(q <m[i][j])
                {
                    m[i][j] = q;
                    s[i][j] = k;
                }
            }
        }
    }
}

void PrintOptimalParens(int s[][LENGTH], int i, int j)
{
    if(i == j)
        printf("A%d", i);
    else
    {
        printf("(");
        PrintOptimalParens(s, i, s[i][j]);
        PrintOptimalParens(s, s[i][j]+1, j);
        printf(")");
    }
}