在《算法导论》第二版中,(P197),关于动态规划的算法,这个动态规划解决的是矩阵链乘法的加括号的问题。给定由n个要相乘的矩阵构成的序列,要计算乘积
最小的相乘次数。算法与《装配线调度》有一些不同点。每一个状态m[i][j]的值不是仅与前两个有关,而是与1~n的值都有关。所以前面所有的值在求最终的值时都会被用到。算法导论中的矩阵看起来真是难受,这个也是因人而异吧,倒过来的矩阵总是看起来怪怪的,还是自己画原始的图看起来舒服。但是作者的意图是为了让人感受到这里面递归和迭代的过程。
我们都是先求数组坐标 (0,1) 然后再求:(0,2) (0,3) (0,4) (0,5)
(1,2) (1,3) (1,4) (1,5)
(2,3) (2,4) (2,5)
(3,4) (3,5)
(4,5)
这样可以看出:每个行都是从零开始,每个列都是依次递增。所以这样可以组织代码:
#include "stdafx.h"
#include <iostream>
using namespace std;
int const NUM_MTRIX = 6; //矩阵的数量
void printArray( int arr[], int length )
{
for ( int i = 0; i < length; ++i )
{
cout<<arr[i]<<" ";
}
cout<<endl;
}
int _tmain(int argc, _TCHAR *argv[])
{
int P[7] = { 30, 35, 15, 5, 10, 20, 25 };
int min[NUM_MTRIX][NUM_MTRIX];
for ( int i = 0; i < NUM_MTRIX; ++i )
{
min[i][i] = 0; //将对角线元素全部赋值为0
}
for ( int col = 1; col < NUM_MTRIX; ++ col )
{
int row = 0;
for ( int tmpCol = col; tmpCol < NUM_MTRIX ; tmpCol++ )
{
min[row][tmpCol] = INT_MAX;
for ( int k = row; k < tmpCol; ++k)
{
int tmpMin = min[row][k] + min[ k +1 ][tmpCol] + P[ row ]* P[ k +1 ]*P[ tmpCol +1 ];
if ( tmpMin < min[row][tmpCol] )
min[row][tmpCol] = tmpMin;
}
++row;
}
}
cout<<min[0][5]<<endl;
}