参考:http://www.cnblogs.com/lidaxin/p/5240220.html

虽然splay会方便很多,但是懒得写,于是写了cdq

首先要想到贪心的思路,因为如果在某天买入是能得到最大收益的,那么应该用所有钱去买,相对的如果在某天卖出能得到最大收益,那么应该全部卖出

方便起见,设\( x[j]=f[j]/(a[j]*rate[j]+b[j])*rate[j] \)表示第j天最多可以拥有的A货币的数量,y[j]=f[j]/(a[j]*rate[j]+b[j])表示第j天最多可以拥有的B货币的数量

于是dp方程就可以写作\( f[i]=max{f[i-1],x[j]*a[i]+y[j]*b[i]} \)

然后按着斜率优化套路推一推就变成了\( y[j]=-a[i]/b[i]*x[j]+f[i]/b[i] \),按照y=kx+b的形式就是y=y[j],k=-a[i]/b[i],x=x[j],b=f[i]/b[i]

然后问题在于这个斜率不是单调的,所以考虑如何将点插入上凸壳

splay显然是可以的,但是写起来太麻烦了

仔细想想,dp的过程就相当于n次询问+插入

这个似乎可以用cdq分治,设f[i]的id为i

首先按照斜率k排序,对于当前要处理的区间(l,r),先把id<=mid的丢到(l,mid)中,把剩下的丢进(mid+1,r),递归处理出(l,mid),然后用处理出的点来构成一个上凸壳,用这个上凸壳来更新(mid+1,r)中的f,因为是按斜率排序的,所以用双指针即可

然后递归处理(mid+1,r),最后把点集按极角排序即可

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=100005,inf=1e9;
int n,s[N];
double f[N];
struct dian
{
double x,y;
bool operator < (const dian &b) const
{
return (x<=b.x)||(x==b.x&&y<=b.y);
}
}p[N],q[N];
struct qwe
{
double q,a,b,r,k;
int id;
}a[N],b[N];
double slv(int i,int j)
{
if(i==0)
return -inf;
if(j==0)
return inf;
if(p[i].x==p[j].x)
return -inf;
return (p[i].y-p[j].y)/(p[i].x-p[j].x);
}
bool cmp(const qwe &a,const qwe &b)
{
return a.k<b.k;
}
void wk(int l,int r)
{
if(l==r)
{
f[l]=max(f[l],f[l-1]);
p[l].y=f[l]/(a[l].a*a[l].r+a[l].b);
p[l].x=p[l].y*a[l].r;
return;
}
int mid=(l+r)>>1,top=0,l1=l,l2=mid+1;
for(int i=l;i<=r;i++)
{
if(a[i].id<=mid)
b[l1++]=a[i];
else
b[l2++]=a[i];
}
for(int i=l;i<=r;i++)
a[i]=b[i];
wk(l,mid);
for(int i=l;i<=mid;i++)
{
while(top>=2&&slv(i,s[top])>=slv(s[top],s[top-1]))
top--;
s[++top]=i;
}
for(int i=r,j=1;i>=mid+1;i--)
{
while(j<top&&a[i].k<=slv(s[j],s[j+1]))
j++;
f[a[i].id]=max(f[a[i].id],p[s[j]].x*a[i].a+p[s[j]].y*a[i].b);
}
wk(mid+1,r);
l1=l,l2=mid+1;
for(int i=l;i<=r;i++)
{
if((p[l1]<p[l2]||l2>r)&&l1<=mid)
q[i]=p[l1++];
else
q[i]=p[l2++];
}
for(int i=l;i<=r;i++)
p[i]=q[i];
}
int main()
{
scanf("%d%lf",&n,&f[0]);
for(int i=1;i<=n;i++)
{
scanf("%lf%lf%lf",&a[i].a,&a[i].b,&a[i].r);
a[i].k=-a[i].a/a[i].b,a[i].id=i;
}
sort(a+1,a+1+n,cmp);
wk(1,n);
printf("%.3f\n",f[n]);
return 0;
}

