1.问题
矩阵链乘法
- 设 A1, A2, … , An 为 n 个矩阵的序列,其中 Ai 为 Pi-1 * Pi阶矩阵,这个矩阵链的输入用向量 P = < P0, P1, … , Pn> 给出。
- 给定向量 P, 确定一种乘法次序,使得基本运算的总次数达到最小。
2.解析
- 动态规划法
Ai…j:表示矩阵乘法的子问题;
m[i,j]:表示得到乘积Ai…j所用最少基本运算次数;
s[i,j]:表示得到乘积Ai…j所用最少基本运算次数最外层括号的位置。
命题 m[i,j] = min{ m[i,k] + m[k,j] + Pi-1PkPj}满足优化原则,即m[i…j]最小值时,m[i,k]和m[k+1,j]也是最小的。
示例:
P=<5,6,2,9,7,6>,n=5
- A1: 5×6
- A2: 6×2
- A3: 2×9
- A4: 9×7
- A5: 7×6
(1)r=1
- m[1,1] = 0;
- m[2,2] = 0;
- m[3,3] = 0;
- m[4,4] = 0;
- m[5,5] = 0;
(2)r=2,i=1,2,3,4; j=2,3,4,5
- m[1,2] = 5×6×2 = 60;
- m[2,3] = 6×2×9 = 108;
- m[3,4] = 2×9×7 = 126;
- m[4,5] = 9×7×6 = 378;
(3)r=3,i=1,2,3;j=3,4,5
- m[1,3]=min{m[1,2]+m[3,3]+(A1A2)A3,m[1,1]+m[2,3]+A1(A2A3)} = 150; s[1,3]=2;
- m[2,4]=min{m[2,3]+m[4,4]+(A2A3)A4,m[2,2]+m[3,4]+A2(A3A4)} = 210; s[2,4]=2;
- m[3,5]=min{m[3,4]+m[5,5]+(A3A4)A5,m[3,3]+m[4,5]+A3(A4A5)} = 210; s[3,5]=4;
(4)r=4,i=1,2;j=4,5
- m[1,4]=min{m[1,1]+m[2,4]+A1(A2A3A4),m[1,2]+m[3,4]+(A1A2)(A3A4),
- m[1,3]+m[4,4]+(A1A2A3)A4} = 326; s[1,4]=2;
- m[2,5]=min{m[2,2]+m[3,5]+A2(A3A4A5),m[2,3]+m[4,5]+(A2A3)(A4A5),
m[2,4]+m[5,5]+(A2A3A4)A5} = 282; s[2,5]=2;
(5)r=5,i=1;j=5
- m[1,5]=min{m[1,1]+m[2,5]+A1(A2A3A4A5),m[1,2]+m[3,5]+(A1A2)(A3A4A5),
m[1,3]+m[4,5]+(A1A2A3)(A4A5),m[1,4]+m[5,5]+(A1A2A3A4)A5} = 330; s[1,5]=2;
s[1,5] = 2: (A1A2)(A3A4A5)
s[3,5] = 4: (A3A4)A5
(A1A2)((A3)A4A5)
3.设计
迭代算法
核心伪代码
MatrixChain(P,n)
输入:矩阵链Ai…j的输入向量P=<Pi-1,Pi,…Pj> 1≤i≤j≤n
输出:计算Ai…j所需最小乘法运算次数m[i…j]和最后一次运算的位置s[i…j]
令所有m[i,j]初值为0,s[i,j]的初值为i 1≤i≤j≤n
For r=2 to n do //r为当前问题规模(长度)
For i=1 to n-r+1 do //i的起点不断变化,各种r长
j=i+r-1 //不同终点
m[i,j]=m[i+1,j]+Pi-1PkPj //划分为Ai(Ai+1…Aj)
s[i,j]=i
For k=i+1 to j-1 do
t =m[i,k]+m[k+1,j]+Pi-1PkPj //不同的划分位置
if t<m[i,j]
Then m[i,j]=t
s[i,j]=k
4.分析
三层循环,算法时间复杂度为O(n³)