动态规划之矩阵连乘问题

一、矩阵连乘

给定n个矩阵｛A1,A2,…,An｝，其中，Ai与Ai+1是可乘的，(i=1,2 ,…,n-1)。用加括号的方法表示矩阵连乘的次序，不同的计算次序计算量（乘法次数）是不同的，找出一种加括号的方法，使得矩阵连乘的次数最小。

要求：
输入 矩阵数，各矩阵行数和列数P(p0,p1,…pn)
输出 矩阵连乘的最优值和最优解

二、代码参考

#include<iostream>
using namespace std;
const int N = 100;
int p[N];
int m[N][N];
int s[N][N];
void MatrixChain(int n) {
	int r,i, j, k;
	for (i = 0; i <= n; i++){
		m[i][i] = 0;
	}
	for (r = 2; r <= n; r++)
	{
		for (i = 1; i <= n - r + 1; i++)
		{
			j = i + r - 1;
			m[i][j] = m[i][i] + m[i + 1][j] + p[i - 1] * p[i] * p[j];
			s[i][j] = i;
			for (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) {
	if (i == j)
	{
		cout << "A" << i << " ";
		return;
	}
	cout << "(";
	Traceback(i, s[i][j]);
	Traceback(s[i][j] + 1, j);
	cout << ")";
}
int main()
{
	int n;
	cout << "请输入矩阵数n:";
	cin >> n;
	int i, j;
	cout << "请输入p(p0,p1,...pn):";
	for (i = 0; i <= n; i++)
	{
		cin >> p[i];
	}
	MatrixChain(n);
	cout << "矩阵连乘最优解为："<<endl;
	Traceback(1, n);
	cout << "\n矩阵连乘最优值为：" << endl;
	cout << m[1][n] << endl;
	return 0;
}

三、运行示例

输入
矩阵行 6
各矩阵行数和列数 30 35 15 5 10 20 25

输出如下