承接着上一篇的装配线调度问题,下一个要解决的问题是矩阵链乘法问题:给定n个要相乘的矩阵构成序列<A1,A2,...,An>,计算乘积A1A2...An。为了得到计算的结果,可以将两个矩阵相乘的算法作为一个子程序,根据括号的顺序做全部的矩阵乘法,矩阵的乘法满足结合律,无论怎么加括号都会得到相同的结果,但注意并不满足交换律。例如,4个矩阵的乘积可用5种不同的括号顺序表示:
但是要注意,尽管矩阵括号的顺序并不影响最后的结果,但对运算的代价有很大的影响,先观察矩阵的标准算法:
如果A是p*q矩阵,B是q*r矩阵,那么结果的矩阵C就是p*r矩阵,计算C的时间由第七行的标量乘法次数决定,这个次数等于p*q*r。为了说明加上不同的括号计算的时间是不同的,假设是3个矩阵进行相乘<A1,A2,A3>,维度分别是10*100,100*5,5*50,假设(A1*A2)*A3进行计算的话,那么次数将是10*100*5+10*5*50=7500,但假如是按照A1*(A2*A3)顺序进行计算的话,那么总的次数将是100*5*50+10*100*50=50000,显然使用前一种括号的顺序可以使得速度快了近10倍。
所以,矩阵链乘法问题可以表述如下:给定n个矩阵构成的一个链<A1,A2,...,An>,其中i=1,2,...,n,矩阵Ai的维数为P(i-1)*Pi,对乘积A1A2...An以一种最小化标量乘法次数进行加括号。