bzoj 1492: [NOI2007]货币兑换Cash cdq分治

题意

题目

分析

有十分显然的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;
}