這個問題是動態規劃的基礎的問題，也是演算法導論中討論過的問題。在這裡先簡單描述一下。假定有一組矩陣需要做乘法操作。但是我們知道首先矩陣乘法滿足了結合律。所以可以按照不同的順序做乘法。而且不同順序做乘法最後的乘法次數是不同的。比如〈A₁, A₂, A₃〉分別是10 × 100, 100 × 5, 和 5 × 50。如果按照((A₁A₂) A₃)的順序來計算,就是7500次，但是如果(A₁ (A₂A₃))這樣的順序，那結果就是75000次！！！這裡有10倍的差距。所以這個題目的意思就是找到一個方案使得需要做乘法的次數最少。

現在我們嘗試用大白話來繼續一下。問題的解決方案側重在兩個方面。一個是最少乘法的次數，另一個就是最小乘法次數的流程。在程式碼中我們分別用m,s來表示。先說一下m的含義。m是個二維陣列，m[i,j]表示從矩陣i到矩陣j的最小乘法次數。也就是A[i],A[i+1]....A[j]的最小乘法次數。可以發現，這個次數應該等於A[i],A[i+1]...A[k] * A[k+1]....A[j-1]A[j]的次數，這裡的k是輸入i到j這個開區間的某個數。找到這個數字，我們也就找到了最小乘法次數。

那怎麼找到這個數字呢？簡單的說就是把所有的可能性都列舉出一遍k=i, i+1, i+2...j都嘗試一遍。比較最小值就可以了。有了這個思路我們可以把遞迴公式寫成:

這裡的p就是你的矩陣維數的陣列。

根據這個公式，我們可以寫出下面的程式碼：

#include "stdafx.h" #include <vector> #include <limits> #include <iostream> typedef std::vector< std::vector<int> > intMatrix; /** * Description: Calulate the matrix chain multiplication order *@param p, order of the matrix demension; *@param m, output value, Save the number of the multiplication for specific i,j. * for instance, the number for n matrix will be stored in m[0][n-1]. *@param s, output value, Save the order of multiplication * for instance, s[i,j]=k, means matrix need to multiple in order m[i...k]*m[k+1....j]; *@return number of the multiplication times. */ int matrixChainOrder (const std::vector<int>& p, std::vector < std::vector<int> >& m, std::vector < std::vector<int> >&s) { int n = (int)p.size()-1; m.resize(n); for (int i=0;i<n;i++) m[i].resize(n); s.resize(n); for (int i=0;i<n;i++) s[i].resize(n); for (int l=2;l<=n;l++) { // length of the chain is l for (int i=0;i<n-l+1;i++) { int j=i+l-1; m[i][j] = std::numeric_limits<int>::max(); for (int k=i;k<j;k++) { int q=m[i][k]+m[k+1][j]+p[i]*p[k+1]*p[j+1]; if(q<m[i][j]) { m[i][j]=q; s[i][j]=k; } } } } return m[0][n-1]; } void printSequence(const intMatrix& s, int i, int j) { if (i<j) { std::cout<<"("; printSequence(s, i, s[i][j]); printSequence(s, s[i][j]+1,j); std::cout<<")"; } else std::cout<< "M["<<i<<"]"; } int _tmain(int argc, _TCHAR* argv[]) { int k[]={30,35,15,5,10,20,25}; std::vector<int> p(k,k+sizeof(k)/sizeof(k[0])); intMatrix m; intMatrix s; // number of multiplication std::cout<<"Minimum multiplication times is "<<matrixChainOrder(p, m, s)<<std::endl; // backtrack the sequence printSequence(s, 0,(int)p.size()-2); getchar(); return 0; }

需要注意的問題是，這個演算法的複雜度在O(n∧3)。而且這個演算法和演算法導論等地方的偽碼比有些細微的差異體現在陣列從0開始還是1開始上。大體上還是一致的。這裡的空間複雜度也到了O(n∧2)。大致的計算順序就是通過這個函式中的l來控制。l=2的時候，依次計算矩陣A[0]*A[1], A[1]*A[2], A[2]*A[3]...的次數。l=3的時候就開始計算A[0]*A[1]*A[2], A[1]*A[2]*A[3]..的次數。用更形象的話來說就是倒三角的順序。下面的這個圖就是對上面流程的描述：

最後說一句，這個矩陣鏈問題是相對獨立的一個問題，和其他的知名的問題都沒有太大的關聯。以後的總結可能有很多不同的變形。