矩陣連乘問題_第1頁
矩陣連乘問題_第2頁
矩陣連乘問題_第3頁
矩陣連乘問題_第4頁
矩陣連乘問題_第5頁
已閱讀5頁,還剩6頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權,請進行舉報或認領

文檔簡介

} if(i==j) return0; intu=LookupChain(i,i)+LookupChain(i+1,j)+p[i-1]*p[i]*p[j]; s[i][j]=i; for(intk=i+1;k<j;k++) { intt=LookupChain(i,k)+LookupChain(k+1,j)+p[i-1]*p[k]*p[j]; if(t<u) { u=t; s[i][j]=k; } } m[i][j]=u; returnu;}intMemoizedMatrixChain(intn,int**m,int**s){ for(inti=1;i<=n;i++) for(intj=i;j<=n;j++)m[i][j]=0; //把?m[i][j]的?矩?陣ó的?上?三▂角?全?部?賦3值μ為a0。£returnLookupChain(1,n);}//動ˉ態(tài)?規(guī)?劃?法ぁ?voidMatrixChain(){ for(inti=1;i<=n;i++) { m[i][i]=0; } for(intr=2;r<=n;r++) { for(inti=1;i<=n-r+1;i++) {intj=r+i-1;m[i][j]=m[i][i]+m[i+1][j]+p[i-1]*p[i]*p[j];s[i][j]=i;for(intk=i+1;k<j;k++) {inttemp=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];if(temp<m[i][j]) {m[i][j]=temp; s[i][j]=k; }} } }}//直±接ó遞蘗歸é法ぁ?intRecurMatrixChain(inti,intj){ if(i==j){return0;} intu=RecurMatrixChain(i,i)+RecurMatrixChain(i+1,j)+p[i-1]*p[i]*p[j];//遞蘗歸é,?p[]為a維?數(shù)簓 s[i][j]=i;//記?錄?加ó括ぁ?號?的?位?置? for(intk=i+1;k<j;k++) { intt=RecurMatrixChain(i,k)+RecurMatrixChain(k+1,j)+p[i-1]*p[k]*p[j];//遞蘗歸é if(t<u){u=t;s[i][j]=k;}//判D斷?哪?個?值μ更ü小?,?選?取?哪?個? } returnu;}voidTraceback(inti,intj){if(i==j)//如?果?只?有瓺一?個?矩?陣ó,?則ò直±接ó輸?出?{cout<<"A"<<i;}elseif(i+1==j)//加ó括ぁ?號?輸?出?{cout<<"(A"<<i<<"A"<<j<<")";}else{cout<<"(";Traceback(i,s[i][j]);Traceback(s[i][j]+1,j);cout<<")";}}MatrixMultiply.cpp#include<iostream>#include"matrix.h"usingnamespacestd;//intn;//定¨義?一?個?全?局?變?量?intmain(){L: intchoice; cout<<"請輸入相乘的矩陣的個數(shù)n="; cin>>n; cout<<"請輸入第1個矩陣行數(shù)和第1個矩陣到第"<<n<<"個矩陣的列數(shù):"; for(inti=0;i<=n;i++) { cin>>p[i]; } cout<<endl; cout<<"請選擇矩陣連乘的算法:"<<endl; cout<<"1.備忘錄算法"<<endl; cout<<"2.直接遞歸算法"<<endl; cout<<"3.動態(tài)規(guī)劃算法"<<endl; cout<<"4.重新輸入矩陣"<<endl; cout<<"5.退出程序"<<endl;cout<<endl; cout<<"請輸入選擇的編號:"; cin>>choice; cout<<endl; while(choice!=5){ switch(choice){ case1:LookupChain(1,n); cout<<"動態(tài)規(guī)劃算法:"<<endl; cout<<"矩陣連乘的最優(yōu)值為:"<<m[1][n]<<endl; cout<<"矩陣連乘的最優(yōu)解為:"; Traceback(1,n); cout<<endl; break; case2: RecurMatrixChain(0,n); cout<<"動態(tài)規(guī)劃算法:"<<endl; cout<<"矩陣連乘的最優(yōu)值為:"<<m[1][n]<<endl; cout<<"矩陣連乘的最優(yōu)解為:"; Traceback(1,n); cout<<endl; break; case3: MatrixChain(); cout<<"動態(tài)規(guī)劃算法:"<<endl; cout<<"矩陣連乘的最優(yōu)值為:"<<m[1][n]<<endl; cout<<"矩陣連乘的最優(yōu)解為:"; Traceback(1,n); cout<<endl; break; case4: gotoL; break; case5: choice=4; break; default: break;} cout<<endl; cout<<"請選擇矩陣連乘的算法:"<<endl; cout<<"1.備忘錄算法"<<endl; cout<<"2.直接遞歸算法"<<endl

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論