矩阵连乘重叠和备忘录.doc_第1页
矩阵连乘重叠和备忘录.doc_第2页
矩阵连乘重叠和备忘录.doc_第3页
矩阵连乘重叠和备忘录.doc_第4页
免费预览已结束,剩余1页可下载查看

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

#include#include/预设宽度setw()函数using namespace std;int m1010,s1010;int p10=30,35,15,5,10,20,25,;/数据预设/函数申明void MatrixMultiply(int a10,int b10,int c10,int ra,int ca,int rb,int cb);void MatrixChain(int p,int n,int m10,int s10);void Traceback(int i,int j,int s10);void Output(int n,int m10,int s10);int RecurMatrixChain(int i,int j);int MemoizedMatrixChain(int n,int m10,int s10);int LookupChain(int i,int j);int main() int ch,c=1;/ch用来标记排序方法的序号 int n=6;/设置为6个矩阵连乘 while(c!=0) coutendl; cout选择解决最优计算次序的方法: endl; cout- 1 动态规划endl; cout- 2 递归直接endl; cout- 3 备忘录endl; coutendlch; coutendl; switch(ch) case 1:MatrixChain(p,n,m,s);/动态规划 Output(n,m,s);/mij和sij Traceback(1,n,s);/最优计算次序 break; case 2:RecurMatrixChain(1,n);/递归直接 /Output(n,m,s); Traceback(1,n,s); break; case 3:MemoizedMatrixChain(n,m,s);/备忘录 /Output(n,m,s); Traceback(1,n,s); break; default: cout选择范围是13,重新输入:endl; return 0;void MatrixMultiply(int a10,int b10,int c10,int ra,int ca,int rb,int cb)if(ca!=rb)/Ai的列Ai+1的行cout矩阵不可乘!;for(int i=0;ira;i+)for(int j=0;jcb;j+)int sum=ai0*b0j;for(int k=1;kca;k+)sum+=aik*bkj;cij=sum;/动态规划void MatrixChain(int p,int n,int m10,int s10)for(int i=1;i=n;i+)mii=0;for(int r=2;r=n;r+)for(int i=1;i=n-r+1;i+)int j=i+r-1;mij=mi+1j+pi-1*pi*pj;sij=i;for(int k=i+1;kj;k+)int t=mik+mk+1j+pi-1*pk*pj;if(tmij)mij=t;sij=k;void Traceback(int i,int j,int s10)/最优计算次序if(i=j)return;/单个矩阵Traceback(i,sij,s);Traceback(sij+1,j,s);coutMultiply Ai,sij;cout and A(sij+1),jendl;/递归式直接计算int RecurMatrixChain(int i,int j)if(i=j)return 0;int u=RecurMatrixChain(i,i)+RecurMatrixChain(i+1,j)+pi-1*pi*pj;sij=i;for(int k=i+1;kj;k+)int t=RecurMatrixChain(i,k)+RecurMatrixChain(k+1,j)+pi-1*pk*pj;if(tu)u=t;sij=k;return u;/备忘录int MemoizedMatrixChain(int n,int m10,int s10)for(int i=1;i=n;i+)for(int j=i;j0)return mij;if(i=j)return 0;int u=LookupChain(i,i)+LookupChain(i+1,j)+pi-1*pi*pj;sij=i;for(int k=i+1;kj;k+)int t=LookupChain(i,k)+LookupChain(k+1,j)+pi-1*pk*pj;if(tu)u=t;sij=k;mij=u;return u;void Output(int n,int m10,int s10)for(int i=1;i=n;i+)for(int j=1;j=n;j+)coutsetw

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论