P3195-[HNOI2008]玩具装箱【斜率优化dp】

正题

题目链接:https://www.luogu.com.cn/problem/P3195


题目大意

nnn个物品,分成若干段,每一段的长度为j−i+∑i=lrCkj-i+\sum_{i=l}^rC_kji+i=lrCk,打包价格为(长度−L)2(长度-L)^2(L)2

求最小价格和。


解题思路

si=∑j=1iCjs_i=\sum_{j=1}^iC_jsi=j=1iCj
fif_ifi表示iii前面的都打包完了,有
fi=fj+(si−sj+i−j−1−L)2f_i=f_j+(s_i-s_j+i-j-1-L)^2fi=fj+(sisj+ij1L)2
fi=fj+(si+i−sj−j−1−L)2f_i=f_j+(s_i+i-s_j-j-1-L)^2fi=fj+(si+isjj1L)2
然后dpdpdp可以做到O(n2)O(n^2)O(n2),因为有平方,考虑斜率优化
定义ai=si+i,bi=sj+j+1+La_i=s_i+i,b_i=s_j+j+1+Lai=si+i,bi=sj+j+1+L

fi=fj+(ai−bj)2f_i=f_{j}+(a_i-b_j)^2fi=fj+(aibj)2
fi=fj+ai2−2aibj+bj2f_i=f_j+a_i^2-2a_ib_j+b_j^2fi=fj+ai22aibj+bj2
2aibj+fi−ai2=fj+bj22a_ib_j+f_i-a_i^2=f_j+b_j^22aibj+fiai2=fj+bj2

对于每一个jjj表示一个点(bj,fj+bj2)(b_j,f_j+b_j^2)(bj,fj+bj2)(后文中称之为决策点)

考虑如何使fif_ifi最小,因为ai2a_{i}^2ai2是定值即让fi−ai2f_{i}-a_{i}^2fiai2最小,那么问题就变为了一条y=2aix+ky=2a_ix+ky=2aix+k的直线,要求经过某个决策点使得kkk最小。

那么显然,可能的点一定是一个下凸壳(相邻的点斜率单调上升),而因为2ai2a_i2ai这个斜率也是单调上升的,我们可以知道答案就是第一个决策点满足与下一个决策点的斜率≥2ai\geq 2a_i2ai

那么我们维护一个单调队列即可。

时间复杂度O(n)O(n)O(n)


codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#define pow2(x) ((x)*(x))
using namespace std;
const int N=5e4+10;
struct node{double x,y;int num;
}q[N];
int n,head,tail,L;
double s[N],a[N],b[N],f[N];
double slope(node x,node y)
{return ((y.y-x.y)/(y.x-x.x));}
int main()
{scanf("%d%d",&n,&L);for(int i=1;i<=n;i++){scanf("%lf",&s[i]);s[i]+=s[i-1];a[i]=s[i]+i;b[i]=s[i]+i+L+1;}b[0]=L+1;head=tail=1;q[1]=(node){b[0],b[0]*b[0],0};for(int i=1;i<=n;i++){while(head<tail&&slope(q[head],q[head+1])<2*a[i])head++;int p=q[head].num;f[i]=f[p]+pow2(a[i]-b[p]);node w=(node){b[i],f[i]+b[i]*b[i],i};while(head<tail&&slope(w,q[tail-1])<slope(q[tail-1],q[tail]))tail--;q[++tail]=w;}printf("%.0lf",f[n]);
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/321405.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【动态规划】书的复制 (ssl 1203)

书的复制书的复制书的复制 Description 现在要把m本有顺序的书分给k个人复制&#xff08;抄写&#xff09;&#xff0c;每个人的抄写速度都一样&#xff0c;一本书不允许分给两个或两个以上的人抄写&#xff0c;分给每个人的书&#xff0c;必须是连续的&#xff0c;比如不能把…

P1975-[国家集训队]排队【树状数组套线段树】

正题 题目链接:https://www.luogu.com.cn/problem/P1975 题目大意 一个nnn个数字的序列&#xff0c;每次交换两个数&#xff0c;求逆序对。 解题思路 显然交换lll和rrr的话对[1..l−1][1..l-1][1..l−1]和[r1,n][r1,n][r1,n]是不会有影响的&#xff0c;所以我们只需要考虑ll…

【动态规划】分组背包 (ssl 2291)

分组背包分组背包分组背包 Description 有N件物品和一个容量为V的背包。第i件物品的费用是c[i]&#xff0c;价值是w[i]。这些物品被划分为若干组&#xff0c;每组中的物品互相冲突&#xff0c;最多选一件。求解将哪些物品装入背包可使这些物品的费用总和不超过背包容量&#…

如何在Visual Studio 2017中使用C# 7+语法

前言之前不知看过哪位前辈的博文有点印象C# 7控制台开始支持执行异步方法&#xff0c;然后闲来无事&#xff0c;搞着&#xff0c;搞着没搞出来&#xff0c;然后就写了这篇博文&#xff0c;不喜勿喷&#xff0c;或许对您有帮助。在Visual Studio 2017配置支持C# 7语法心想都VS20…

11、mybatis的功能架构分析

1、Mybatis功能架构 1&#xff09;API接口层&#xff1a;提供给外部使用的接口API&#xff0c;开发人员通过这些本地API来操纵数据库。接口层一接收到调用请求就会调用数据处理层来完成具体的数据处理。 2&#xff09;数据处理层&#xff1a;负责具体的SQL查找、SQL解析、SQL执…

我的ACM之路-写于南宁站后

先说说我的ACM之路吧。自大一下学期开学开始&#xff0c;受到OI大佬舍友的影响&#xff0c;遂入坑。开始啃白书和紫书&#xff0c;然而发现好难啊。意识到大概是自己算法和数据结构知识欠缺的太厉害了&#xff0c;于是去网上买了本数据结构的教材&#xff08;隐约记得是清华大学…

P1903-[国家集训队]数颜色/维护队列【带修莫队】

正题 题目链接:https://www.luogu.com.cn/problem/P1903 题目大意 要求支持两个操作 QLR:Q\ \ L\ \ R:Q L R:询问L,RL,RL,R之间有多少个不同的数RPCol:R\ \ P\ \ Col:R P Col:将PPP修改为ColColCol 解题思路 莫队中多加一个时间维度ttt表示前面有多少个修改操作&#xf…

【动态规划】方格取数 (ssl 1010)

方格取数方格取数方格取数 Description 设有N*N的方格图(N<10,我们将其中的某些方格中填入正整数,而其他的方格中则放入数字0。如下图所示&#xff08;见样例&#xff09;&#xff1a; 某人从图的左上角的A 点出发&#xff0c;可以向下行走&#xff0c;也可以向右走&#…

ASP.NET Core 使用UrlFirewall对请求进行过滤

一. 前言UrlFirewall 是一个开源、轻便的对http请求进行过滤的中间件&#xff0c;可使用在webapi或者网关&#xff08;比如Ocelot&#xff09;,由我本人编写&#xff0c;并且开源在github&#xff1a;https://github.com/stulzq/UrlFirewall 欢迎star.二.UrlFirewall 介绍UrlFi…

1、设计模式和原则总述

目录 一、什么是设计模式&#xff1f;有什么用&#xff1f; 二、设计模式 三、设计原则 一、什么是设计模式&#xff1f;有什么用&#xff1f; 设计模式是一套代码设计的经验总结&#xff0c;使用设计模式可以提高代码的重用性、可靠性&#xff0c;提交代码内聚&#xff0c;降…

P3600-随机数生成器【dp,数学期望】

正题 题目链接:https://www.luogu.com.cn/problem/P3600 题目大意 nnn个数的序列&#xff0c;每个数是[1..x][1..x][1..x]中的一个&#xff0c;有qqq个区间[l..r][l..r][l..r]&#xff0c;求所有区间最小值的最大值的期望。 解题思路 首先如果一个区间包含别的区间&#xff0…

Visual Studio 15.7预览版4改进Git、C++支持

对于即将到来的Visual Studio 2017 15.7&#xff0c;微软已经发布了多个新的预览版本。这些版本的变更很有限&#xff0c;似乎离正式发布不远了。通常&#xff0c;变更的涉及面很广——因此&#xff0c;不管是用什么语言&#xff0c;开发人员都可以看到一些好处。第一次&#x…

【动态规划】农田个数 (ssl 1633)

农田个数农田个数农田个数 Description 你的老家在河北农村。过年时&#xff0c;你回老家去拜年。你家有一片NM农田&#xff0c;将其看成一个NM的方格矩阵&#xff0c;有些方格是一片水域。你的农村伯伯听说你是学计算机的&#xff0c;给你出了一道题&#xff1a; 他问你…

VB程序逆向常用的函数

转自: http://www.cnblogs.com/bbdxf/p/3780187.html VB程序逆向常用的函数 1) 数据类型转换: a) __vbaI2Str 将一个字符串转为8 位&#xff08;1个字节&#xff09;的数值形式(范围在 0 至 255 之间) 或2 个字节的数值形式(范围在 -32,768 到 32,767 之间)。 b)__vbaI4St…

P2827-蚯蚓【队列】

前言 早年一直拿堆过不了&#xff0c;结果发现要用队列 正题 题目链接:https://www.luogu.com.cn/problem/P2827 题目大意 有nnn条蚯蚓&#xff0c;每次选取最长的一条&#xff0c;切成⌊x∗p⌋\lfloor x*p\rfloor⌊x∗p⌋和x−⌊x∗p⌋x-\lfloor x*p\rfloorx−⌊x∗p⌋的两…

.net core DI 注册 Lazy 类型

当我们在 .net core (2.1) 中运行如下代码注入 Lazy<T> 变量的时候&#xff1a;public AccountService(Lazy<IHttpContextAccessor> httpContextAccessor) { }可能会遇到这样的错误提示&#xff1a;InvalidOperationException: Unable to resolve service for type…

操作系统复习笔记 02-03 OS Structure 操作系统结构

02-03操作系统结构&#xff08;OS_Structure&#xff09;[]陈述&#xff1a;1.IO设备与CPU可并行运行。2.每一个设备控制器负责一个设备类型。3.每一个设备控制器有一个局部缓存。4.CPU通过局部缓存与主存交换数据。5.设备控制器通过引起中断通知CPU操作已完成。[]中断机制&…

【深搜】 棋盘 【NOIp普及组 2017 第三题】 (luogu 3956/ssl 2851)

棋盘棋盘棋盘 题目大意&#xff1a; 有一个M*M的棋盘&#xff0c;要从&#xff08;1&#xff0c;1&#xff09;到&#xff08;m,m&#xff09;&#xff0c;中间有n个有颜色的格子&#xff0c;只能踩在有颜色的格子上&#xff0c;跳到不同颜色的格子要花费1元&#xff0c;可以将…

P3332-[ZJOI2013]K大数查询【树套树】

正题 题目链接&#xff1a;https://www.luogu.com.cn/problem/P3332 题目大意 开始nnn个可以重复的集合&#xff0c;要求支持操作 1lrc:1\ l\ r\ c:1 l r c:将ccc加入集合l∼rl\sim rl∼r中2lrk:2\ l\ r\ k:2 l r k:查询l∼rl\sim rl∼r的并集中第kkk大的数 解题思路 此题考…

别跟我谈EF抵抗并发,敢问你到底会不会用EntityFramework

前言一直以来写的博文都是比较温婉型的博文&#xff0c;今天这篇博文算是一篇批判性博文&#xff0c;有问题欢迎探讨&#xff0c;如标题&#xff0c;你到底会不会用EntityFramework啊。你到底会不会用EntityFramework啊面试过三年至六年的同行&#xff0c;作为过面试者到如今作…