矩阵链乘法

1.问题

设A1,A2,…An为n个矩阵的序列，其中Ai为Pi-1×Pi阶矩阵，这个矩阵链的输入用向量P=<P0,P1,…,Pn>给出

给定向量P，确定一种乘法次序，使得基本运算的总次数达到最小

例如P=<10,20,30,50>,则A1:10×20，A2:20×30，A3:30×50
（1） (A1A2)A3:10×20×30+10×30×50=21000
（2）A1(A2A3):10×20×50+20×30×50=40000

2.解析

蛮力法
枚举所有可能的乘法次序，针对每种次序计算基本运算的次数，从中找出具有最小 运算次数的乘法次序，每一种乘法次序对应了一种在 n 个项中加 n-1 对括号。

动态规划法


m[i,j]=min{m[i,k]+m[k+1,j]+Pi=1PkPj}满足优化原则。也就是说，m[i,j]最小值时，m[i,k],m[k+1,j]也是最小的。

3.设计

迭代算法

MatrixChain(P,n)
输入:矩阵链Ai….j的输入向量P=<Pi-1,Pi,…,Pj>  1≤i≤j≤n
输出:计算Ai….j所需最小乘法运算次数m[i..j]和最后一次运算的位置s[i..j] (s[i..j]= k,(Ai…Ak)(Ak+1…Aj))
1.令所有 m[i,j]初值为0, s[i,j]的初值为i， 1≤i≤ j≤n
2. For r=2 to n do   //r为当前问题规模(长度)
3.	For i=1 to n-r+1 do   //i的起点不断变化，各种r长
4.		j=i+r-1  //不同终点
5.		m[i,j]=m[i+1,j]+Pi-1PkPj //划分为4.4..)
6.		s[i,j]=i
7.		For k=i+1 to j-1 do	 //不同的划分位置
8.			t= m[i,k]+m[k+l,j]+ Pi-1PkPj
9.			if t<m[i,j]
10.			Then m[i,j]=t
11.			s[i,j]=k

4.分析

主要时间复杂度在主函数的三重循环中，所以时间复杂度为：T=O(n^3)

5.源码

矩阵链乘法代码