矩阵连乘问题

矩阵连乘问题

       给定n个矩阵｛A1,A2,…,An｝，其中Ai与Ai+1是可乘的，i=1,2 ,…,n-1。如何确定计算矩阵连乘积的计算次序，使得依此次序计算矩阵连乘积需要的数乘次数最少。例如：

　　A1={30x35} ; A2={35x15} ;A3={15x5} ;A4={5x10} ;A5={10x20} ;A6={20x25} ;

最后的结果为：((A1(A2A3))((A4A5)A6))  最小的乘次为15125。

样例输入

6

30 35 15 5 10 20 25

样例输出

15125

解法一（动态规划）：

矩阵连乘问题(动态规划)：
两个矩阵分别为：x*y,y*z则两个矩阵的乘法次数为（x*y*z）
所以多个矩阵的连乘分成，前k个的乘积和k~j的乘积相乘
使用dp[i][j]代表i~j矩阵之间相乘最少乘法次数 
r代表做乘法的矩阵数量，i开始矩阵，j为结束矩阵
初始就是，i矩阵与i+1~j矩阵做乘法
然后开始循环，k代表中间矩阵位置。
i~j -> i~k*k+1~j
则此时的乘法数目为：dp[i][k]+dp[k+1][j]+a[i-1]*a[i]*a[j] 
a[i-1]*a[i]*a[j]代表本次乘法的次数 

#include <stdio.h>
#include <string.h>
int main()
{
	int n;
	int r,i,j,k;
	int a[100];
	int dp[100][100];
	memset(dp,0,sizeof(dp));
	scanf("%d",&n); 
	for(i=0;i<=n;i++)
	{
		scanf("%d",&a[i]);
	}
	for(r=1;r<=n;r++)
	{
		for(i=1;i<=n-r;i++)
		{
			j=i+r;
			dp[i][j]=dp[i][i]+dp[i+1][j]+a[i-1]*a[i]*a[j];
			for(k=i;k<=j;k++)
			{
				int temp=dp[i][k]+dp[k+1][j]+a[i-1]*a[k]*a[j];
				if(temp<dp[i][j]){
					dp[i][j]=temp;
				} 
			}
		}
	}
	printf("%d\n",dp[1][n]);
	return 0;
 }