在Mars星球上,每个Mars人都随身佩带着一串能量项链。在项链上有N颗能量珠。能量珠是一颗有头标记与尾标记的珠子,这些标记对应着某个正整数。并且,对于相邻的两颗珠子,前一颗珠子的尾标记一定等于后一颗珠子的头标记。因为只有这样,通过吸盘(吸盘是Mars人吸收能量的一种器官)的作用,这两颗珠子才能聚合成一颗珠子,同时释放出可以被吸盘吸收的能量。如果前一颗能量珠的头标记为m,尾标记为r,后一颗能量珠的头标记为r,尾标记为n,则聚合后释放的能量为m*r*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*2*3=60。
这一串项链可以得到最优值的一个聚合顺序所释放的总能量为
((4⊕1)⊕2)⊕3)=10*2*3+10*3*5+10*5*10=710。
第一行是一个正整数N(4≤N≤100),表示项链上珠子的个数。第二行是N个用空格隔开的正整数,所有的数均不超过1000。第i个数为第i颗珠子的头标记(1≤i≤N),当i<N< span>时,第i颗珠子的尾标记应该等于第i+1颗珠子的头标记。第N颗珠子的尾标记应该等于第1颗珠子的头标记。
至于珠子的顺序,你可以这样确定:将项链放到桌面上,不要出现交叉,随意指定第一颗珠子,然后按顺时针方向确定其他珠子的顺序。
只有一行,是一个正整数E(E≤2.1*109),为一个最优聚合顺序所释放的总能量。
4
2 3 5 10
710
d[i][j]表示合并区间[i,j]得到的最大能量。
d[i][j] = max(d[i][k] + d[k+1][j] + a[i]*a[k+1]*a[j+1]) ,其中a[i]表示第i颗珠子的头标记(注意搞清楚这几个头尾的关系,差点被搞晕)
假设以1号珠子为起点,则最大能量为d[1][n]。
现在没说起点是哪个,最直接想到的方法是枚举。可是,跑一次的复杂度是O(n^3),算上枚举的话就是O(n^4),对于n=100,很有可能会TLE,得有个改进的方法。(天梯上面的数据没卡时间,我枚举法交上去居然A了)
(以下思路来源http://blog.sina.com.cn/s/blog_4c396f4301000bol.html)
假设有五颗珠子,编号分别是1 2 3 4 5,则以2为起点时的顺序是2 3 4 5 1。
两次计算中重复计算的量有[2,3],[3,4],[4,5],[2,4],[3,5],[2,5] (方括号中的数字代表编号)
我们可以考虑将数组扩大一倍,令a[n+1]=a[1],a[n+2]=a[2]...a[n+n]=a[n].(最后一个没用,因为是循环)
这样一i为起点的答案就是d[i][i+n-1],扫一遍取最大就行了。
这样就将复杂度降为O(n^3)了。
#include<iostream>
#include<cassert>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<string>
#include<iterator>
#include<cstdlib>
#include<vector>
#include<stack>
#include<map>
#include<set>
using namespace std;
#define rep(i,f,t) for(int i = (f),_end_=(t); i <= _end_; ++i)
#define rep2(i,f,t) for(int i = (f),_end_=(t); i < _end_; ++i)
#define dep(i,f,t) for(int i = (f),_end_=(t); i >= _end_; --i)
#define dep2(i,f,t) for(int i = (f),_end_=(t); i > _end_; --i)
#define clr(c, x) memset(c, x, sizeof(c) )
typedef long long int64;
const int INF = 0x5f5f5f5f;
const double eps = 1e-8;
//*****************************************************
int a[210];
int d[210][210];
int n;
int main()
{
int mv = 0;
scanf("%d",&n);
for(int i = 1; i <= n; ++i){
scanf("%d",a+i);
a[i+n] = a[i];
}
for(int i = n+n-1; i >= 1; --i)
{
for(int j = i+1; j < min(n+n,i+n); ++j)
{
d[i][j] = max(d[i][j-1]+a[i]*a[j]*a[j+1], d[i+1][j]+a[i]*a[i+1]*a[j+1]);
for(int k = i+1; k < j-1; ++k){
d[i][j] = max(d[i][j], d[i][k]+d[k+1][j]+a[i]*a[k+1]*a[j+1]);
}
}
}
int ans = 0;
for(int i = 1; i <= n; ++i)ans = max(ans,d[i][i+n-1]);
cout<<ans<<endl;
return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。