[NOIp2006] 能量项链

Description

在(Mars)星球上，每个(Mars)人都随身佩带着一串能量项链。在项链上有(N)颗能量珠。能量珠是一颗有头标记与尾标记的珠子，这些标记对应着某个正整数。并且，对于相邻的两颗珠子，前一颗珠子的尾标记一定等于后一颗珠子的头标记。因为只有这样，通过吸盘（吸盘是(Mars)人吸收能量的一种器官）的作用，这两颗珠子才能聚合成一颗珠子，同时释放出可以被吸盘吸收的能量。如果前一颗能量珠的头标记为(m)，尾标记为(r)，后一颗能量珠的头标记为r，尾标记为(n)，则聚合后释放的能量为(m imes r imes n)（(Mars)单位），新产生的珠子的头标记为(m)，尾标记为(n)。 需要时，(Mars)人就用吸盘夹住相邻的两颗珠子，通过聚合得到能量，直到项链上只剩下一颗珠子为止。显然，不同的聚合顺序得到的总能量是不同的，请你设计一个聚合顺序，使一串项链释放出的总能量最大。 例如：设(N=4)，(4)颗珠子的头标记与尾标记依次为((2,3) (3,5) (5,10) (10,2))。我们用记号⊕表示两颗珠子的聚合操作，((j)⊕(k))表示第(j,k)两颗珠子聚合后所释放的能量。则第(4)、(1)两颗珠子聚合后释放的能量为： ((4)⊕(1))(=10 imes 2 imes 3=60)。 这一串项链可以得到最优值的一个聚合顺序所释放的总能量为： (((4)⊕(1))⊕(2))⊕(3)）=(10 imes 2 imes 3+10 imes 3 imes 5+10 imes 5 imes 10=710)。

Input

第一行是一个正整数(N(4≤N≤100))，表示项链上珠子的个数。第二行是(N)个用空格隔开的正整数，所有的数均不超过(1000)。第(i)个数为第(i)颗珠子的头标记((1≤i≤N))，当(i<N)时，第(i)颗珠子的尾标记应该等于第(i+1)颗珠子的头标记。第(N)颗珠子的尾标记应该等于第(1)颗珠子的头标记。

至于珠子的顺序，你可以这样确定：将项链放到桌面上，不要出现交叉，随意指定第一颗珠子，然后按顺时针方向确定其他珠子的顺序。

Output

一个正整数(E(E≤2.1 imes (10)^9))，为一个最优聚合顺序所释放的总能量。

Sample Input

4
2 3 5 10

Sample Output

710

题解

/*
这是一个区间合并DP
由于是圆环，所以我们的第一步是将圆展开成线
接下来，设f[i][j]表示合并闭区间[i,j]所获得的最大能量
接下来推一下状态
我们令待合并的两个区间为f[l1][r1],f[l2][r2]
很显然l2=r1+1
于是f[l1][r2]=max{f[l1][r1]+f[r1+1][r2]+a[l1]*a[r1+1]*a[r2+1]}
再考虑一下r1的取值范围
因为一个闭区间[l,r]最少包括一个元素
所以l1<=r1,r1+1<=r2
变形并合并得l1<=r1<r2
*/
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
using namespace std;

const int N=250;
int n,a[N],f[N][N];

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;++i) scanf("%d",&a[i]);
	int n2=n+n;//减少常数计算
	for(int i=1;i<=n;++i) a[i+n]=a[i];
	a[n2+1]=a[1];
	for(int i=n2-1;i>=1;--i)//i即l1
		for(int j=i+1;j<=min(i+n-1,n2);++j)
		//j即r2，此处有一个优化，因为最多只有n个珠子
			for(int k=i;k<j;++k)//i即r1
				f[i][j]=max(f[i][j],f[i][k]+f[k+1][j]+a[i]*a[k+1]*a[j+1]);
	int Ans=0;
	for(int i=1;i<=n;++i) Ans=max(Ans,f[i][i+n-1]);
	printf("%d
",Ans);
	return 0;
}

---

/*
思路二
f[][]数组的含义同上
我们算f[][]数组可以根据区间的长度
状态转移方程同上
*/
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
using namespace std;

const int N=250;
int n,a[N],f[N][N];

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;++i) scanf("%d",&a[i]);
	int n2=n+n;//减少常数计算
	for(int i=1;i<=n;++i) a[i+n]=a[i];
	a[n2+1]=a[1];
	for(int i=1;i<n;++i)
		for(int j=n2-i;j>=1;--j)
			for(int k=j;k<j+i;++k)
				f[j][j+i]=max(f[j][j+i],f[j][k]+f[k+1][j+i]+a[j]*a[k+1]*a[j+i+1]);
	int Ans=0;
	for(int i=1;i<=n;++i) Ans=max(Ans,f[i][i+n-1]);
	printf("%d
",Ans);
	return 0;
}