矩阵链乘问题
对于给定的n个矩阵,M1, M2 ,…, Mn,其中矩阵Mi 和Mj 是可乘的,要求确定计算矩阵连乘积 ( M1M2 …Mn )的计算次序,使得按照该次数计算 矩阵连乘积时需要的乘法次数最少
1、描述最优解结构
目标:
求出矩阵链乘Mi Mi+1 ┅Mj-1 Mj(i<j) 所需的最少乘法次数
分析:
共有j-i+1个矩阵,按照做最后一次乘法的位置进行划分,矩阵链乘一共可分为j-i种情况
最优解可表示为:
(Mi ┅ Mk) (Mk+1 ┅ Mj)(i≤k<j)
2、递归地定义最优值(写出动态规划方程)
记第t个矩阵Mt的列数为rt
令r0为矩阵 M1的行数 ,则第t个矩阵的行数为rt-1
故Mi ┅ Mk连乘所得是ri-1rk维矩阵 ,Mk+1 ┅ Mj连乘所得是rkrj维矩阵
两这两个矩阵相乘做ri-1rkrj次乘法
记(Mi ┅ Mk)和(Mk+1 ┅ Mj)的最少乘法次数,记为m[i][k]和m[k+1][j]
故(Mi … Mk)(Mk+1 … Mj)的最少乘法次数为m[i][k]+m[k+1][j]+ri-1rkrj
将满足i≤k<j 的共j-i种情况逐一进行比较,
可知m[i][j]=min{m[i][k]+m[k+1][j]+ri-1rkr