参考: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. Qt 安装与配置记录

    一 安装的时候得选一个Qt安装啊!!不要忘了展开这一项,而只安装Qt creator 展开之后会发现有很多版本,为了方便,选自带编译器mingw,就不需要麻烦的配置了 二 打开Qt creator 后 ...

  2. xtu summer individual 3 F - Opening Portals

    Opening Portals Time Limit: 2000ms Memory Limit: 262144KB This problem will be judged on CodeForces. ...

  3. IIS 和 ASP.NET ISAPI

    前几天有一个朋友在MSN上问我“ASP.NET 从最初的接收到Http request到最终生成Response的整个流程到底是怎样的?”我觉得这个问题涉及到IIS和ASP.NETASP.NET Ru ...

  4. linux shell & man chmod

    man chomd MBP xgqfrms-mbp:~ xgqfrms-mbp$ man chmod and entries will be removed regardless of their i ...

  5. HDU——2063 过山车

    过山车 Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submis ...

  6. 洛谷 P1023 税收与补贴问题

    P1023 税收与补贴问题 题目背景 每样商品的价格越低,其销量就会相应增大.现已知某种商品的成本及其在若干价位上的销量(产品不会低于成本销售),并假设相邻价位间销量的变化是线性的且在价格高于给定的最 ...

  7. TCP/IP协议族-----22、万维网和HTTP

  8. Django学习系列之Form表单结合ajax

      Forms结合ajax Forms的验证流程: 定义用户输入规则的类,字段的值必须等于html中name属性的值(pwd= forms.CharField(required=True)=<i ...

  9. System.AccessViolationException”类型的未经处理的异常在 System.Data.dll 中发生。其它信息:尝试读取或写入受保护的内存。这通常指示其它内存已损坏。

    错误背景: 操作系统:编程环境:VS2013.  语言:VB.net:  数据库:SQLserver2008 做数据库连接时.发生的错误: 错误提示为: 说明:用VB.net连接SQLServer数据 ...

  10. Excel数据字典转换为PDM(且显示表名、字段相应的中文描写叙述)

    在工作中遇到了一个问题就是把Excel数据字典转换为PDM. 可是转换完毕了全是英文,原来对表名.字段名的中文描写叙述就没有了. 且对于这个问题在powerdesigner15.2以后能够直接完毕.可 ...