1. 问题
设A1,A2,…,An为n个矩阵的序列,其中Ai为Pi-1×Pi阶矩阵,这个矩阵链的输入用向量P=<P0,P1,…,Pn>给出。
给定向量P,确定一种乘法次序,使得基本运算的总次数达到最小。
例如,P=<10,100,5,50>,则A1:10×100,A2:100×5,A3:5×50
1)(A1A2)A3:10×100×5+10×5×50=7500
2)A1(A2A3):10×100×50+100×5×50=75000
2. 解析
动态规划法
Ai…j:表示矩阵链相乘的子问题AiAi+1…Aj;
m[i…j]:表示得到乘积Ai…j所用的最少基本运算次数;
假定,最后一次相乘发生在矩阵链Ai…k和Ak+1…j之间,即
AiAi+1…Aj=(AiAi+1…Ak)×(Ak+1Ak+2…Aj) k=I,i+1,…,j-1
e.g. Ai(Ai+1…Aj);(AiAi+1)(Ai+2…Aj);…;(AiAi+1Ai+2…)(Aj-1Aj);(AiAi+1Ai+2…Aj-1)Aj
0 i=j
m[i…j]={
min{m[I,k]+m[k+1,j]+Pi=1PkPj}(i<=k<j) i<j
其中Ai=Pi-1×Pi,Ak=Pk-1×Pk,Aj=Pj-1×Pj
AiAi+1…Aj=(AiAi+1…Ak)(Ak+1Ak+2…Aj)
→(Pi-1Pk)(PkPj)
→Pi-1PkPj
命题m[i…j]= min{m[I,k]+m[k+1,j]+Pi=1PkPj}(i<=k<j) 满足优化原则,即m[i…j]最小值时,m[i,k]和m[k+1,j]也是最小的。
3.设计
MatrixChain(P,n)
输入:矩阵链A1..</