（POJ - 1651）Multiplication Puzzle（区间DP）

题目链接：1651 -- Multiplication Puzzle

题意：给出一组N个数，每次从中抽出一个数（第一和最后一个不能抽），该次的得分即为抽出的数与相邻两个数的乘积。直到只剩下首尾两个数为止。问最小得分是多少？

分析：这道题显然是一道区间dp题，设dp[i][j]表示把从i到j的数全部抽光所能得到的最小得分 ，由区间dp的模板我们知道，现在已经遍历到了区间长度为j-i+1的地方了，所以区间长度小于j-i+1的都已经更新完成，我们就可以从i到j枚举最后一个抽的数来更新dp[i][j]，也就是：

for(int len=1;len<=n-2;len++)
    {
        for(int i=2;i+len-1<=n-1;i++)
        {
            int j=i+len-1;
            for(int k=i;k<=j;k++)//枚举中间值
                dp[i][j]=min(dp[i][j],dp[i][k-1]+dp[k+1][j]+a[k]*a[i-1]*a[j+1]);
        }
    }

这个代码应该不难理解，下面我再针对这道题说一下一些细节：我们知道一开始dp应该初始化为正无穷大，但是我们在写转移方程时发现会出现i大于k-1的情况，如果这个值也被初始化为正无穷大的话那么动态转移就无法继续进行，我们就需要进行分类讨论，但是如果我们一开始令有效值为正无穷大，令无效值为0那我们就避免了这些分类讨论，有效值是指dp[i][j](i<=j)

下面分别附上分类讨论和不分类讨论的代码（区别就在于初始化，分类讨论不会用到无效值，所以无效值初始化就可以忽略不计，而不分类讨论就需要先将无效值初始化为0，防止对状态转移造成影响）

分类讨论代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>
#include<map>
#include<cmath>
#include<queue>
using namespace std;
const int N=103;
long long a[N],dp[N][N];//dp[i][j]表示把从i到j的数全部抽光所能得到的最小得分 
int main()
{
	int n;
	cin>>n;
	for(int i=1;i<=n;i++)
		scanf("%lld",&a[i]);
	memset(dp,0x3f,sizeof dp);
	for(int i=2;i<=n-1;i++)
		dp[i][i]=a[i]*a[i-1]*a[i+1];
	for(int len=2;len<=n-2;len++)
	{
		for(int i=2;i+len-1<=n-1;i++)
		{
			int j=i+len-1;
			for(int k=i;k<=j;k++)//枚举中间值 
			{
				if(k==i) dp[i][j]=min(dp[i][j],dp[k+1][j]+a[k]*a[i-1]*a[j+1]);
				else if(k==j) dp[i][j]=min(dp[i][j],dp[i][k-1]+a[k]*a[i-1]*a[j+1]);
				else dp[i][j]=min(dp[i][j],dp[i][k-1]+dp[k+1][j]+a[k]*a[i-1]*a[j+1]);
			}
		}
	}
	printf("%lld",dp[2][n-1]);
	return 0;
}

不分类讨论代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>
#include<map>
#include<cmath>
#include<queue>
using namespace std;
const int N=103;
long long a[N],dp[N][N];//dp[i][j]表示把从i到j的数全部抽光所能得到的最小得分 
int main()
{
	int n;
	cin>>n;
	for(int i=1;i<=n;i++)
		scanf("%lld",&a[i]);
	for(int i=2;i<=n-1;i++)//初始化 
	for(int j=i;j<=n-1;j++)
		dp[i][j]=0x3f3f3f3f3f3f3f3f;
	for(int len=1;len<=n-2;len++)
	{
		for(int i=2;i+len-1<=n-1;i++)
		{
			int j=i+len-1;
			for(int k=i;k<=j;k++)//枚举中间值
				dp[i][j]=min(dp[i][j],dp[i][k-1]+dp[k+1][j]+a[k]*a[i-1]*a[j+1]);
		}
	}
	printf("%lld",dp[2][n-1]);
	return 0;
}