bzoj1492/luogu4027 货币兑换 (斜率优化+cdq分治)
设f[i]是第i天能获得的最大钱数,那么
f[i]=max{在第j天用f[j]的钱买,然后在第i天卖得到的钱,f[i-1]}
然后解一解方程什么的,设$x[j]=\frac{F[j]}{A[j]*Rate[j]+B[j]}$,$y[j]=Rate[j]*x[j]$的话,就能得到$f[i]=max\{y[j]*A[i]+x[j]*B[i],f[i-1]\}$
然后再推一波斜率优化的式子,就可以得到,当j1比j2优时,$\frac{y[j1]-y[j2]}{x[j1]-x[j2]}<-\frac{B[i]}{A[i]}$
然而$-\frac{B[i]}{A[i]}$这玩意并不单调,所以需要用splay或者cdq分治维护一个斜率递减的凸包,来查询满足上式的斜率最大的那个点。
然而懒得写splay了,所以用cdq分治来做(wa了几发以后才发现还不如去写splay呢写splay就不只是wa了)
我们来分治对所有点的询问。当然需要按照i从小到大来做。
然后在做一个区间的时候,我是想用它的左子区间(保证这些的值已经求完了)去更新右子区间的值。
也就是说,需要先做左子区间,然后做当前区间,最后再做右子区间(原理和寻找宝藏是一样的)
做的时候,左边按x值排序,右边按$-\frac{B[i]}{A[i]}$排序,然后给左半边用栈维护一个凸包,再用右半边的去扫,更新结果就行了。
(别忘了f[i]可以等于f[i-1],具体做法是在l=r的时候更新一下)
(常数好像写得很大..不管了.....)
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<set>
#include<ctime>
#define LL long long
using namespace std;
const int maxn=; LL rd(){
LL x=;char c=getchar();int neg=;
while(c<''||c>''){if(c=='-') neg=-;c=getchar();}
while(c>=''&&c<='') x=x*+c-'',c=getchar();
return x*neg;
} int N,stk[maxn],arr[maxn],tmp[maxn];
double f[maxn],A[maxn],B[maxn],R[maxn]; inline double gety(int i){return R[i]*f[i]/(A[i]*R[i]+B[i]);}
inline double getx(int i){return f[i]/(A[i]*R[i]+B[i]);}
inline double gett(int i){return -B[i]/A[i];}
inline double getk(int i,int j){
double dx=getx(i)-getx(j);if(fabs(dx)<=1e-) dx=1e-;
return (gety(i)-gety(j))/dx;
}
inline bool cmpl(int a,int b){return getx(a)<getx(b);}
inline bool cmpr(int a,int b){return gett(a)<gett(b);} void cdq(int l,int r){
if(l>=r){f[arr[l]]=max(f[arr[l]],f[arr[l]-]);return;}
int m=l+r>>,t=;
cdq(l,m);sort(arr+l,arr+m+,cmpl);
memcpy(tmp+m+,arr+m+,*(r-m));sort(arr+m+,arr+r+,cmpr);
for(int p=l;p<=m;p++){
while(t>=&&getk(stk[t],stk[t-])<getk(stk[t],arr[p])) t--;
stk[++t]=arr[p];
}for(int q=m+;q<=r;q++){
while(t>=&&getk(stk[t],stk[t-])<gett(arr[q])) t--;
int j=stk[t],i=arr[q];
f[i]=max(f[i],gety(j)*A[i]+getx(j)*B[i]);
}
memcpy(arr+m+,tmp+m+,*(r-m));cdq(m+,r);
} int main(){
//freopen("testdata.in","r",stdin);
//freopen(".out","w",stdout);
int i,j,k;
N=rd(),f[]=rd();
for(i=;i<=N;i++){
scanf("%lf%lf%lf",&A[i],&B[i],&R[i]);
arr[i]=i;
}cdq(,N);
double ans=;
for(i=;i<=N;i++){
ans=max(ans,f[i]);
}printf("%.3lf\n",ans);
return ;
}
bzoj1492/luogu4027 货币兑换 (斜率优化+cdq分治)的更多相关文章
- LUOGU P4027 [NOI2007]货币兑换 (斜率优化+CDQ分治)
传送门 解题思路 题目里有两句提示一定要看清楚,要不全买要不全卖,所以dp方程就比较好列,f[i]=max(f[j]*rate[j]*a[i])/(rate[j]*a[j]+b[j])+(f[j]*b ...
- BZOJ_3963_[WF2011]MachineWorks_斜率优化+CDQ分治
BZOJ_3963_[WF2011]MachineWorks_斜率优化+CDQ分治 Description 你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM) ...
- [BZOJ1492][NOI2007]货币兑换Cash(斜率优化+CDQ分治)
1492: [NOI2007]货币兑换Cash Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 5838 Solved: 2345[Submit][Sta ...
- 【BZOJ-1492】货币兑换Cash DP + 斜率优化 + CDQ分治
1492: [NOI2007]货币兑换Cash Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 3396 Solved: 1434[Submit][Sta ...
- 【BZOJ1492】[NOI2007]货币兑换Cash 斜率优化+cdq分治
[BZOJ10492][NOI2007]货币兑换Cash Description 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和 B纪念券(以下简称B券).每 ...
- BZOJ.1492.[NOI2007]货币兑换(DP 斜率优化 CDQ分治/Splay)
BZOJ 洛谷 如果某天能够赚钱,那么一定会在这天把手上的金券全卖掉.同样如果某天要买,一定会把所有钱花光. 那么令\(f_i\)表示到第\(i\)天所拥有的最多钱数(此时手上没有任何金券),可以选择 ...
- [Noi2014]购票 BZOJ3672 点分治+斜率优化+CDQ分治
Description 今年夏天,NOI在SZ市迎来了她30周岁的生日.来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会.全国的城市构成了一棵以SZ市为根的有根树,每个城市与它的 ...
- 洛谷.4655.[CEOI2017]Building Bridges(DP 斜率优化 CDQ分治)
LOJ 洛谷 \(f_i=s_{i-1}+h_i^2+\min\{f_j-s_j+h_j^2-2h_i2h_j\}\),显然可以斜率优化. \(f_i-s_{i-1}-h_i^2+2h_ih_j=f_ ...
- BZOJ3963 WF2011MachineWorks(动态规划+斜率优化+cdq分治)
按卖出时间排序后,设f[i]为买下第i台机器后的当前最大收益,则显然有f[i]=max{f[j]+gj*(di-dj-1)+rj-pi},且若此值<0,应设为-inf以表示无法购买第i台机器. ...
随机推荐
- nginx location url解析过程
- iOS开发简记(3):tips提示
我有一个需求:在点击或长按某个按钮时,需要显示提示,包括简单的文字提示,还有复杂一点的图片甚至是动态图的提示(可能还要加上文字). (1)文字tips 使用之前介绍的qmuikit里面的QMUITip ...
- Shell编程基础篇-上
1.1 前言 1.1.1 为什么学Shell Shell脚本语言是实现Linux/UNIX系统管理及自动化运维所必备的重要工具, Linux/UNIX系统的底层及基础应用软件的核心大都涉及Shell脚 ...
- windows平台下编辑的内容传到linux平台出现中文乱码的解决办法
现象说明:在windows下编辑的内容,上传到linux平台下出现中文乱码.如下: 在windows平台编写haha.txt文件,内容如下: 上传到linux平台,出现中文乱码,如下: 基本上面出现的 ...
- Jenkins部署Python项目实战
一.背景 我们工作中常用Jenkins部署Java代码,因其灵活的插件特性,例如jdk,maven,ant等使得java项目编译后上线部署一气呵成,同样对于脚本语言类型如Python上线部署,利用Je ...
- cf166e 在四面体上寻找路线数 递推,取模
来源:codeforces E. Tetrahedron You are given a tetrahedron. Let's mark its vertices ...
- 结对项目——Core设计与实现
写在前面:关于结对编程 结对编程我一直认为是一种非常好的合作方式,他的形式主要是由一个人负责代码编写,另一个人则在一旁即时对写下的代码进行审查,这样可以大大减少代码实现方面的错误. 这次我的结对伙伴是 ...
- 《Linux内核分析》第八周学习笔记
<Linux内核分析>第八周学习笔记 进程的切换和系统的一般执行过程 郭垚 原创作品转载请注明出处 <Linux内核分析>MOOC课程http://mooc.study.163 ...
- “数学口袋精灵”第二个Sprint计划(第六~八天)
“数学口袋精灵”第二个Sprint计划----第六天~第八天进度 任务分配: 冯美欣:欢迎界面的背景音乐完善 吴舒婷:游戏界面的动作条,选择答案后的音效 林欢雯:代码算法设计 第六天: 进度: 冯美欣 ...
- weex 开发踩坑日记--环境配置、安卓运行、adb、开发
环境配置方面 1.需要安装java和android环境,java的话一定要下载jdk而不是jre. 在"系统变量"新建一个变量名为JAVA_HOME的变量,变量值为你本地java的 ...