动态规划：石子归并(区间型动态规划)

题目描述 Description

有n堆石子排成一列，每堆石子有一个重量w[i], 每次合并可以合并相邻的两堆石子，一次合并的代价为两堆石子的重量和w[i]+w[i+1]。问安排怎样的合并顺序，能够使得总合并代价达到最小。

输入描述 Input Description

第一行一个整数n（n<=100）

第二行n个整数w1,w2...wn  (wi <= 100)

输出描述 Output Description

一个整数表示最小合并代价

样例输入 Sample Input

4

4 1 1 4

样例输出 Sample Output

18

算法分析：

首先要对问题进行划分，解决子问题往往更容易下手；

假设区间为（i，j）,设最优归并代价为：f（i，j）,设区间石子总和为sum（i，j）,那么

f( i , j ) = min( f( i , k ) + f( k+1, j) + sum( i , j ) ) 其中的  i <= k <= j，这是第一层分析为状态转移方程；

f( i , i) = 0，对于所有的区间长度为 0 的区间求解，其最优代价都是 0 ，

当f ( i , j ) 的区间长度为2时，即有f ( i , j ) = f ( i , i ) + f ( j , j ) + sum( i , j );

当f ( i , j ) 的区间长度为3时，即有f ( i , j ) = f ( i , i ) + f ( i , j ) + sum( i , j ) 和  f ( i , i+1 ) + f ( i+2 , j ) + sum( i , j )

和  f ( i , j ) + f ( j , j ) + sum( i , j ) 中的最小值；

以此类推...

代码实现（含注释）：

#include <iostream>   
#include <stdio.h>
using namespace std;
#define M 101
#define INF 1000000000
int n,f[M][M],sum[M][M],stone[M];
int main()
{
	int i,j,k,t;
	cin>>n;
	for(i=1;i<=n;i++)
		scanf("%d",&stone[i]);

	for(i=1;i<=n;i++)
	{
		f[i][i]=0;
		sum[i][i]=stone[i];//题目特点，每次归并都要加上归并区间中石子的总代价
		for(j=i+1;j<=n;j++)
			sum[i][j]=sum[i][j-1]+stone[j];
	}
	for(i=2;i <= n;i++)
	{
		for(j=1;j <= n;j++)
		{
			f[j][j+i-1] = INF;
			for(k=j;k < j+i-1;k++)
			{
				if(f[j][j+i-1] > f[j][k] + f[k+1][j+i-1] + sum[j][j+i-1])
                //关键的转态转移方程：区间的每一个划分子问题集中最优的就对应该问题的最优解
				f[j][j+i-1] = f[j][k] + f[k+1][j+i-1] + sum[j][j+i-1];
			}
		}
	}
	printf("%d\n",f[1][n]); 
	return 0;
}