矩阵连乘问题

问题叙述： 

        给定n个矩阵（A1,A2,A3.....An}，其中Ai与Ai+1是可乘的，i=1,2,...n-1。考察n个矩阵的连乘积A1A2A3,....An。由于矩阵乘法满足结合律，故计算矩阵的连乘积可以有许多不同的计算次序，这种计算次序可以用加括号的方式来确定。加括号的方式决定了整个计算量（指的是乘法调用的次数）。所以自然会提出矩阵连乘积的最优计算次序问题。

#include<iostream>
using namespace std;
#define N 7 

//m[][]存放数乘次数 
//s[][]存放断开位置比如S[1][5]=3可表示：(A1A2A3)(A4A5)，当然内部断链还会继续划分A1A2A3 
void MatrixChain(int *p,int n,int m[][N],int s[][N]){
	for(int i=1;i<=n;i++){
		m[i][i]=0;      //本身不需要数乘设置为0 
	}
	for(int r=2;r<=n;r++){    //r指矩阵链长度，从2开始 
		for(int i=1;i<=n-r+1;i++){    //i指断开的位置 
			int j=i+r-1;   //矩阵链的末尾矩阵   
			m[i][j]=m[i][i]+m[i+1][j]+p[i-1]*p[i]*p[j];
			s[i][j]=i;
			//改变断点，试探出最小的情况(相当于内部断开查找(A1A2A3)可以断开的位置)
			for(int k=i+1;k<j;k++){ //这里面将断链点从i+1开始，可以断链的点直到j-1为止
				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;  //记录内部的断点位置
				}
			}
		}
	}
} 

void Traceback(int i,int j,int s[][N]){
	if(i==j)       //回归条件
    {
        cout<<"A"<<i;
    }
    else       //按照最佳断点一分为二，接着继续递归
    {
        cout<<"(";
        //判断 1,6,s→1,3,s →1,1,s
        Traceback(i,s[i][j],s);
        // 判断2,3,s →2,2,s
        Traceback(s[i][j]+1,j,s);
        cout<<")";
    }
}

int main(){
	// A1=30*35  A2=35*15  A3=15*5...... 
	int p[N]={30,35,15,5,10,20,25};

	int m[N][N],s[N][N];
	MatrixChain(p,N-1,m,s);
	Traceback(1,6,s); 
}