矩阵连乘问题

一、问题描述：

给定n个矩阵｛A1,A2,…,An｝，其中Ai与Ai+1是可乘的，i=1，2…，n-1。如何确定计算矩阵连乘积的计算次序，使得依此次序计算矩阵连乘积需要的数乘次数最少。例如，给定三个连乘矩阵{A1，A2，A3}的维数分别是10*100，100*5和5*50，采用（A1A2）A3，乘法次数为10*100*5+10*5*50=7500次，而采用A1（A2A3），乘法次数为100*5*50+10*100*50=75000次乘法，显然，最好的次序是（A1A2)A3，乘法次数为7500次。

要求：给定由n个矩阵构成的序列｛A1，A2，...，An}，对乘积A1A2...An，找到最小化乘法次数的加括号方法。

二、问题分析

将矩阵连乘积Ai，Ai+1，Ai+2……Aj，简记为A[i:j]，这里i<=j。

如果要得到A[i:j]的最优计算次序，计算的最后-步一定是经过若 千次计算后形成的两个矩阵相乘，如下式所示: .(AiAi+1+... Ak)(Ak+1Ak+2……Aj)。

即这个计算次序的最后一步在矩阵Ak和Ak+1之间分割为两部分进行矩阵相乘，i

在这j-i种分割中，一定有一种是导致最优结果的。

如果计算A[i:j]的最优次序时，从A[k]分割开是最优那么计算矩阵子链A[i:k] 和A[k+1:j]的次序也一定是最优的。因此我们可以递归的求出所有j-i种分割所导致的最优计算量，通过比较选出导致最少乘法次数的切割方式，将问题解决。

最终乘法计算次数=得到A[i:k]时产生的乘法计算次数+得到A[k+1:j]时产生的乘法计算次数+A[i:k]和A[k+1:j]相乘的乘法计算次数。

代码如下：

#include

#include

#include

#include

using namespace std;

/*

*矩阵连乘（备忘录方法：自顶向下递归）

*/

vector > m;//m[i][j]表示矩阵Ai连乘到Aj的最少运算次数

vector > s;//s[i][j]记录矩阵Ai和矩阵Aj之间的分割点

//计算该连乘式子的最佳结合方式

int MatrixChain(vector& p,int beg, int end)

{

   if(m[beg][end]>0) return m[beg][end];

   if(beg==end) return 0;

   int u = MatrixChain(p,beg,beg) +MatrixChain(p,beg+1,end)+p[beg-1]*p[beg]*p[end];

   s[beg][end] = beg;

   for (int K = beg+1; K

   {

      int t = MatrixChain(p,beg,K) + MatrixChain(p,K+1,end) + p[beg-1]*p[K]*p[end];

      if (t

      {

          u = t;s[beg][end] = K;

      }

   }

   m[beg][end] = u;

   return u;

}

//输出该连乘式子的最佳结合方式

void PrintMatrixChain(int n,int m)

{

   if(n==m)

   {

      cout<<"A"<

      return;

   }

   int k = s[n][m];

   if(n==k)

      PrintMatrixChain(n,k);

   else

   {

      cout<<"(";

      PrintMatrixChain(n,k);

      cout<<")";

   }

   if(k+1==m)

      PrintMatrixChain(k+1,m);

   else

   {

      cout<<"(";

      PrintMatrixChain(k+1,m);

      cout<<")";

   }

}

int main()

{

   vector vec;

copy(istream_iterator(cin),istream_iterator(),back_inserter(vec));

   int n = vec.size()-1;//一共有n个矩阵相乘

   m = vector >(n+1,vector(n+1,0));//0行0列空余

   s = vector >(n+1,vector(n+1,0));//0行0列空余

   //初始化m数组

   for(int i = 0;i<=n;i++) m[i][i] = 0;

   int u = MatrixChain(vec,1,n);

   cout<<"最优解为计算"<

   PrintMatrixChain(1,vec.size()-1);

Return 0；

}