MCM(矩阵链乘法)
这是《算法导论》动态规划中的一个问题。问题简述如下:我们在求解矩阵相乘时通常会有一个最优括号方案来对矩阵进行顺序相乘,这样会减少大量的计算时间。
我们知道矩阵A、B相乘,只能是当矩阵A的列数等于矩阵B的行数时才能进行相乘,且假设A·B = C,A为p×q规模的矩阵,B为q×r的矩阵,则乘积C的规模为p×r的矩阵。计算A·B所需的时间由乘法次数决定,即pqr。
例如:有三个矩阵的规模分别为:A = 10×100,B = 100×5,C = 5×50。如果按顺序(AB)C计算,则需要10×100×5 + 10×5×50 = 7500次乘法,如果按顺序A(BC)计算,则需要100×5×50 + 10×100×50 = 75000次乘法...所以,按第一种顺序计算矩阵链乘积要比第二种快非常多。
矩阵链乘法问题可描述如下:给定n个矩阵的链<A1,A2,...,An>,矩阵Ai的规模为pi-1×pi(1 ≤ i ≤ n),求完全括号化方案,使得计算乘积A1A2...An所需标量乘法次数最少。
分析问题描述可知道,我们假设在矩阵链乘中找到一个合适的切割点进行括号化,可得到问题的最优括号方案,然后对切割点的两边的子问题用相同的方法独立求解,不断进行下去,则最终将得到原问题的最优解。这就是该问题的最优子结构。
构造递归公式:令m[i][j] 表示计算矩阵Ai...j所需标量乘法次数的最小值,那么,原问题的最优解 —— 计算A1...n所需的最低代价就是m[1][n]。
对于i = j 时,矩阵链只包含唯一的矩阵 Ai...i = Ai ,不需要做任何标量乘法运算。若 i < j,利用最优子结构来计算m[i][j]。如前面假设Ai,Ai + 1,...,Aj的最优括号化方案的分割点在矩阵的Ak和Ak + 1之间,其中i ≤ k < j。那么,m[i][j]就等于计算Ai...k和Ak + 1...j的代价加上两者相乘的代价的最小值。由于矩阵Ai的大小为pi×pi+1,易知Ak+1...j相乘的代价为pi·pk+1·pj+1次标量乘法运算。因此,得到:
m[i][j] = m[i][k] + m[k + 1][j] + p[i]·p[k + 1]·p[j + 1]
由于k值取法只有 j - i 种可能的取值,即k = i, i + 1, ..., j。检查所有可能的情况,找到最优,因此递归公式变为:
( 0, i = j;
m[i][j] = {
( min(m[i][k] + m[k + 1][j] + p[i]·p[k + 1]·p[j + 1]), i < j.
O(n^3)迭代解法:
#include <iostream>
#include <vector> class Solution {
public:
std::vector<std::vector<int> > MatrixChainOrder(std::vector<int>& p)
{
const int n = p.size() - 1;
std::vector<std::vector<int> > m(n, std::vector<int> (n)); for(int i = 0; i < n; i++)
{
m[i][i] = 0;
}
for(int l = 1; l < n; l++)
{
for(int i = 0; i < n - l; i++)
{
int j = i + l;
m[i][j] = INT_MAX - 1;
for(int k = i; k < j; k++)
{
int q = m[i][k] + m[k + 1][j] + p[i]*p[k + 1]*p[j + 1];
if(q < m[i][j])
{
m[i][j] = q;
}
}
}
} return m;
}
}; int main()
{
std::vector<int> p{30, 35, 15, 5, 10, 20, 25};
Solution solve;
std::vector<std::vector<int> > res = solve.MatrixChainOrder(p); for(int i = res.size() - 1; i >= 0; i--)
{
for(int j = 0; j <= i; j++)
std::cout << res[j][i] << ' ';
std::cout << std::endl;
} return 0;
}
递归版感觉有点问题...所以再研究一段时间。
参考资料:以上皆为《算法导论》 — 动态规划原文
MCM(矩阵链乘法)的更多相关文章
- POJ1260 Pearls(dp,矩阵链乘法)
题目链接. 题目大意: 给定一个n,和两个序列a[i], p[i]. a[i] 表示需要购买 i品质 的数量,p[i] i 等级的价格. 1.每个品质都会有不同的价格,价格依据品质上升而上升 2.买一 ...
- Algorithm --> 矩阵链乘法
动态规划--矩阵链乘法 1.矩阵乘法 Note:只有当矩阵A的列数与矩阵B的行数相等时A×B才有意义.一个m×r的矩阵A左乘一个r×n的矩阵B,会得到一个m×n的矩阵C. #include ...
- CODEVS 3546 矩阵链乘法
http://codevs.cn/problem/3546/ 题目 给定有n个要相乘的矩阵构成的序列(链)<A1,A2,A3,.......,An>,要计算乘积A1A2.....An.一组 ...
- 【动态规划】XMU 1029 矩阵链乘法
题目链接: http://acm.xmu.edu.cn/JudgeOnline/problem.php?id=1029 题目大意: 题同乘法难题.给n+1个数,头尾不能动,中间的数可取出,取出时代价是 ...
- 【CODEVS】3546 矩阵链乘法
[算法]区间DP [题解] 注意先输出右括号后输出左括号. f[i][i+x-1]=min(f[i][i+x-1],f[i][j]+f[j+1][i+x-1]+p[i]*p[j+1]*p[i+x]) ...
- 51 Nod 1013 3的幂的和 矩阵链乘法||逆元+快速幂
这道题我写了两种写法 一种利用逆元 a/b%mod=a*c%mod; (c是b的逆元)易得2的逆元就是5~~~04: 一种是矩阵快速幂 利用递推式得出结论 #include<cstdio> ...
- 【UVa-442】矩阵链乘——简单栈练习
题目描述: 输入n个矩阵的维度和一些矩阵链乘表达式,输出乘法的次数.如果乘法无法进行,输出error. Sample Input 9 A 50 10 B 10 20 C 20 5 D 30 35 E ...
- COJ 0016 20603矩阵链乘
传送门:http://oj.cnuschool.org.cn/oj/home/solution.htm?solutionID=35454 20603矩阵链乘 难度级别:B: 运行时间限制:1000ms ...
- (最大矩阵链乘)Matrix-chain product
Matrix-chain product. The following are some instances. a) <3, 5, 2, 1,10> b) < ...
随机推荐
- Navicat Premium 12安装、激活
Navicat Premium 12安装 Navicat Premium 12激活
- LoadRunner通过验证参数判断事物的成功与失败
if(atoi(lr_eval_string("{Param_DiscountID}")) > 0){ //lr_message("多机多酒:%s",lr ...
- Linux02——解压和压缩
压缩文件①(tar)(zip) tar -zcvf (压缩名字).tar.gz (要压缩的内容,中间可用空格表示) -c 创建一个新归档 -v 显示详细信息 -f 指定压缩后的文件名 -z 通过gzi ...
- java截取小数点后两位
String a = "123.3445776";int i = a.indexOf(".");System.out.println(a.substring(0 ...
- 【转】通过python调用jenkins 常用api操作
原:https://www.cnblogs.com/L-O-N/p/11608220.html # -*- coding: utf- -*- import jenkins class TestJenk ...
- springboot笔记-1.自动化配置的关键
最近发现看过的东西容易忘,但是写一遍之后印象倒是会深刻的多. 总所周知springboot极大的简化了java开发繁琐性,而其最大的优势应该就是自动化配置了.比如要使用redis,我们直接引入相关的包 ...
- netty(三)---NioEventLoop分析
问题 : NioEventLoop 作用到底是什么?是在哪里用到的? NioEventLoop 和我们开头创建的 ServerBootstrap 和 EventLoopGroup 是什么关系 ? Ni ...
- Springboot学习:底层依赖与自动配置的原理
springboot依赖的父项目 我们在创建springboot项目的时候,设置了一个父项目: 这个项目可以点进去,可以发现它依赖于另一个父项目 再次点进去,发现没有依赖父项目了 观察这个项目的pom ...
- 《梳理业务的三个难点》---创业学习---训练营第二课---HHR---
一,<开始学习> 1,融资的第一步:把业务一块一块的梳理清楚. 2,预热思考题: (1)投资人会问能做多大,天花板怎么算?你的答案可以得到大家的认同吗?(四种方法,直接引用,自顶向下,自底 ...
- idea中springboot静态资源及页面跳转问题
1,静态资源放在resources/static下,html页面放在resources/templates下 2,在html中引入静态资源时,不用带static(对于路径来说是透明的) 3, 配置ht ...