【BZOJ 1492】 [NOI2007]货币兑换Cash

1492: [NOI2007]货币兑换Cash

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 1948  Solved: 879
[Submit][Status][Discuss]

Description

Input

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

Output

只有一个实数MaxProfit，表示第N 天的操作结束时能够获得的最大的金钱 数目。答案保留3 位小数。

Sample Input

3 100
1 1 1
1 2 2
2 2 3

Sample Output

225.000

HINT

测试数据设计使得精度误差不会超过10-7。
对于40%的测试数据，满足N ≤ 10；
对于60%的测试数据，满足N ≤ 1 000；
对于100%的测试数据，满足N ≤ 100 000；

dp+斜率优化+CDQ分治。

f[i]表示第i天后的最大获利。

f[i]=max(f[i-1],xj*ai+yj*bi)

xj,yj表示把第j天的最大获利f[j]都买成纪念券，能买到的A券和B券的数量。

设p=xj*ai+yj*bi

那么yj=-ai/bi*xj+p/bi

把（xj,yj）看成点，p是截距，就是要让一条斜率为-ai/bi的直线从无穷高的地方向下平移，第一个碰到的点(xj,yj)就是使得答案最优的点。

可以发现最优的点在之前的(xj,yj)点形成的上凸壳上，且斜率为-ai/bi的后继的边的靠左的点最优。

显然可以用splay来维护，但是CDQ分治要好写的多。

把斜率降序排序，在递归解决完(l,mid)之后，计算(l,mid)对(mid+1,r)的影响：

将(l,mid)的点建成上凸壳，因为(mid+1,r)是按照斜率降序排序的，所以只要扫一次凸壳就可以把(mid+1,r)全部更新。

注意在递归计算完(mid+1,r)之后，要按照x升序排序，因为当前的(l,r)全部算完了，他的用处就是建凸壳了。

#include <iostream>#include <algorithm>#include <cstdio>#include <cstring>#include <cstdlib>#include <cmath>#define eps 1e-9#define M 100005using namespace std;double f[M];int n,s[M];struct data{double x,y,a,b,k,rate;int id;}p[M],q[M];bool cmp(data a,data b){return a.k>b.k;}double Getk(int x,int y){if (!y) return -1e20;if (fabs(p[x].x-p[y].x)<eps) return 1e20;return (p[x].y-p[y].y)/(p[x].x-p[y].x);}void CDQ(int l,int r){if (l==r){f[l]=max(f[l],f[l-1]);p[l].y=f[l]/(p[l].a*p[l].rate+p[l].b);p[l].x=p[l].y*p[l].rate;return;}int l1=l,mid=(l+r)>>1,l2=mid+1;for (int i=l;i<=r;i++)if (p[i].id<=mid) q[l1++]=p[i];    else q[l2++]=p[i];    for (int i=l;i<=r;i++)p[i]=q[i];CDQ(l,mid);int top=0;for (int i=l;i<=mid;i++){while (top>1&&Getk(s[top-1],s[top])<eps+Getk(s[top],i))top--;s[++top]=i;}s[++top]=0;int j=1;for (int i=mid+1;i<=r;i++){while (j<top&&Getk(s[j],s[j+1])+eps>p[i].k)j++;f[p[i].id]=max(f[p[i].id],p[s[j]].x*p[i].a+p[s[j]].y*p[i].b);}CDQ(mid+1,r);l1=l,l2=mid+1;for (int i=l;i<=r;i++)if ((p[l1].x<p[l2].x||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",&p[i].a,&p[i].b,&p[i].rate);p[i].k=-p[i].a/p[i].b;p[i].id=i;}sort(p+1,p+n+1,cmp);CDQ(1,n);printf("%.3lf\n",f[n]);return 0;}