bzoj 1492: [NOI2007]货币兑换Cash【贪心+斜率优化dp+cdq】的更多相关文章

  1. BZOJ 1492 [NOI2007]货币兑换Cash:斜率优化dp + cdq分治

    传送门 题意 初始时你有 $ s $ 元,接下来有 $ n $ 天. 在第 $ i $ 天,A券的价值为 $ A[i] $ ,B券的价值为 $ B[i] $ . 在第 $ i $ 天,你可以进行两种操 ...

  2. bzoj 1492 [NOI2007]货币兑换Cash(斜率dp+cdq分治)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=1492   [题意] 有AB两种货币,每天可以可以付IPi元,买到A券和B券,且A:B= ...

  3. 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的数量. 然后就很明显了...平衡 ...

  4. ●BZOJ 1492 [NOI2007]货币兑换Cash

    题链: http://www.lydsy.com/JudgeOnline/problem.php?id=1492 题解: 斜率优化DP,CDQ分治 定义$DP[i]$为第i天结束后的最大收益. 由于题 ...

  5. BZOJ 1492: [NOI2007]货币兑换Cash [CDQ分治 斜率优化DP]

    传送门 题意:不想写... 扔链接就跑 好吧我回来了 首先发现每次兑换一定是全部兑换,因为你兑换说明有利可图,是为了后面的某一天两种卷的汇率差别明显而兑换 那么一定拿全利啊,一定比多天的组合好 $f[ ...

  6. bzoj 1492: [NOI2007]货币兑换Cash

    Description 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和 B纪念券(以下 简称B券).每个持有金券的顾客都有一个自己的帐户.金券的数目可以是一个 ...

  7. BZOJ 1492: [NOI2007]货币兑换Cash 斜率优化 + splay动态维护凸包

    Description 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和 B纪念券(以下 简称B券).每个持有金券的顾客都有一个自己的帐户.金券的数目可以是一个 ...

  8. 斜率优化(CDQ分治,Splay平衡树):BZOJ 1492: [NOI2007]货币兑换Cash

    Description Input 第一行两个正整数N.S,分别表示小Y 能预知的天数以及初始时拥有的钱数. 接下来N 行,第K 行三个实数AK.BK.RateK,意义如题目中所述 Output 只有 ...

  9. BZOJ 1492 [NOI2007]货币兑换Cash (CDQ分治/splay 维护凸包)

    题目大意:太长了略 splay调了两天一直WA弃疗了 首先,我们可以猜一个贪心,如果买/卖,就一定都买/卖掉,否则不买/卖 反正货币的行情都是已知的,没有任何风险,所以肯定要选择最最最优的方案了 容易 ...

随机推荐

  1. Linux学习笔记01

    1.Linux不靠扩展名区分文件类型2.存储设备必须先挂载才能使用3.Windows下的程序不能直接在Linux中安装和运行 一.服务器的管理预配置Linux的目录的作用:/bin/存放系统命令的目录 ...

  2. Linux下汇编语言学习笔记17 ---

    这是17年暑假学习Linux汇编语言的笔记记录,参考书目为清华大学出版社 Jeff Duntemann著 梁晓辉译<汇编语言基于Linux环境>的书,喜欢看原版书的同学可以看<Ass ...

  3. QT-Embedded-4.5.3在海思35xx上移植

    QT4.5.3在海思3520A上移植步骤-修订版 2015年3月29日星期日, 16:59:03 1.首先要保证已经安装了海思的交叉编译器: #arm-hi  + Tab key to show wh ...

  4. 【进击后端】mongodb入门

    环境:ubuntu 重要字:数据库,集合,记录,字段 api地址:https://docs.mongodb.com/manual/reference/method/db.collection.upda ...

  5. [bzoj3513][MUTC2013]idiots_FFT

    idiots bzoj-3513 MUTC-2013 题目大意:给定$n$根木棍,问随机选择三根能构成三角形的概率. 注释:$1\le n\le 3\cdot 10^5$,$1\le a_i\le 1 ...

  6. mysql常用jar包

    连接Mysql数据库: 常用的连接池有两种 DBCP连接池 C3P0连接池 Apache的commons组件 -- DBCP连接池: commons-dbutils-1.4.jar 封装并简化了JDB ...

  7. mysql你确定掌握的那些sql语句

    1.创建表 create table test(uid int not null,create_time timestamp default current_timestamp); 即:没有双引号,单 ...

  8. HDU 5319 Painter(枚举)

    Painter Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Su ...

  9. android自己定义控件系列教程----视图

    理解android视图 对于android设备我们所示区域事实上和它在底层的绘制有着非常大的关系,非常多时候我们都仅仅关心我们所示,那么在底层一点它究竟是怎么样的一个东西呢?让我们先来看看这个图. w ...

  10. android意图传參数(四)

    一.依照向导创建一个project,layout的activity_main.xml文件内容例如以下: <RelativeLayout xmlns:android="http://sc ...