矩阵连乘问题
给定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;
}