题意
分析
有十分显然的dp做法,设f[i]表示到第i天能得到的最大现金。
不难得到
f[i]=max(f[j]∗rate[j]/(b[j]+a[j]∗rate[j])∗a[i]+f[j]/(b[j]+a[j]∗rate[j])∗b[i])
设
x[j]=f[j]∗rate[j]/(b[j]+a[j]∗rate[j]),y[j]=f[j]/(b[j]+a[j]∗rate[j])
那么有
f[i]=max(x[j]∗a[i]+y[j]∗b[i])
移项后可得
f[j]=f[i]/b[i]−x[j]∗a[i]/b[i]
问题就转化成了给出若干个点
(x[j],y[j])
,求一条斜率为
−a[i]/b[i]
且过其中一个点的直线,使得其截距最大。
显然答案只可能在上凸壳上。有一种方法是用splay来维护上凸壳,还可以用cdq分治来做。
这里推荐比较好写的cdq分治做法。
代码
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=100005;
const double eps=1e-10;
int n,q[N],s;
double a[N],b[N],k[N],f[N];
struct data{double x,y;int id;}p[N],tmp[N],line[N];
bool cmp(data a,data b)
{
return a.x<b.x||fabs(a.x-b.x)<=eps&&a.y>b.y;
}
double get_k(data p,data q)
{
if (p.x==q.x) return -1e15;
return (p.y-q.y)/(p.x-q.x);
}
void solve(int l,int r)
{
if (l==r) return;
int mid=(l+r)/2;
solve(l,mid);
double mx=0;
for (int i=l;i<=mid;i++) mx=max(mx,f[i]);
for (int j=mid+1;j<=r;j++) f[j]=max(f[j],mx);
int tot=0;
for (int i=l;i<=mid;i++) tmp[++tot].x=f[i]*p[i].x,tmp[tot].y=f[i]*p[i].y;
sort(tmp+1,tmp+tot+1,cmp);
int head=1,tail=1;q[1]=1;
for (int i=2;i<=tot;i++)
{
while (head<tail&&get_k(tmp[q[tail-1]],tmp[q[tail]])<get_k(tmp[q[tail]],tmp[i])) tail--;
q[++tail]=i;
}
int cnt=0;
for (int i=mid+1;i<=r;i++) line[++cnt].x=a[i]/b[i],line[cnt].id=i;
sort(line+1,line+cnt+1,cmp);
for (int i=1;i<=cnt;i++)
{
int pos=line[i].id;
while (head<tail)
{
double w1=tmp[q[head]].x*a[pos]+tmp[q[head]].y*b[pos],w2=tmp[q[head+1]].x*a[pos]+tmp[q[head+1]].y*b[pos];
if (w1<w2||fabs(w1-w2)<=eps) head++;
else break;
}
f[pos]=max(f[pos],tmp[q[head]].x*a[pos]+tmp[q[head]].y*b[pos]);
}
solve(mid+1,r);
}
int main()
{
scanf("%d%d",&n,&s);
for (int i=1;i<=n;i++)
{
scanf("%lf%lf%lf",&a[i],&b[i],&k[i]);
p[i].x=(double)k[i]/(b[i]+a[i]*k[i]);
p[i].y=(double)1/(b[i]+a[i]*k[i]);
f[i]=s;
}
solve(1,n);
double ans=0;
for (int i=1;i<=n;i++) ans=max(ans,f[i]);
printf("%.3lf",ans);
return 0;
}