矩阵相乘问题

矩阵相乘问题

ID:70615626

大小:611.10 KB

页数:4页

时间:2021-11-23

矩阵相乘问题_第1页
矩阵相乘问题_第2页
矩阵相乘问题_第3页
矩阵相乘问题_第4页
资源描述:

《矩阵相乘问题》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、符号三角形问题一、问题描述:给定n个矩阵{A1,A2,...,An},其中Ai与Ai+1是可乘的,i=1,2,...,n-1,求矩阵相乘至少需要多少次数乘?二、输入:A110*100,A2100*5,A35*50三、输出:图一为A1*A2*A3至少需要的次数乘数值三、算法描述:寻找最优子结构此问题最难的地方在于找到最优子结构。对乘积A1A2...An的任意加括号方法都会将序列在某个地方分成两部分,也就是最后一次乘法计算的地方,我们将这个位置记为k,也就是说首先计算A1...Ak和Ak+1...An,然后

2、再将这两部分的结果相乘。最优子结构如下:假设A1A2...An的一个最优加括号把乘积在Ak和Ak+1间分开,则前缀子链A1...Ak的加括号方式必定为A1...Ak的一个最优加括号,后缀子链同理。一开始并不知道k的确切位置,需要遍历所有位置以保证找到合适的k来分割乘积。一、算法设计:构造递归解设m[i,j]为矩阵链Ai...Aj的最优解的代价,则构建辅助表,解决重叠子问题从第二步的递归式可以发现解的过程中会有很多重叠子问题,可以用一个nXn维的辅助表m[n][n]s[n][n]分别表示最优乘积代价及其分

3、割位置k。辅助表s[n][n]可以由2种方法构造,一种是自底向上填表构建,该方法要求按照递增的方式逐步填写子问题的解,也就是先计算长度为2的所有矩阵链的解,然后计算长度3的矩阵链,直到长度n;另一种是自顶向下填表的备忘录法,该方法将表的每个元素初始化为某特殊值(本问题中可以将最优乘积代价设置为一极大值),以表示待计算,在递归的过程中逐个填入遇到的子问题的解。二、程序:#include<iostream>usingnamespacestd;constintN=4;//p为矩阵链

4、,p[0],p[1]代表第一个矩阵的行数和列数,p[1],p[2]代表第二个矩阵的行数和列数......length为p的长度//所以如果有六个矩阵,length=7,m为存储最优结果的二维矩阵,s为存储选择最优结果路线的//二维矩阵voidMatrixChainOrder(int*p,intm[N][N],ints[N][N],intlength){intn=length-1;intl,i,j,k,q=0;//m[i][i]只有一个矩阵,所以相乘次数为0,即m[i][i]=0;for(i=1;i&am

5、p;lt;length;i++){m[i][i]=0;}//l表示矩阵链的长度//l=2时,计算m[i,i+1],i=1,2,...,n-1(长度l=2的链的最小代价)for(l=2;l<=n;l++){for(i=1;i<=n-l+1;i++){j=i+l-1;//以i为起始位置,j为长度为l的链的末位,m[i][j]=0x7fffffff;//k从i到j-1,以k为位置划分for(k=i;k<=j-1;k++){q=m[i][k]+m[k+1][j]+p[i

6、-1]*p[k]*p[j];if(q<m[i][j]){m[i][j]=q;s[i][j]=k;}}}}cout<<m[1][N-1]<<endl;}voidPrintAnswer(ints[N][N],inti,intj){if(i==j){cout<<"A"<<i;}else{cout<<"(&

7、amp;quot;;PrintAnswer(s,i,s[i][j]);PrintAnswer(s,s[

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。