题目:给定n个矩阵{A1,A2,...,An},其中Ai与Ai+1是可乘的(i=1,2,...,n-1),考察这n个矩阵的连乘积的最优计算次序
举例说明矩阵结合方式对数乘次数的影响:
矩阵连乘积A1A2A3,3个矩阵的维数分别为10x100,100x5和5x50,连乘时加括号的方式有:
((A1*A2)*A3) 数乘次数:10*100*5+10*5*50=7500
(A1*(A2*A3)) 数乘次数:100*5*50+10*100*50=75000
以矩阵链ABCD为例
按照矩阵链长度递增计算最优值
矩阵链长度为1时,分别计算出矩阵链A、B、C、D的最优值
矩阵链长度为2时,分别计算出矩阵链AB、BC、CD的最优值
矩阵链长度为3时,分别计算出矩阵链ABC、BCD的最优值
矩阵链长度为4时,计算出矩阵链ABCD的最优值
比如有六个矩阵
思路步骤为:
- 分析最优解的结构
- 建立递归关系
- 计算最优值
以下是本人手写推算过程
详细代码(malloc动态分配内存)
#include <stdio.h>
#include <limits.h>
#include <stdlib.h> // 需要包含这个头文件来使用malloc
// 定义最大矩阵数(为了简化示例,我们假设有一个固定的最大矩阵数)
#define MAX_MATRICES 100
int MatrixChain(int **m,int **s,int *p,int n)
{
// m是最优值二位数组,s是最优值断开点的二维数组
// p是储存矩阵行列数的数组,用于计算数乘次数
// n是矩阵个数
for(int i=1;i<=n;i++)
{
m[i][i]=0;
}
for(int r=2;r<=n;r++)//对角线循环,r=2,3,4,5,6
{
for(int i=1;i<=n-r+1;i++)//行循环,i=1~5到i=1~1
{
int j=i+r-1;//控制列,每一个i对应一个j,m[i][j]
m[i][j]=m[i+1][j]+p[i-1]*p[i]*p[j];//默认以i为断电时为最优值
s[i][j]=i;//默认断电索引为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];
// 尝试其他断点从i+1到j-1,逐个遍历
if(t<m[i][j])
{
m[i][j]=t;//更新最优值
s[i][j]=k;//更新最优质的断点
}
}
}
}
}
int main() {
int n, ra, ca;
int p[MAX_MATRICES + 1];
int **s;
int **m;
scanf("%d", &n);
m = (int **) malloc((n + 1) * sizeof(int *));
s = (int **) malloc((n + 1) * sizeof(int *));
for (int i = 0; i <= n; i++) {
m[i] = (int *) malloc((n + 1) * sizeof(int));
s[i] = (int *) malloc((n + 1) * sizeof(int));
for (int j = 0; j <= n; j++) {
m[i][j] = INT16_MAX;//初始化一个较大的值
s[i][j] = 0;//断开点初始化为0
}
}
for (int i = 1; i <= n; i++) {
scanf("%d %d", &ra, &ca);
p[i - 1] = ra;//i从1开始,但数组从0开始
p[i] = ca;
}
// for(int i=0;i<=n;i++)
// printf("%d ",p[i]);
//调用函数
MatrixChain(m, s, p, n);
//输出最优解
printf("%d", m[1][n]);
// 释放分配的内存
for (int i = 0; i <= n; ++i) {
free(m[i]);
free(s[i]);
}
free(m);
free(s);
return 0;
}