3002 石子归并 3

题目描述 Description

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

输入描述 Input Description

第一行一个整数n(n<=3000)

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

输出描述 Output Description

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

样例输入 Sample Input

4

4 1 1 4

样例输出 Sample Output

18

数据范围及提示 Data Size & Hint

数据范围相比“石子归并” 扩大了

这个题和石子归并1唯一的区别就是数据范围变大了

于是用基本的做法写出来就TLE,只能得50分

这道题正确的解法是四边形不等式优化dp,为此初步了解四边形不等式优化方法

通俗的说,就是多了一个s[l][r]数组,用以记录得到l到r区间的最优解用的是哪个点作为断点

关于s[][]的正确性的证明我还没能弄懂,不过其原理很显然,在石子归并问题中其断点随区间向右移动,是有单调性的,因此有s(i,j-1)≤s(i,j)≤s(i+1,j)

所以在区间l到r内枚举断点时只需要枚举s[l][r-1]~s[l+1][r]之间的点

复杂度降为n^2

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=,inf=;
int f[maxn][maxn],s[maxn][maxn];
int n,a[maxn],w[maxn][maxn];
int main(){
scanf("%d",&n);
for(int i=;i<=n;i++){
scanf("%d",&a[i]);
w[i][i]=a[i];
}
for(int i=;i<=n;i++)
for(int j=i+;j<=n;j++)
w[i][j]=w[i][j-]+a[j];
for(int i=;i<=n;i++)s[i][i]=i;
for(int p=;p<n;p++){
for(int i=;i<=n-p;i++){
int j=i+p;
f[i][j]=inf;
for(int k=s[i][j-];k<=s[i+][j];k++){
if(f[i][j]>f[i][k]+f[k+][j]+w[i][j]){
f[i][j]=f[i][k]+f[k+][j]+w[i][j];
s[i][j]=k;
}
}
}
}
printf("%d",f[][n]);
}

最新文章

  1. Service and controller in angularJs
  2. svn 更新命令(冲突时使用theirs)
  3. IT公司100题-26-左旋转字符串
  4. 有序列表和无序列表、流、格式布局:position
  5. ubuntu下QT输出程序控制台界面难看的解决方法
  6. 【ELK_Log4net】.net Core重写一个TcpAppender
  7. easyui datagrid 右边框被隐藏
  8. 【转贴】一次 JDBC 与 MySQL 因 “CST” 时区协商误解导致时间差了 14 或 13 小时的排错经历
  9. 借助ssh隧道和中间主机,使本地主机可以直连远程主机
  10. excel vlookup简易样例【原】
  11. 使用JAVA数组实现顺序栈
  12. Redis中存入存储的编码方式不一致解决问题
  13. 转载:return *this和return this
  14. 序列化ADODataSet, ADOQuery
  15. 23种设计模式之命令模式(Command)
  16. shell打印彩色输出
  17. CSS常见简写规则整理
  18. Python爬虫教程-00-写在前面
  19. hadoop2.x常用端口、定义方法及默认端口、hadoop1.X端口对比
  20. 【luogu P1939 【模板】矩阵加速(数列)】 题解

热门文章

  1. 从动态库的def文件生成lib文件
  2. [STM8L]基于STM8L152的TAB段式LCD液晶驱动的分析 - 单片机干货 - 中国电子技术论坛 - 最好最受欢迎电子论坛!
  3. java入门了解05
  4. 51nod 1610
  5. CSS3实现自定义Checkbox动画
  6. ThreadPoolExecutor线程池进阶使用
  7. OpenCV——PS 滤镜算法之平面坐标到极坐标的变换
  8. Could not load the &quot;xxx.png&quot; image referenced from a nib in the bundle with identifier &quot;com.xxxx&quot;
  9. ActivityMq的使用(小例子)
  10. 每天一个Linux命令(2):ls命令