BZOJ 1492: [NOI2007]货币兑换Cash

强大的cdq+斜率优化 推荐一个PPT 传送门
斜率优化的式子在这里就不推了 上面的PPT写的也很好
按照原顺序 维护凸包 再用斜率去切凸包 大概就是这样？
可以用归并排序来对斜率排序 这样的时间复杂度是$O(nlogn)$
第一次写 特别蛋疼 所以几乎是直接搬了黄学长的代码。。
如果你不会写 可以先看代码 再打

#include<bits/stdc++.h>
using namespace std;
const int N=100002;
const double eps=1e-9;
struct node{double x,y,a,b,k,ra;int w,id;}p[N],t[N];
bool Cmp(node a,node b){return a.k>b.k;}
double f[N]; int n,tp,st[N];
double getk(int a,int b)
{
    if(fabs(p[a].x-p[b].x)<eps)return 1e20;
    return (p[b].y-p[a].y)/(p[b].x-p[a].x);
}
void solve(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].ra+p[l].b);
        p[l].x=p[l].y*p[l].ra;
        return;
    }
    int mid=(l+r)>>1,l1=l,l2=mid+1,j=1;
    for(int i=l;i<=r;i++)
        if(p[i].id<=mid)t[l1++]=p[i];
        else t[l2++]=p[i];
    for(int i=l;i<=r;i++)p[i]=t[i];
    solve(l,mid); tp=0;
    for(int i=l;i<=mid;i++)
    {
        while(tp>1 && getk(st[tp-1],st[tp])<getk(st[tp-1],i)+eps)tp--;
        st[++tp]=i;
    }
    for(int i=mid+1;i<=r;i++)
    {
        while(j<tp && getk(st[j],st[j+1])+eps>p[i].k)j++;
        f[p[i].id]=max(f[p[i].id],p[st[j]].x*p[i].a+p[st[j]].y*p[i].b);
    }
    solve(mid+1,r);
    l1=l,l2=mid+1;
    for(int i=l;i<=r;i++)
    {
        if( ( (p[l1].x<p[l2].x|| (fabs(p[l1].x-p[l2].x)<eps&&p[l1].y<p[l2].y) )|| l2>r) && l1<=mid)
            t[i]=p[l1++];
        else t[i]=p[l2++];
    }
    for(int i=l;i<=r;i++)p[i]=t[i];
}
int main()
{
    int i;
    scanf("%d%lf",&n,&f[0]);
    for(i=1;i<=n;i++)
    {
        scanf("%lf%lf%lf",&p[i].a,&p[i].b,&p[i].ra);
        p[i].k=-p[i].a/p[i].b,p[i].id=i;
    }
    sort(p+1,p+1+n,Cmp);
    solve(1,n);
    printf("%.3lf\n",f[n]);
    return 0;
}