算法设计与分析 动态规划 3.1矩阵连乘问题（2022.03.18）

实现第五版，王晓东编著的第三章3.1矩阵连乘问题的代码。课本46页。

#include<iostream>
using namespace std;

void MatrixChain(int *p,int n,int **m,int **s);
void Traceback(int i,int j,int **s);
int main(){
	int n;
	cout<<"请输入矩阵个数:";
	scanf("%d",&n);
	int **m=new int *[n];//动态二维数组的创建，好东西，好好背
	int **s=new int *[n];
	for(int i=0;i<n+1;i++){
		m[i]=new int[i];//二维数组的初始化
		s[i]=new int[i]; 
	} 
	int *p=new int [n+1];//一维动态数组
	cout<<"请输入矩阵参数:"; 
	for(int i=0;i<=n;i++)
	scanf("%d",&p[i]);
	MatrixChain(p,n,m,s);
	cout<<"最小计算量为："<<m[1][n]<<endl;
	Traceback(1,n,s);
	return 0;
} 
void MatrixChain(int *p,int n,int **m,int **s){
	for(int i=1;i<=n;i++) m[i][i]=0;
	for(int r=2;r<=n;r++){
		for(int i=1;i<=n-r+1;i++){
			int j=i+r-1;
			m[i][j]=m[i+1][j]+p[i-1]*p[i]*p[j];
			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){
	if(i==j) return;
	Traceback(i,s[i][j],s);
	Traceback(s[i][j]+1,j,s);
	cout<<"Multiply A"<<i<<","<<s[i][j]<<" ";
	cout<<"and A"<<(s[i][j]+1)<<","<<j<<endl; 
}

其实也就是具体化了主函数，学会了动态数组的建立，两个void函数都是课本上的具体代码