算法导论D 沙子的质量(矩阵连乘换皮)

题目：

题目描述

设有N堆沙子排成一排，其编号为1，2，3，…，N（N< =300）。每堆沙子有一定的数量，可以用一个整数来描述，现在要将N堆沙子合并成为一堆，每次只能合并相邻的两堆，合并的代价为这两堆沙子的数量之和，合并后与这两堆沙子相邻的沙子将和新堆相邻，合并时由于选择的顺序不同，合并的总代价也不相同，如有4堆沙子分别为1 3 5 2我们可以先合并1、2堆，代价为4，得到4 5 2又合并1，2堆，代价为9，得到9 2，再合并得到11，总代价为4+9+11=24，如果第二步是先合并2，3堆，则代价为7，得到4 7，最后一次合并代价为11，总代价为4+7+11=22；问题是：找出一种合理的方法，使总的代价最小。输出最小代价。

输入

第一行一个数N表示沙子的堆数N。 第二行N个数，表示每堆沙子的质量。 a[i]< =1000。

输出

合并的最小代价。

样例输入 Copy

4

1 3 5 2

样例输出 Copy

22

解析：

矩阵连乘的换皮，关键点也是找到一个最佳的拆分位置k

合并从i到j的最小代价就是合并从i到k的最小代价加上合并从k+1到j的最小代价再加上把这两堆合起来所需要的代价

dp[i][j] = dp[i][k] + dp[k+1][j] + sum,其中sum时把这两堆合并起来的代价，也就是从i到j的沙子的总质量。

代码：

#include<iostream>
using namespace std;
const int N = 305;
int dp[N][N];
int a[N];
int main()
{
	int n;
	cin>>n;
	for(int i=1;i<=n;i++)
		cin>>a[i];
	for(int i=1;i<=n;i++)
		dp[i][i] = 0;
	for(int r=2;r<=n;r++)
	{
		for(int i=1;i<=n-r+1;i++)
		{
			int j = i+r-1;
			int sum=0;
			for(int k=i;k<=j;k++)
				sum+=a[k];
			dp[i][j] = dp[i][i] + dp[i+1][j] + sum;
			for(int k=i+1;k<j;k++)
			{
				int t = dp[i][k] + dp[k+1][j] + sum;
				if(t<dp[i][j])
					dp[i][j] = t;
			}
		}
	}
	cout<<dp[1][n];
}