BZOJ_1492_[NOI2007]货币兑换Cash_CDQ分治+斜率优化
Description
![](https://www.lydsy.com/JudgeOnline/upload/201604/dd%281%29.png)
![](https://www.lydsy.com/JudgeOnline/upload/201604/dd%282%29.png)
Input
Output
只有一个实数MaxProfit,表示第N天的操作结束时能够获得的最大的金钱数目。答案保留3位小数。
Sample Input
1 1 1
1 2 2
2 2 3
Sample Output
HINT
首先贪心的想,每天的操作只有三种可能,什么也不干,全部买入和全部卖出。
设x[i]为第i天能获得A卷的数量,y[i]为第i天能获得B卷的数量。
y[i]=f[i]/(rate[i]*A[i]+B[i]),x[i]=y[i]*rate[i]。
那么f[i]=max(f[i-1],x[j]*A[i]+y[j]*B[i]),其中我们需要找到一个已经确定F值得j来转移这个过程。
F[i]=X[j]*A[i]+Y[j]*B[i]
F[i]/A[i]=Y[j]*(B[i]/A[i])+X[j]。
把Y[j]当做斜率。那么我们会发现不仅Y[j]不单调,每次切的横坐标(B[i]/A[i])也不单调。
这时候我们使用CDQ分治。
本题的思想:对左边进行处理,处理左边对右边的影响(F的转移),再对右边进行处理。
先按(B[i]/A[i])排序,然后每步再按斜率归并上去。
也就是说我们到了一个局面,满足左边的斜率单调,右边每次切的横坐标单调。
这个就很好办了,对左边单调栈求凸包,然后维护一个指针扫一遍凸包。
正确性:所有能更新F[i]的F[j]在更新F[i]之前都已经更新完毕了,同时F[i]被1~i-1中的每个F[j]都更新到了。
代码:
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;
typedef double f2;
#define N 100050
#define eps 1e-6
f2 xx[N],yy[N],f[N],A[N],B[N],rate[N],K[N],pos[N];
int n,t[N],tmp[N],S[N];
#define Y(i,j) (K[j]*pos[i]+xx[j])
bool cmp(int x,int y) {return pos[x]<pos[y];}
bool judge(int p1,int p2,int p3) {
return (yy[p2]-yy[p3])*(xx[p1]-xx[p2])<=(yy[p1]-yy[p2])*(xx[p2]-xx[p3]);
}
void solve(int l,int r) {
if(l==r) {
f[l]=max(f[l],f[l-1]);
yy[l]=f[l]/(rate[l]*A[l]+B[l]);
xx[l]=yy[l]*rate[l];
K[l]=yy[l];
return ;
}
int mid=(l+r)>>1;
int i,j=l,k=mid+1;
for(i=l;i<=r;i++) {
if(t[i]<=mid) tmp[j++]=t[i];
else tmp[k++]=t[i];
}
for(i=l;i<=r;i++) t[i]=tmp[i];
solve(l,mid);
int top=0;
for(i=l;i<=mid;i++) {
while(top>1&&judge(S[top-1],S[top],t[i])) top--;
S[++top]=t[i];
}
for(j=1,i=mid+1;i<=r;i++) {
while(j<top&&Y(t[i],S[j+1])>=Y(t[i],S[j])) j++;
f[t[i]]=max(f[t[i]],Y(t[i],S[j])*A[t[i]]);
}
solve(mid+1,r);
i=j=l,k=mid+1;
while(j<=mid&&k<=r) {
if(K[t[j]]<=K[t[k]]) tmp[i++]=t[j++];
else tmp[i++]=t[k++];
}
while(j<=mid) tmp[i++]=t[j++];
while(k<=r) tmp[i++]=t[k++];
for(i=l;i<=r;i++) t[i]=tmp[i];
}
int main() {
scanf("%d%lf",&n,&f[0]);
int i,j;
for(i=1;i<=n;i++) {
scanf("%lf%lf%lf",&A[i],&B[i],&rate[i]);
t[i]=i; pos[i]=B[i]/A[i];
}
sort(t+1,t+n+1,cmp);
solve(1,n);
printf("%.3f\n",f[n]);
}