矩阵连乘问题 动态规划

题目：给定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的最优值

比如有六个矩阵

思路步骤为：

1. 分析最优解的结构
2. 建立递归关系
3. 计算最优值

以下是本人手写推算过程

详细代码（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;
}