/*
矩阵链乘法:
给定一个n个矩阵的序列<A1,A2,...,An>,矩阵Ai的规模为Pi-1*Pi(1<=i<=n),计算乘积A1*A2*...*An
设Ai,j表示AiAi+1...Aj乘积的结果矩阵
设m[i,j]表示计算矩阵Ai,j乘法次数最小值
那么我们需要求的结果就是m[1,n]
设AiAi+1...Aj的最优括号方案的分割点在矩阵Ak和Ak+1之间,其中i<=k<j。
那么m[i,j]等于计算Ai,k和Ak+1,j的代价加上两者相乘的代价的最小值
设矩阵Ai的大小为Pi-1*Pi,那么Ai,k与Ak+1,j相乘的代价是Pi-1PiPj
其中k属于[i,j-1]
当m[i,i]时,表示只有一个矩阵,所以相乘的次数为0
m[i,j]= { 0 , i = j
{ min i<=k<j { m[i,k] + m[k+1, j] + Pi-1*Pk*Pj } , i < j
我们用s[i,j]保存AiAi+1...Aj的最优分割点k
采用自底向上的方法,应该按照长度顺序求解矩阵链括号化问题
输入第一行的第一个数n为矩阵的个数,第二行有n+1个数,分别为n个矩阵各自的行数+最后第n个矩阵的列数
输出最小乘法次数
输入:
6
30 35 15 5 10 20 25
输出:
15125
((A1(A2A3))((A4A5)A6))
*/
#include <iostream>
#include <string.h>
using namespace std;
const int MAXSIZE = 10000;
int m[MAXSIZE][MAXSIZE]; //m[i][j]表示矩阵AiAi+1...Aj需要的最少乘法次数
int s[MAXSIZE][MAXSIZE];//保存最优分割点,s[1][n]计算了A1...An的最优分割点,
int p[MAXSIZE]; //保存矩阵的行数
int matrixChain(int n)
{
//初始化只有1个矩阵时的最小乘法次数为0
for(int i = 0 ; i < MAXSIZE ; i++)
{
m[i][i] = 0;
}
//逐渐将矩阵链长度从2累加到 n + 1 - l ,之所以 - l 是因为我们需要使j = i + l,
//将矩阵链长度从2累加到n
for(int l = 2 ; l <= n ; l++)
{
//i表示起始的矩阵下标,因为 终点下标 j = i + l - 1 <= n ,所以推得 i 的范围最多 = n - (l - 1) = n + 1 -l
for(int i = 1 ; i <= n - l + 1 ; i++)
{
int j = i + l - 1;//第一趟循环求m[i,i+1],第二趟求m[i,i+2]
int max = -1;
m[i][j] = INT_MAX;
for(int k = i ; k < j ; k++)
{
int value = m[i][k] + m[k+1][j] + p[i-1]*p[k]*p[j];
if(value < m[i][j])
{
m[i][j] = value;
s[i][j] = k;
}
}
}
}
return m[1][n];
}
void printResult(int i,int j)
{
if(i == j)
{
cout << "A" << i ;
}
else
{
cout << "(" ;
printResult(i , s[i][j]);//递归求取最优分割点A1...s[1,n] As[1,n]+1...n
printResult(s[i][j] + 1 , j);
cout << ")";
}
}
void process()
{
int n;
while(cin >> n)
{
for(int i = 0 ; i <= n ; i++)
{
cin >> p[i] ;
}
int iResult = matrixChain(n);
cout << iResult << endl;
printResult(1,n);
}
}
int main(int argc, char* argv[])
{
process();
getchar();
return 0;
}
算法导论:第15章 动态规划_2矩阵链乘法
最新推荐文章于 2021-02-22 21:58:11 发布