最优矩阵连乘

描述

   一个n*m矩阵由n行m列共n*m个数排列而成。两个矩阵A和B可以相乘当且仅当A的列数等于B的行数。一个N*M的矩阵乘以一个M*P的矩阵等于一个N*P的矩阵，运算量为nmp。
   矩阵乘法满足结合律，A*B*C可以表示成(A*B)*C或者是A*(B*C)，两者的运算量却不同。例如当A=2*3 B=3*4 C=4*5时，(A*B)*C=64而A*(B*C)=90。显然第一种顺序节省运算量。
   现在给出N个矩阵，并输入N+1个数，第i个矩阵是a[i-1]*a[i]

输入格式

第一行n(n<=100)
第二行n+1个数

输出格式

最优的运算量

测试样例1

输入

3 
2 3 4 5

输出

64

思路：

dp[i][j]表示i和j之间的最小代价。

dp[i][j] = min(dp[i][k] + dp[k][j] + a[i] * a[k] * a[j]); i < k < j;

先计算距离为2的，然后距离为3的，渐次增加。

#include <iostream>  
#include <cstdio>  
#include <cstring>  
using namespace std;  
  
const int MAXN = 110;
const int INF = 0x3f3f3f3f;
int n, a[MAXN], dp[MAXN][MAXN];

int main()  
{  
	scanf("%d", &n);
	for (int i = 0; i <= n; i++)
	{
		scanf("%d", &a[i]);
	}

	memset(dp, 0, sizeof(dp));
	for (int r = 2; r <= n; r++)
	{
		for (int i = 0; i  + r <= n; i++)
		{
			int j = i + r;
			dp[i][j] = INF;
			for (int k = i + 1; k < j; k++)
			{
				dp[i][j] = min(dp[i][j], dp[i][k] + dp[k][j] + a[i] * a[k] * a[j]);
			}
		}
	}

	printf("%d\n", dp[0][n]);

    return 0;  
}