矩阵链相乘(动态规划法)

问题分析

矩阵链相乘问题是一个经典的动态规划问题。给定一系列矩阵,目标是找到一种最优的乘法顺序,使得所有矩阵相乘所需的标量乘法次数最少。矩阵链相乘问题的关键在于利用动态规划来避免重复计算子问题。

算法设计

  • 定义子问题:设 𝑤[𝑖,𝑗]表示计算矩阵 𝐴𝑖 到 𝐴𝑗的最小乘法次数。
  • 递归关系:根据递归关系式:

    其中 𝑝 是矩阵的维度数组,𝑝𝑖表示第 𝑖个矩阵的行数,​ p_{i+1}表示第 𝑖个矩阵的列数。   
  •  w[i,j]  表示第 i个矩阵到第 j个矩阵,路径上的矩阵链的最小代价
  •  p_{i-1}p_{k}p_{j} 代表从最优子解向它的最优父解构造时需要的代价
  • 边界条件:当 𝑖=𝑗 时,𝑤[𝑖,𝑗]=0,因为单个矩阵不需要乘法。
  • 计算顺序:从小规模问题开始,逐步计算较大规模问题。

其他的情况以此类推.....

伪代码

function MatrixChainOrder(p):
    n = length(p) - 1
    let w be a 2D array of size n x n
    for i = 1 to n:
        w[i, i] = 0
    for l = 2 to n:  // l is the chain length
        for i = 1 to n - l + 1:
            j = i + l - 1
            w[i, j] = infinity
            for k = i to j - 1:
                q = w[i, k] + w[k + 1, j] + p[i - 1] * p[k] * p[j]
                if q < w[i, j]:
                    w[i, j] = q
    return w[1, n]

C实现

#include <stdio.h>
#include <limits.h>

void MatrixChainOrder(int p[], int n) {
    int w[n][n];
    int i, j, k, l, q;

    for (i = 1; i < n; i++)
        w[i][i] = 0;

    for (l = 2; l < n; l++) {
        for (i = 1; i < n - l + 1; i++) {
            j = i + l - 1;
            w[i][j] = INT_MAX;
            for (k = i; k <= j - 1; k++) {
                q = w[i][k] + w[k + 1][j] + p[i - 1] * p[k] * p[j];
                if (q < w[i][j])
                    w[i][j] = q;
            }
        }
    }

    printf("Minimum number of multiplications is %d\n", w[1][n - 1]);
}

int main() {
    int arr[] = {1, 2, 3, 4};
    int size = sizeof(arr) / sizeof(arr[0]);

    MatrixChainOrder(arr, size);

    return 0;
}

结果如下:

总结

时间复杂度:算法的时间复杂度描述了算法执行所需要的时间与输入规模之间的关系。对于矩阵链相乘问题的动态规划算法,其时间复杂度通常为 𝑂(𝑛3)O(n3),其中 𝑛n 表示矩阵链中矩阵的个数。这样的时间复杂度在实际应用中通常是可以接受的。

空间复杂度:算法的空间复杂度描述了算法所需的额外空间与输入规模之间的关系。对于矩阵链相乘问题的动态规划算法,通常需要创建一个二维数组来存储子问题的解,因此空间复杂度为 𝑂(𝑛2)O(n2)。空间复杂度较低,适用于大部分计算机内存空间充足的情况。

初来乍到,蠢蠢小白,欢迎各路大神批评指正!!

相关推荐

  1. 动态规划矩阵

    2024-06-06 23:02:01       17 阅读
  2. 动态规划

    2024-06-06 23:02:01       6 阅读
  3. 动态规划学习

    2024-06-06 23:02:01       4 阅读
  4. [动态规划]矩阵取数游戏

    2024-06-06 23:02:01       29 阅读

最近更新

  1. TCP协议是安全的吗?

    2024-06-06 23:02:01       16 阅读
  2. 阿里云服务器执行yum,一直下载docker-ce-stable失败

    2024-06-06 23:02:01       16 阅读
  3. 【Python教程】压缩PDF文件大小

    2024-06-06 23:02:01       15 阅读
  4. 通过文章id递归查询所有评论(xml)

    2024-06-06 23:02:01       18 阅读

热门阅读

  1. 力扣linkedlist

    2024-06-06 23:02:01       5 阅读
  2. 【高频】如何保证缓存和数据库一致

    2024-06-06 23:02:01       6 阅读
  3. git本地仓库与远程仓库关联

    2024-06-06 23:02:01       10 阅读
  4. 如何重新设置路由器密码

    2024-06-06 23:02:01       8 阅读
  5. 华为防火墙基础配置

    2024-06-06 23:02:01       10 阅读
  6. 整理好了!2024年最常见 20 道 Kafka面试题(九)

    2024-06-06 23:02:01       12 阅读
  7. git随记

    git随记

    2024-06-06 23:02:01      7 阅读
  8. C++位运算

    2024-06-06 23:02:01       9 阅读
  9. springboot 解耦、隔离、异步的原则以及实战

    2024-06-06 23:02:01       9 阅读
  10. 图论第三天

    2024-06-06 23:02:01       6 阅读