CUMTOJ:algorithm-沙子的质量(题解）

题目描述

设有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。 

输出

合并的最小代价

样例输入

4
1 3 5 2

样例输出

22

算法描述

沙子每一步合并的总代价等于到上一步为止所需代价的最小值加上本次合并的代价，具有最优子结构性质；同时一个中间状态下一步可到达多个后继状态，具有重叠子问题性质。本题可采用动态规划算法求解。

C++代码

#include<iostream>
using namespace std;
int main()
{
	int n;
	cin >> n;
	int* a = new int[n+1];
	for (int i = 1; i < n+1; i++)
		cin >> a[i];
	int** b = new int* [n+1];
	for (int i = 0; i < n + 1; i++)
		b[i] = new int[n + 1];
	for (int i = 1; i < n+1; i++)
		for (int j = i; j < n + 1; j++)
		{
			int temp = 0;
			for (int s = i; s <= j; s++)
				temp += a[s];
			b[i][j] = temp;//b[i][j]存放从i堆到j堆沙子的质量之和
		}
	int** p = new int* [n + 1];
	for (int i = 0; i < n + 1; i++)
		p[i] = new int[n + 1];
	for (int i = 1; i < n + 1; i++)
		p[i][i] = 0;
	for (int k = 1; k < n; k++)
	{
		for (int i = 1; i <= n - k; i++)
		{
			int temp = 10000000;
			for (int j = i; j < i + k; j++)
				if (p[i][j] + p[j + 1][i + k] < temp)
					temp = p[i][j] + p[j + 1][i + k];
			p[i][i + k] = temp+b[i][i+k];//p[i][j]存放将i到j的所有堆沙子合并所需要的总代价
		}                                //总代价等于前过程所需要的代价加上i到j堆沙子的质量和
	}
	cout << p[1][n] << endl;
	return 0;
}