Input:输入两行标准输入,第一行表示矩阵的个数,第二行是这些矩阵的维数。
Output:输出两行标准输出,第一行为这些矩阵的最小数乘次数,第二行为这些矩阵加括号后的表达式。
输入样例1:
6
30 35 15 5 10 20 25
输出样例1:
15125
((A1(A2A3))((A4A5)A6))
具体使用的数据结构:
输入n表示矩阵的个数
int *p = new int[n+1]用来表示矩阵的维数。例第一个矩阵的维数是(p[0],p[1]),第i个矩阵的维数是(p[i-1],p[i])
二维数组m[i][j]用来表示第i个矩阵到第j个矩阵所需要的最小数乘次数
二维数组s[i][j]=k用来表示第i个矩阵到第j个矩阵在第k的矩阵后断开。
子问题的递归关系:
#include <iostream>
#include <stdio.h>
#include <stdlib.h>
using namespace std;
int m[101][101];//m[i][j]表示矩阵i到j需要的最小数乘次数
int s[101][101];//s[i][j]表示矩阵i到j从哪里断开
void MatrixChain(int n,int *p)
{
for(int i=1;i<=n;i++){
m[i][i] = 0;//第一次遍历,对角线元素全部置0
}
for(int r=1;r<n;r++)
for(int i=1;i<=n-r;i++){
int j=i+r;
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 temp = 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;}
}
}
}
//递归打印出最优解表达式
void PrintAnswer(int i,int j)
{
if(i==j) cout<<"A"<<i;
else{
cout<<"(";
PrintAnswer(i,s[i][j]);
PrintAnswer(s[i][j]+1,j);
cout<<")";
}
}
int main()
{
int n;//表示有n个矩阵相乘
cin>>n;
int *p = new int[n+1];//数组下标从0开始,n个矩阵有n+1个维度参数
for(int i=0;i<=n;i++){
cin>>p[i];//依次输入矩阵维度
}
MatrixChain(n,p);
cout<<m[1][n]<<endl;
PrintAnswer(1,n);
return 0;
}
运行结果:
总结: