动态规划 | 矩阵链乘法

问题

设A1, A2, …, An为n个矩阵的序列，其中Ai为Pi-1 * Pi阶矩阵，这个矩阵链的输入用向量P=<P0, P1, …, Pn>给出。给定向量P，确定一种乘法次序，使得基本运算的总次数达到最小。

解析

示例：n=6, P=<10, 5, 15, 20, 4, 6, 8>
A1: 105
A2: 515
A3: 1520
A4: 204
A5: 46
A6: 68

(1) r=1
m[1][1]=0
m[2][2]=0
…

(2) r=2 i=1,…,5 j=2,…,6
m[1][2]=10515=750
m[2][3]=51520=1500
m[3][4]=15204=1200
m[4][5]=2046=480
m[5][6]=468=192

(3) r=3 i=1,…,4 j=3,…,6
m[1][3]=min{m[1][1]+m[2][3]+A1(A2A3),m[1][2]+m[3][3]+(A1A2)A3}=min(2500, 3750)=2500
m[2][4]=min{m[2][2]+m[3][4]+A2(A3A4),m[2][3]+m[4][4]+(A2A3)A4}=min(1500, 1900)=1500
m[3][5]=min{m[3][3]+m[4][5]+A3(A4A5),m[3][4]+m[5][5]+(A3A4)A5}=min(2280,1560)=1560
m[4][6]=min{m[4][4]+m[5][6]+A4(A5A6),m[4][5]+m[6][6]+(A4A5)A6}=min(832,1440)=832

(4) r=4 i=1,…,3 j=4,…,6
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}=min(1700,3900,3300)=1700
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}=min(2010,3960,1620)=1620
m[3][6]=min{m[3][3]+m[4][6]+A3(A4A5A6),m[3][4]+m[5][6]+(A3A4)(A5A6),m[3][5]+m[6][6]+(A3A4A5)A6}=min(3232,2784,2280)=2280

(5) r=5 i=1,…,2 j=5,…,6
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}=min(1920,3420,4180,1940)=1920
m[2][6]=min{m[2][2]+m[3][6]+A2(A3A4A5A6),m[2][3]+m[4][6]+(A2A3)(A4A5A6),m[2][4]+m[5][6]+(A2A3A4)(A5A6),m[2][5]+m[6][6]+(A2A3A4A5)A6}=min(2520,3132,2012,1860)=1860

(6) r=6 i=1 j=6
m[1][6]=min{m[1][1]+m[2][6]+ A1(A2A3A4A5A6), m[1][2]+m[3][6]+(A1A2)(A3A4A5A6), m[1][3]+m[4][6]+(A1A2A3)(A4A5A6), m[1][4]+m[5][6]+(A1A2A3A4)(A5A6),m[1][5]+m[6][6]+(A1A2A3A4A5)A6}=min(2260,4230,4932,2212,2400)=2212

设计

算法核心伪代码：

for(int r=2;r<=n;++r){
	for(int i=1;i<=n-r+1;++i){
		int j=i+r-1;
		m[i][j]=m[i+1][j]+P[i-1]*P[i]*P[j];
		s[i][j]=i;
		for(int k=i+1;k<=j-1;++k){
			int t=m[i][k]+m[k+1][j]+ P[i-1]*P[k]*P[j];
			if(t<m[i][j]){
				m[i][j]=t;
				s[i][j]=k;
			}
		}
	}
}

分析

时间复杂度：T(n)=O(n^3)

源码

GitHub