【题目链接】 http://acm.hdu.edu.cn/showproblem.php?pid=3507

【题目大意】

  将长度为n的数列分段,最小化每段和的平方和。

【题解】

  根据题目很容易得到dp[j]=min(dp[k]+(s[j]-s[k])2),因为是从前往后转移,且决策单调,因此在CDQ分治的同时进行分治DP即可。

【代码】

#include <cstdio>
typedef long long LL;
const int N=500005;
int n,M,t;
LL f[N],g[N],a[N],s[N],INF=1LL<<60;
void DP(int l,int r,int dl,int dr){
int m=(l+r)>>1,i,dm=0;
LL &ret=g[m]; ret=INF;
for(i=dl;i<=dr&&i<m;i++){
LL t=f[i]+(s[m]-s[i])*(s[m]-s[i])+M;
if(t<ret)ret=t,dm=i;
}if(l<m)DP(l,m-1,dl,dm);
if(r>m)DP(m+1,r,dm,dr);
}
void CDQ(int l,int r){
if(l==r)return;
int mid=(l+r)>>1;
CDQ(l,mid);
DP(mid+1,r,l,mid);
for(int i=r;i>mid;i--)if(g[i]<f[i])f[i]=g[i];
CDQ(mid+1,r);
}
int main(){
while(~scanf("%d%d",&n,&M)){
for(int i=1;i<=n;i++){scanf("%lld",&s[i]);s[i]+=s[i-1];}
for(int i=1;i<=n;i++)f[i]=s[i]*s[i]+M;CDQ(0,n);
printf("%lld\n",f[n]);
}return 0;
}

  

最新文章

  1. TemplateMethod(模块方法模式)
  2. OAF_开发系列15_实现OAF组件重用和继承(案例)
  3. SQL2005四个排名函数(row_number、rank、dense_rank和ntile)的比较
  4. &lt;转载&gt;NPOI Excel 单元格背景颜色对照表
  5. ICA
  6. 使用NSData处理数据
  7. ios网站,博客
  8. maven+jetty项目在tomcat部署
  9. qt动态更新界面的菜鸟代码,请指出
  10. CentOS 在同一窗口打开文件夹
  11. OC可点击的两种轮播图效果
  12. Odoo 10的Linux安装
  13. 【转】Spring Boot 构建应用——快速构建 Spring Boot 应用
  14. backbond Model方法(set)
  15. aop日志记录
  16. springboot(十二):springboot单元测试、打包部署
  17. TCP3次握手和4次挥手及其为什么
  18. PyCharm 2019 最新激活方式总结(最新最全最有效!!!
  19. Beta阶段第五篇Scrum冲刺博客-Day4
  20. postman 断点

热门文章

  1. Java格式化输出
  2. 防抖(Debouncing)和节流(Throttling)
  3. 使用verilog实现4选1数据选择器的几种方法
  4. 如何获取浏览器的DNS解析时间
  5. E=MC2 - 搜搜百科
  6. redis存储session配制方法
  7. 现在网页中流行的css3样式
  8. ORACLE case when then
  9. Java环境的配置
  10. BZOJ 4199: [Noi2015]品酒大会( 后缀数组 + 并查集 )