矩阵连乘问题C++实现
矩阵连乘问题C++
- 1.认真审阅题目,明确题目的已知条件和求解的目标;
- 2.问题建模
- 3.算法设计;
- 4.编码实现;
1.认真审阅题目,明确题目的已知条件和求解的目标;
给定n个矩阵{A1,A2,A3……,An},其中Ai与Ai+1(i = 1,2 ,3,4……n-1)是可乘的,加括号的方法表示矩阵连乘的次序,不同加括号的方法所对应的计算次序是不同的。
2.问题建模
【例4-2】三个矩阵 A1 A2 A3 连乘,用加括号的方法表示其计算次序。
3个矩阵相乘,其加括号的方法一共有两种,具体如下:
【例4-3】4个矩阵连乘,用加括号的方法表示其计算次序。
4个矩阵连乘,其加括号的方法共有5种,具体如下:
不同加括号的方法所对应的计算量也是不同的,甚至差别很大。由于在矩阵相乘的过程中,仅涉及加法和乘法两种基本运算,乘法耗时远远大于加法耗时,故采用矩阵连乘所需乘法的次数来对不同计算次序的计算量进行衡量。
【例4-4】三个矩阵 A1 ,A2 ,A3 的行列分别为10×100、100×5、5×50,求例4-2中的两种加括号方法所需要乘法的次数。
两种加括号方法所需要乘法的次数分别为
那么,矩阵连乘问题就是对于给定n个连乘的矩阵,找出一种加括号的方法,使得矩阵连乘的计算量最小。
容易想到的解决方法是穷举法,即对n个矩阵连乘的每一种加括号方法进行乘法次数的统计,从中找出最小的计算量所对应的加括号方法。这种方法的复杂性取决于加括号的方法的种数。对于n个矩阵连乘,其加括号的方法有多少种呢?
考查矩阵连乘,不管哪种加括号的方法,最终都归结为两部分结果矩阵相乘,这两部分从n个连乘矩阵中的哪个矩阵处分开呢?设可能从A k 和A k+1 处将n个矩阵分成两部分,其中k=1,2,…,n-1。令P(n)代表n个矩阵连乘不同的计算次序,即不同加括号的方式,则n个矩阵连乘加括号的方式可通过两步操作来实现:①分别完成对两部分加括号;②对所得的结果加括号。由此
解此递推方程可得P(n)实际上是Catalan数,即P(n)=C(n-1),其中 。故穷举法的复杂性非常高,是n的指数级函数,显然,该方法不可行。
3.算法设计;
采用自低向上的方法求解最优质的具体的求解步骤设计如下:
步骤1:确定合适的数据结构,采用二维数组m来存放各个子问题的最优值,二维数组来存放来存放各个子问题的最优策略,(如果s[i][j]=k,则最优加括号的 方法为(Ai……Ak)(Ak+1……,Aj),一维数组P
步骤2:初始化。令m[i][j] = 0 , s[i][j] = 0,其中 i = 1, 2, ……
步骤3:循环阶段
步骤3-1:按照递归关系式计算两个矩阵AiAi+1,相乘时的最优值,并将其存入m[i][i+1],同时将最优决策计入s[i][i+1],i = 1, 2, ……
步骤3-2:按照递归关系式计算3个矩阵AiAi+1Ai+2相乘时的最优值并将其存入m[i][i+2],同时最优决策计入s[1][i+2],i = 1, 2, ……, n – 2
以此类推
步骤3-(n-1):按照递归关系式计算n个矩阵A1,A2,……An相乘时的最优质的并将其存入m[1][n],同时将最优决策计入s[1][n]
至此,m[1][n]即为原问题的最优值
步骤4:根据二维数组s记录的最优决策信息来构造最优解
步骤4-1:递归构造A1,……,As[1][n]的最优解,直到包含一个矩阵结束
步骤4-2:递归构造As[1][n]+1……An的最优解,直到包含一个矩阵结束
步骤4-3:将步骤4-1和步骤4-2递归的结果加括号
4.编码实现;
```cpp
#include <iostream>
using namespace std;
void MatrixChain(int *p, int n, int * * m, int * * s)//用于求最优值
{//初始化 for(int i = 1; i <= n; i ++){m[i][i] = 0;s[i][i] = 0;}for(int r = 2; r <= n; r++){//不同规模的子问题 for(int i = 1; i <= n - r + 1; i++){//每一个规模为r的矩阵连乘序列的首矩阵Ai int j = i + r - 1;//每个规模为r的矩阵连乘序列的尾矩阵为Ajm[i][j] = m[i+1][j] + p[i-1] * p[i]*p[j];//决策为k=i的乘法次数s[i][j] = 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];if(t < m[i][j]){m[i][j] = t;s[i][j] = k;}} }}
}
void Traceback(int i, int j, int ** s){//s[i][j]记录了断开的位置,即计算A[i:j]的加括号方式为(A[i:s[i][j]]) * A [s[i][j] + 1: j])if(i == j) return;Traceback(i, s[i][j], s);//递归打印A[i:s[i][j]的加括号方式Traceback(s[i][j] + 1, j, s);//递归打印A[s[i][j]+1:j ]的加括号的方式cout << "A" << "[" << i << ":" << s[i][j] << "]" << "乘以" << "A""[" << (s[i][j] + 1) << ":" << j << "]"<< endl;
}
int main(){return 0;
}
## 5.算法分析;显然,语句int t = m[i][k] + m[k+1][j] + p [i-1= * p[k]*p[j]为算法MatriChain的基本语句,最坏的情况下该语句的执行次数为O(n的三次方)故散发的最坏时间复杂度为O(n的三次方)
9.总结。
矩阵连乘问题C++实现相关推荐
- 给定一个 n X n 的矩阵 和 R,求旋转 90 X R 度以后的矩阵
1. 问题 [例题]循环输入,每组数据先输入一个 n(n<=10),然后 n 行 n 列数据代表一个矩阵, 每个数据为一个整数,然后再给出一个数字 R,输出这个矩阵按照顺时针旋转 90 X R ...
- 顺时针打印矩阵 python
前言 题目:输入一个矩阵,按照从外向里以顺时针的顺序依次打印出每一个数字,例如,如果输入如下矩阵: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 则依次打印出数字1,2 ...
- 矩阵乘以其矩阵转置求导-数学
20210703 矩阵论 https://zhuanlan.zhihu.com/p/288541909?utm_source=wechat_session 矩阵运算法则 20210529 https: ...
- 矩阵的卷积核运算(一个简单小例子的讲解)深度学习
卷积运算:假设有一个卷积核h,就一般为3*3的矩阵: 有一个待处理矩阵A: h*A的计算过程分为三步 第一步,将卷积核翻转180°,也就是成为了 第二步,将卷积核h的中心对准x的第一个元素,然后对应元 ...
- LeetCode简单题之二进制矩阵中的特殊位置
题目 给你一个大小为 rows x cols 的矩阵 mat,其中 mat[i][j] 是 0 或 1,请返回 矩阵 mat 中特殊位置的数目 . 特殊位置 定义:如果 mat[i][j] == 1 ...
- LeetCode简单题之距离顺序排列矩阵单元格
题目 给定四个整数 row , cols , rCenter 和 cCenter .有一个 rows x cols 的矩阵,你在单元格上的坐标是 (rCenter, cCenter) . 返回矩阵中的 ...
- LeetCode简单题之判断矩阵经轮转后是否一致
题目 给你两个大小为 n x n 的二进制矩阵 mat 和 target .现 以 90 度顺时针轮转 矩阵 mat 中的元素 若干次 ,如果能够使 mat 与 target 一致,返回 true : ...
- LeetCode简单题之矩阵中的幸运数
题目 给你一个 m * n 的矩阵,矩阵中的数字 各不相同 .请你按 任意 顺序返回矩阵中的所有幸运数. 幸运数是指矩阵中满足同时下列两个条件的元素: 在同一行的所有元素中最小 在同一列的所有元素中最 ...
- LeetCode简单题之矩阵对角线元素的和
题目 给你一个正方形矩阵 mat,请你返回矩阵对角线元素的和. 请你返回在矩阵主对角线上的元素和副对角线上且不在主对角线上元素的和. 示例 1: 输入:mat = [[1,2,3], [4,5,6], ...
- LeetCode简单题之托普利茨矩阵
题目 给你一个 m x n 的矩阵 matrix .如果这个矩阵是托普利茨矩阵,返回 true :否则,返回 false . 如果矩阵上每一条由左上到右下的对角线上的元素都相同,那么这个矩阵是 托普利 ...
最新文章
- 微服务架构下的静态数据通用缓存机制!
- hbase shell
- 数据结构实验之二叉树四:还原二叉树
- 即时通讯 我穿上球鞋
- 客户端javascript对象的几何属性(获得大小及坐标)
- pt5 mysql预处理_技术分享 | MySQL 监控利器之 Pt-Stalk
- java关闭applet_java – Applet会自动关闭
- Teradata SQL tips
- android 爱加密 脱壳,简单尝试脱“爱加密”官网加固的DEX壳
- 指数退避和AIMD为什么都青睐数字2
- ntoskrnl.exe 占用80端口 apache 无法启动
- Oracle 利用储存过程插入循环插入大量数据方法
- 龙芯中科官方宣布,龙芯中科企业办公信息化平台全面完成国产化替代
- 上周工作总结及本周工作安排
- 2021年目标,我打算这样去实现
- Win7运行程序时弹出“windows防火墙安全警报“的解决方法
- iPhone4S安装Linux系统,Absinthe 2.0.4 官网Windows/Mac/Linux原版下载—iPhone4S和iPad2完美越狱工具...
- 大众点评团购退款的小伎俩,折射出全商金融的事实?
- Johnson 全源最短路径算法
- python crop