矩阵连乘

---

问题：

---

给定n个矩阵：A1,A2,…,An，其中Ai与Ai+1是可乘的，i=1，2…，n-1。确定计算矩阵连乘积的计算次序，使得依此次序计算矩阵连乘积需要的数乘次数最少。输入数据为矩阵个数和每个矩阵规模，输出结果为计算矩阵连乘积的计算次序和最少数乘次数。

---

解析：

***由于矩阵乘法满足结合律，故计算矩阵的连乘积可以有许多不同的计算次序。这种计算次序可以用加括号的方式来确定。若一个矩阵连乘积的计算次序完全确定，也就是说该连乘积已完全加括号，则可以依此次序反复调用2个矩阵相乘的标准算法计算出矩阵连乘积。

---

具体代码

---

public class Matrix {
	public static void main(String[] args) {
		int []p = {30,35,15,5,10,20,25};
		int [][]m = new int [p.length+1][p.length+1];
		int [][]s = new int [p.length+1][p.length+1];
		MatrixChain(p,m,s);
	    System.out.println("乘法的最优次序：");
	    traceback(s,1,p.length-1);
	}
	public static void MatrixChain(int []p,int [][]m,int [][]s){
		
		for(int i=1;i<=p.length;i++){
			m[i][i]=0;
		}		
		for(int r = 2; r<= p.length;r++){
			for(int i=1;i<=p.length-r;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;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;
				   }
			   }
			}	
		}
	}
	private static void traceback(int[][] s, int i, int j) {
				if(i==j) 
					System.out.print("A"+i);    
				else if(i+1==j) 
					System.out.print(" (A"+i+" * "+" A"+j+") ");   
				else{     
					System.out.print(" (");   
					traceback(s,i,s[i][j]);   
					traceback(s,s[i][j]+1,j);      
					System.out.print(") ");    
				}		
	}
}