BZOJ 1492 [NOI2007]货币兑换Cash (CDQ分治/splay 维护凸包)
题目大意:太长了略
splay调了两天一直WA弃疗了
首先,我们可以猜一个贪心,如果买/卖,就一定都买/卖掉,否则不买/卖
反正货币的行情都是已知的,没有任何风险,所以肯定要选择最最最优的方案了
容易得到方程
$dp[i]=max(dp[i-1],a[i]*\frac{dp[j]*rate[j]}{rate[j]*a[j]+b[j]}+b[i]*\frac{dp[j]}{rate[j]*a[j]+b[j]})$
显然是要用凸优化了
splay非常无脑,splay维护此题的凸包,需要找前驱,删前驱,找后继,删后继,一大堆特判...绝对恶心到吐
所以这是一篇$CDQ$分治题解
令$x=\frac{dp[j]}{rate[j]*a[j]+b[j]},y=x*rate[j]$
移项,可得
$dp[i]-b[i]*x=a[i]*y$
$y=\frac{dp[i]}{a[i]}-\frac{b[i]}{a[i]}x$
发现斜率$k=-\frac{b[i]}{a[i]}$是一定的,我们在外层把斜率k从小到大排序,可以优化掉一个$log$,递归时按$x$从小到大排序
这样,递归时,每一层内部都是按$k$有序的,把这一层按照时间分为左右两个部分(不要破坏$k$的有序状态)
先递归处理左半个区间,回溯后,左半部分的答案已知,且不会被右半部分的答案所影响
且左半部分按$x$从小到大排序,右半部分按斜率$k$从小到大排序,取最小值,由于$k<0$,用单调栈维护一个上凸包即可
处理完了左边对右边的贡献,递归处理右半部分
回溯时,先处理$dp[i]=dp[i-1]$的情况,再按$x$排序,回溯到上一层
#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N1 101000
#define M1 205
#define ll long long
#define dd double
#define uint unsigned int
#define inf 233333333
#define il inline
using namespace std; const dd eps=(1e-);
int n,m;
int stk[N1];
dd A[N1],B[N1],r[N1],X[N1],Y[N1],K[N1],f[N1];
//struct node{dd x,y,k,ans;int id;};
int cmp1(int s1,int s2){return K[s1]-K[s2]<;}
int id[N1],tmp[N1];
dd get_slope(int s1,int s2){
if(!s2) return inf;
return (Y[s1]-Y[s2])/(X[s1]-X[s2]);
}
void CDQ(int L,int R)
{
if(R-L<=) return;
int M=(L+R)>>;
int tp=,i,j,pl=L,pr=M,k,cnt;
for(int i=L;i<R;i++)
if(id[i]<M) tmp[pl++]=id[i];
else tmp[pr++]=id[i];
for(int i=L;i<R;i++)
id[i]=tmp[i];
CDQ(L,M);
for(i=L;i<M;i++)
{
k=id[i];
if(tp>&&fabs(X[stk[tp]]-X[k])<eps&&Y[k]-Y[stk[tp]]<eps) continue;
while(tp>&&get_slope(stk[tp],stk[tp-])<=get_slope(k,stk[tp-]))
tp--;
stk[++tp]=k;
}
for(i=M;i<R;i++)
f[i]=max(f[i-],f[i]);
for(i=M;i<R;i++)
{
while(tp>&&get_slope(stk[tp],stk[tp-])<=K[id[i]])
tp--;
k=id[i],j=stk[tp];
f[k]=max(f[k],A[k]*Y[j]+B[k]*X[j]);
X[k]=f[k]/(A[k]*r[k]+B[k]);
Y[k]=r[k]*X[k];
}
CDQ(M,R);
i=L,cnt=L,j=M;
while(i<M&&j<R){
if(X[id[i]]<X[id[j]])
tmp[cnt++]=id[i],i++;
else
tmp[cnt++]=id[j],j++;
}
while(i<M) tmp[cnt++]=id[i],i++;
while(j<R) tmp[cnt++]=id[j],j++;
for(i=L;i<R;i++)
id[i]=tmp[i],f[id[i]]=max(f[id[i]],f[id[i]-]);
};
dd S;
int main()
{
//freopen("t1.in","r",stdin);
scanf("%d%lf",&n,&S);
for(int i=;i<=n;i++)
{
scanf("%lf%lf%lf",&A[i],&B[i],&r[i]);
K[i]=-B[i]/A[i];id[i]=i;
}
f[]=S,X[]=S/(A[]*r[]+B[]),Y[]=X[]*r[];
sort(id+,id+n+,cmp1);
CDQ(,n+);
dd ans=;
for(int i=;i<=n;i++)
ans=max(ans,f[i]);
printf("%.3lf\n",ans);
return ;
}
BZOJ 1492 [NOI2007]货币兑换Cash (CDQ分治/splay 维护凸包)的更多相关文章
- BZOJ 1492: [NOI2007]货币兑换Cash [CDQ分治 斜率优化DP]
传送门 题意:不想写... 扔链接就跑 好吧我回来了 首先发现每次兑换一定是全部兑换,因为你兑换说明有利可图,是为了后面的某一天两种卷的汇率差别明显而兑换 那么一定拿全利啊,一定比多天的组合好 $f[ ...
- BZOJ 1492: [NOI2007]货币兑换Cash 斜率优化 + splay动态维护凸包
Description 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和 B纪念券(以下 简称B券).每个持有金券的顾客都有一个自己的帐户.金券的数目可以是一个 ...
- [BZOJ1492] [NOI2007] 货币兑换Cash(cdq分治+斜率优化)
[BZOJ1492] [NOI2007] 货币兑换Cash(cdq分治+斜率优化) 题面 分析 dp方程推导 显然,必然存在一种最优的买卖方案满足:每次买进操作使用完所有的人民币:每次卖出操作卖出所有 ...
- bzoj1492[NOI2007]货币兑换Cash cdq分治+斜率优化dp
1492: [NOI2007]货币兑换Cash Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 5541 Solved: 2228[Submit][Sta ...
- BZOJ 1492: [NOI2007]货币兑换Cash( dp + 平衡树 )
dp(i) = max(dp(i-1), x[j]*a[i]+y[j]*b[i]), 0<j<i. x, y表示某天拥有的最多钱去买金券, 金券a和金券b的数量. 然后就很明显了...平衡 ...
- bzoj 1492 [NOI2007]货币兑换Cash(斜率dp+cdq分治)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=1492 [题意] 有AB两种货币,每天可以可以付IPi元,买到A券和B券,且A:B= ...
- bzoj 1492: [NOI2007]货币兑换Cash【贪心+斜率优化dp+cdq】
参考:http://www.cnblogs.com/lidaxin/p/5240220.html 虽然splay会方便很多,但是懒得写,于是写了cdq 首先要想到贪心的思路,因为如果在某天买入是能得到 ...
- ●BZOJ 1492 [NOI2007]货币兑换Cash
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=1492 题解: 斜率优化DP,CDQ分治 定义$DP[i]$为第i天结束后的最大收益. 由于题 ...
- bzoj 1492: [NOI2007]货币兑换Cash
Description 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和 B纪念券(以下 简称B券).每个持有金券的顾客都有一个自己的帐户.金券的数目可以是一个 ...
随机推荐
- IOS - 零碎
---恢复内容开始--- 1.模拟器目录: ProjectNameApk.documents.library(cache.preference.cookies).temp 2.Edit-Refacto ...
- nyoj56-阶乘因式分解(一)
56-阶乘因式分解(一) 内存限制:64MB时间限制:3000msSpecial Judge: No accepted:6submit:7 题目描述: 给定两个数m,n,其中m是一个素数. 将n(0& ...
- Python for json
# Json- json格式是一个键值对形式的数据集 - key:字符串 - value:字符串,数字,列表,json - json使用大括号包裹 - 键值对直接用逗号隔开 - json和Python ...
- $_SERVER 详解
$_SERVER['HTTP_ACCEPT_LANGUAGE']//浏览器语言 $_SERVER['REMOTE_ADDR'] //当前用户 IP . $_SERVER['REMOTE_HOST'] ...
- codevs——T3657 括号序列
http://codevs.cn/problem/3657/ 时间限制: 1 s 空间限制: 256000 KB 题目等级 : 黄金 Gold 题解 题目描述 Description ...
- POJ 1320
作弊了--!该题可以通过因式分解得到一个佩尔方程....要不是学着这章,估计想不到.. 得到x1,y1后,就直接代入递推式递推了 x[n]=x[n-1]*x[1]+d*y[n-1]*y[1] y[n] ...
- 剑指Offer读书笔记(持续更新中)
(1)定义一个空的类型,里面没有不论什么成员变量和成员函数,对该类型求sizeof,得到的结果是多少? 答案是1.空类型的实例中不包括不论什么信息,本来求sizeof应该是0,可是当我们声明该类型实例 ...
- TT流程随笔
细节: 如果本地可以自动登录, 先实现本地登录,发送事件通知,再请求登录服务器 如果本地不可以登录(第一次或退出后),直接请求登录服务器 登录服务器返回消息服务器ip port / 文件服务器 链接消 ...
- 《深入理解Android 卷III》第六章 深入理解控件(ViewRoot)系统
<深入理解Android 卷III>即将公布,作者是张大伟.此书填补了深入理解Android Framework卷中的一个主要空白,即Android Framework中和UI相关的部分. ...
- 由GPS坐标计算半径
在实际应用当中,一般是通过一个个体的编码来查找该编码对应的地区中心的经纬度,然后再根据这些经纬度来计算彼此的距离,从而估算出某些群体之间的大致距离范围(比如酒店旅客的分布范围-各个旅客的邮政编码对应的 ...