[bzoj1492]货币兑换Cash

[bzoj1492]货币兑换Cash

---

不难推出式子 

dp[i]=max(A[j]*a[i]+B[j]*b[i])
//A[j]=f[j]/(a[j]+b[j]/rate[j]);
//B[j]=f[j]/(rate[j]*a[j]+b[j]);
//分别表示全部存入能得到的A,B金券数

cdq分治统计l~mid对mid+1~r的贡献

将左边的按A[i]的大小排序，维护下凸壳(A[i],B[i])，我们发现斜率的相对大小是不会随b[i]改变的，所以对于右边的所有点都可以在这个凸壳上跑(斜率乘上b[i])就行了

• 代码

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+5;
int n,s;
long double eps=1e-18;
long double rate[N],a[N],b[N];
long double f[N];
struct node{
    long double A,B;
    bool operator < (const node b)const{
        if(fabs(A-b.A)<=eps)return B>b.B;
        return A<b.A;
    }
}nd[N];

int id[N];
int st[N*2],top;

long double calcx(int l,int r){
    return nd[r].A-nd[l].A;
}
long double calcy(int l,int r){
    return nd[r].B-nd[l].B;
}

inline void merge(int l,int r){
    if(l==r)return;
    int mid=(l+r)>>1;
    merge(l,mid);
    long double mx=0;
    for(int i=l;i<=mid;i++){
        nd[i].A=f[i]/(a[i]+b[i]/rate[i]);
        nd[i].B=f[i]/(rate[i]*a[i]+b[i]);
        mx=max(mx,f[i]);
    }
    sort(nd+l,nd+mid+1);
    top=0;
    for(int i=l;i<=mid;i++){
        while(top>1&&
        calcy(st[top-1],i)*calcx(st[top-1],st[top])+eps > calcy(st[top-1],st[top])*calcx(st[top-1],i) )top--;
        st[++top]=i;
    }
    for(int i=mid+1;i<=r;i++)f[i]=max(f[i],mx);
    for(int i=mid+1;i<=r;i++){
        long double k=a[i];
        int l=1,r=top;
        while(l<r){
            int mid=(l+r)>>1;
            f[i]=max(f[i],nd[st[mid]].A*a[i]+nd[st[mid]].B*b[i]);
            if( calcy(st[mid],st[mid+1])*b[i] < eps-k*calcx(st[mid],st[mid+1]) )
                r=mid;
            else
                l=mid+1;
        }
        f[i]=max(f[i],nd[st[l]].A*a[i]+nd[st[l]].B*b[i]);
    }
    merge(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],&rate[i]);
    f[1]=f[0]=(long double)s;
    merge(1,n);
    printf("%.3Lf\n",f[n]);
}