矩阵连乘问题-动态规划-计算机算法

问题描述

问题分析

建立递归关系

代码实现

#include <iostream>
#include <algorithm>
#include <iomanip>
using namespace std;
#define MAX 0xfffffff
#define N 100
int n;
int p[N];
int s[N][N],dp[N][N];  //S存储切割位置，dp存储最优值 

void MatricChain(){
	for(int i=0;i<n+1;i++)//赋初值 
		for(int j=0;j<n+1;j++){
			dp[i][j]=MAX;
			s[i][j]=0;
		}
		
	for(int i=0;i<=n;i++){
		cin>>p[i];
		dp[i][i]=0;//只有一个矩阵时不能相乘
	}
 
	for(int L=2;L<=n;L++){//相乘矩阵的个数 
		for(int i=1;i<=n;i++){ 
			int j=i+L-1;
			if(j>n) break; 
			for(int k=i;k<j;k++){//遍历切割位置
				int min_=dp[i][j];
				int temp=dp[i][k]+dp[k+1][j]+p[i-1]*p[k]*p[j];
				if(temp<min_){
					dp[i][j]=temp;
					min_=temp;
					s[i][j]=k;
				}
			//	dp[i][j]=min(dp[i][j],dp[i][k]+dp[k+1][j]+p[i-1]*p[k]*p[j]);
			}
		}
	}
}

void Traceback(int i,int j){
	if(i==j) return;
	int k=s[i][j];
	Traceback(i,k);
	Traceback(k+1,j);
	cout<<"A["<<i<<":"<<k<<"]*A["<<k+1<<":"<<j<<"]"<<endl;
}

int main()
{
	cin>>n;
	
	MatricChain();
	
	
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++)
			cout<<setw(10)<<dp[i][j]<<" ";
		cout<<endl;
	}

	
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++)
			cout<<setw(5)<<s[i][j]<<" ";
		cout<<endl;
	}
	 
	cout<<"连乘的最少次数是"<<dp[1][n]<<"次。"<<endl;
	Traceback(1,n);
}
/*
6
2 7 5 4 2 3 8
5
30 35 15 5 10 20
*/

时间复杂度

O（n^3）

具体讲解

讲解