\(NOI2007\)货币兑换
辣鸡题目毁我青春
刚看到这题的时候我一脸懵逼啊，问了一波\(fdf\)才把六十分推出来，然后强行理解了一波斜率和横坐标都不单调的斜率优化。理解了一下发现这不是傻逼题吗，然后我调斜率计算调到怀疑人生。
\(Part1\) 关于斜率不存在
这个很显然啊，但是我之前居然都没被卡过。
根据两个点纵坐标来判断到底是返回\(inf\)还是\(-inf\)
\(Part2\) 关于两点重合
两点的纵坐标也相等时，是返回\(-inf\)的！别问我为什么浮点数比较用等于号是有用的，估计他们掉精度是有一个法则的吧，掉成一样的了。
\(Part3\) 关于弹栈
浮点数比较真的真的不要忽略等于号，在弹栈的时候，斜率相等的情况是要弹掉的，不然的话会\(WA\)，感觉可能是二分的时候会挂掉，这个等于号是很有用的。（我记得凸包也不能有三点共线吧）
挺裸的一道题，但是被上面的几个细节坑到了，积累一下，以后就不会出错了，挺好。
下面是\(CDQ\)写法的代码。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<set>
#include<queue>
#include<map>
#include<vector>
#define N 101000
#define EPS 1e-7
#define ll long long
#define RG register
#define inf 1e18+1
using namespace std;
inline int read(){
RG int x=0,o=1; RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
if(ch=='-') o=-1,ch=getchar();
while(ch>='0'&&ch<='9') x=((x<<3)+(x<<1))+ch-'0',ch=getchar();
return x*o;
}
int n,S,Q[N]; double f[N],a[N],b[N],r[N],ans;
struct mona { int id; double K,A,B; } s[N],p[N];
inline double Slope(const int &i,const int &j){
if(fabs(s[j].A-s[i].A)<EPS) return s[i].B<s[j].B?inf:-inf;
return (s[i].B-s[j].B)/(s[i].A-s[j].A);
}
inline void CDQ(int l,int r){
if(l==r){
f[l]=max(f[l],f[l-1]);
s[l].A*=f[l],s[l].B*=f[l]; return ;
} int mid=l+r>>1;
CDQ(l,mid); RG int top=1; Q[1]=0;
for(RG int i=l;i<=mid;++i){
while(top>1&&Slope(Q[top-1],Q[top])<=Slope(Q[top-1],i)) --top;
Q[++top]=i;
}
for(RG int i=mid+1;i<=r;++i){
int L=2,R=top,ans=1;
while(L<=R){
int Mid=L+R>>1;
if(Slope(Q[Mid-1],Q[Mid])>s[i].K) L=Mid+1,ans=Mid;
else R=Mid-1;
} f[i]=max(f[i],(-s[Q[ans]].A*s[i].K+s[Q[ans]].B)*b[i]);
} CDQ(mid+1,r);
for(RG int l1=l,l2=mid+1,len=l;len<=r;++len)
if((s[l1].A<s[l2].A&&l1<=mid)||l2>r) p[len]=s[l1++];
else p[len]=s[l2++];
for(RG int i=l;i<=r;++i) s[i]=p[i];
}
int main(){
n=read(),S=read(),f[1]=S;
for(RG int i=1;i<=n;++i) scanf("%lf%lf%lf",&a[i],&b[i],&r[i]);
for(RG int i=1;i<=n;++i)
s[i]=(mona) { i,-a[i]/b[i],r[i]/(a[i]*r[i]+b[i]),1.0/(a[i]*r[i]+b[i]) };
CDQ(1,n); printf("%.3lf\n",f[n]);
}