HDU 3824/ BZOJ 3963 [WF2011]MachineWorks (斜率优化DP+CDQ分治维护凸包)
题面
BZOJ传送门(中文题面但是权限题)
HDU传送门(英文题面)
分析
定义f[i]f[i]f[i]表示在iii时间(离散化之后)卖出手上的机器的最大收益.转移方程式比较好写f[i]=max{f[j]−p[j]+r[j]+(d[i]−d[j]−1)∗g[j]}f[i]=max\{f[j]-p[j]+r[j]+(d[i]-d[j]-1)*g[j]\}f[i]=max{f[j]−p[j]+r[j]+(d[i]−d[j]−1)∗g[j]}
显然可以斜率优化,移项之后得到(f[j]−p[j]+r[j]−d[j]∗g[j]−g[j])=(−d[i]∗g[j])+(f[i])(f[j]-p[j]+r[j]-d[j]*g[j]-g[j])=(-d[i]*g[j])+(f[i])(f[j]−p[j]+r[j]−d[j]∗g[j]−g[j])=(−d[i]∗g[j])+(f[i])
也就是y=kx+by=kx+by=kx+b的形式,我们要让f[i]f[i]f[i]最大,也就是截距最大.那么维护一个上凸包就行了.为了求凸包插入点的xxx坐标单增,于是用CDQCDQCDQ分治来转移,转移时用归并排序.
- 啊啊啊啊啊啊啊啊啊啊啊啊啊 把whilewhilewhile写成ififif然后WA爆了
CODE
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
template<typename T>inline void read(T &num) {
char ch; int flg=1;
while((ch=getchar())<'0'||ch>'9')if(ch=='-')flg=-flg;
for(num=0;ch>='0'&&ch<='9';num=num*10+ch-'0',ch=getchar());
num*=flg;
}
const int MAXN = 100005;
const double eps = 1e-10;
struct Vector {
LL x, y;
Vector(){}
Vector(LL _x, LL _y):x(_x), y(_y){}
inline Vector operator -(const Vector &o) { return Vector(x-o.x, y-o.y); }
inline double operator *(const Vector &o) { return (double)x*o.y - (double)y*o.x; } //会炸longlong,只有用double
}p[MAXN], Q[MAXN];
LL f[MAXN];
struct Node {
LL p, r, d, g, x, y; int id;
inline void init() {
read(d), read(p), read(r), read(g);
x = g, y = -p+r-d*g-g;
}
}a[MAXN], tmp[MAXN];
inline bool cmpD(const Node &A, const Node &B) { return A.d < B.d; }
inline bool cmp(const Node &A, const Node &B) {
return A.x == B.x ? A.y + f[A.id] < B.y + f[B.id] : A.x < B.x;
}
inline int dcmp(double x) {
return x < -eps ? -1 : x < eps ? 0 : 1;
}
inline LL line(const Vector &P, LL k) { return P.y - k*P.x; } //求截距
void CDQ(int l, int r) {
if(l == r) { f[l] = max(f[l], f[l-1]); return; }
int mid = (l + r) >> 1;
CDQ(l, mid);
int n = 0, s = 0, t = 0;
for(int i = l; i <= mid; ++i) if(f[a[i].id] >= a[i].p)
p[++n] = Vector(a[i].x, a[i].y + f[a[i].id]);
for(int i = 1; i <= n; ++i) {
while(s+1 < t && dcmp((Q[t-1]-Q[t-2]) * (p[i]-Q[t-2])) >= 0) --t;
Q[t++] = p[i];
}
for(int i = mid+1; i <= r; ++i) {
while(s+1 < t && line(Q[s+1], -a[i].d) >= line(Q[s], -a[i].d)) ++s;
if(s < t) f[i] = max(f[i], line(Q[s], -a[i].d));
}
CDQ(mid+1, r);
int cur1 = l, cur2 = mid+1;
for(int i = l; i <= r; ++i) {
if(cur2 > r || (cur1 <= mid && cmp(a[cur1], a[cur2]))) tmp[i] = a[cur1++];
else tmp[i] = a[cur2++];
}
for(int i = l; i <= r; ++i) a[i] = tmp[i];
}
int n, D, kase;
int main () {
while(read(n), read(f[0]), read(D), n+f[0]+D) {
for(int i = 1; i <= n; ++i) a[i].init();
a[++n].d = D + 1;
sort(a + 1, a + n + 1, cmpD);
for(int i = 1; i <= n; ++i) a[i].id = i, f[i] = 0;
CDQ(1, n);
printf("Case %d: %lld\n", ++kase, f[n]);
}
}
//f[i] = f[j] - p[j] + r[j] + (D[i]-D[j]-1)*G[j]
//f[i] + (-D[i])*G[j] = f[j] - p[j] + r[j] - (D[j]+1)*G[j]
HDU 3824/ BZOJ 3963 [WF2011]MachineWorks (斜率优化DP+CDQ分治维护凸包)的更多相关文章
- BZOJ 3963: [WF2011]MachineWorks 斜率优化 + splay动态维护凸包
Description 你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理,公司使用更加先进的机械设备生产先进的机器.原来的那一台生产机器已经坏了,所以 ...
- 【bzoj3672】[Noi2014]购票 斜率优化dp+CDQ分治+树的点分治
题目描述 给出一棵以1为根的带边权有根树,对于每个根节点以外的点$v$,如果它与其某个祖先$a$的距离$d$不超过$l_v$,则可以花费$p_vd+q_v$的代价从$v$到$a$.问从每个点到1花费 ...
- BZOJ1492:[NOI2007]货币兑换 (CDQ分治+斜率优化DP | splay动态维护凸包)
BZOJ1492:[NOI2007]货币兑换 题目传送门 [问题描述] 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和B纪念券(以下简称B券).每个持有金券的 ...
- P4027 [NOI2007]货币兑换(斜率优化dp+cdq分治)
P4027 [NOI2007]货币兑换 显然,如果某一天要买券,一定是把钱全部花掉.否则不是最优(攒着干啥) 我们设$f[j]$为第$j$天时用户手上最多有多少钱 设$w$为花完钱买到的$B$券数 $ ...
- 【BZOJ2149】拆迁队(斜率优化DP+CDQ分治)
题目: 一个斜率优化+CDQ好题 BZOJ2149 分析: 先吐槽一下题意:保留房子反而要给赔偿金是什么鬼哦-- 第一问是一个经典问题.直接求原序列的最长上升子序列是错误的.比如\(\{1,2,2,3 ...
- BZOJ 1492 [NOI2007]货币兑换Cash:斜率优化dp + cdq分治
传送门 题意 初始时你有 $ s $ 元,接下来有 $ n $ 天. 在第 $ i $ 天,A券的价值为 $ A[i] $ ,B券的价值为 $ B[i] $ . 在第 $ i $ 天,你可以进行两种操 ...
- ●BZOJ 3963 [WF2011]MachineWorks
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3963 题解: 斜率优化DP,CDQ分治. 先按时间排序.(规定以下内容的第i台机器的卖出时间 ...
- BZOJ 3963: [WF2011]MachineWorks [CDQ分治 斜率优化DP]
传送门 当然了WF的题uva hdu上也有 你的公司获得了一个厂房N天的使用权和一笔启动资金,你打算在这N天里租借机器进行生产来获得收益.可以租借的机器有M台.每台机器有四个参数D,P,R,G.你可以 ...
- HDU 2993 MAX Average Problem(斜率优化DP)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2993 题目大意:给定一个长度为n(最长为10^5)的正整数序列,求出连续的最短为k的子序列平均值的最大 ...
随机推荐
- 登陆Linux服务器时触发邮件提醒
目前,客户只能在发现数据或者虚拟机被恶意侵入或者用户的误操作导致了数据的丢失之后,采取善后的手段,但是并没法做到提前的预警.那么通过 PAM 模块,就可以实现用户登录及获取root 权限时,通过邮件的 ...
- [转帖]Cacls和ICacls
Cacls和ICacls https://www.cnblogs.com/Aley/p/11089538.html Need Study 解释: Cacls:显示或修改文件的访问控制列表(ACL) ...
- C++之父给 C 程序员的建议
1. 在 C++中几乎不需要用宏, 用 const 或 enum 定义显式的常量, 用 inline 避免函数调用的额外开销,用模板去刻画一族函数或类型,用 namespace 去避免命名冲突. 2. ...
- 数据分析之--Mataplotlib入门
目录 Mataplotlib Seaborn 绘制线性图 图片的标题 点和线的样式 X和Y轴可读的映射 直方图 柱状图 条件性柱状图 饼图 箱图 散步图 3D图 Excel数据导入数据库 Matapl ...
- python中的 __inti__ 和 __new__ 方法的区别
这个要从Python的面向对象实例化的过程说起 类名() 之后,开辟一块内存空间,然后调用__init__把空间的内存地址作为self的参数传递到函数的内部,所有和self有关的参数,属性都会和sel ...
- Swarm系列7--存储介绍
存储介绍 1. 存储使用 与docker一样,在使用swarm服务级别的时候可以定义服务的存储需求, docker存储介绍参考: Docker之应用数据管理(volume/bind mount/tmp ...
- python __enter__ 与 __exit__的作用,以及与 with 语句的关系(转)
https://blog.csdn.net/xc_zhou/article/details/80810111 python __enter__ 与 __exit__的作用,以及与 with 语句的关系
- (九)springmvc之json的数据请求(客户端发送json数据到服务端)
index.jsp <%@ page language="java" contentType="text/html; charset=UTF-8" pag ...
- Python练习_函数进阶_day10
1. 1.作业 1,写函数,接收n个数字,求这些参数数字的和.(动态传参) 2,读代码,回答:代码中,打印出来的值a,b,c分别是什么?为什么? a=10 b=20 def test5(a,b): p ...
- Keras 笔记
1. 从 meta 模型恢复graph, 修改node 并保存 from __future__ import absolute_import from __future__ import div ...