[NOI2007]货币兑换

货币兑换

题解

忽然发现以前写过一篇的，但这次的做法与以前有所不同，常数也优秀了许多。

首先明确一个结论，最优的做法股票一定是一次买完一次卖完的，毕竟它又不收税。
如果我只卖一部分或只买一部分就能够赚钱的话，那我全部买完卖完就一定能赚更多。

有了上面的结论，我们可以考虑dp。
设在第$i$天最多能得到$f_i$的钱，答案就是$f_n$。
第$i$天能买A卷，$x_i=\frac{f_i{r}_i}{a_i{r}_i+b_i}$。
第$i$天能买B卷，$y_i=\frac{f_i}{a_i{r}_i+b_i}$。
可得状态转移方程
$$f_i=\underset{j\in [1,j)}{\max }{a}_i{x}_j+b_i{y}_j$$
$$\frac{f_i}{b_i}=\underset{j\in [1,j)}{\max }{x}_j\frac{a_i}{b_i}+y_j$$
这样就很明显是一个二次函数的形式了，由于$a_i,b_i$都是定值，所以我们要求的是在$\frac{a_i}{b_i}$处最大的一条线段$x_{j}x+y_j$，其中$x_j$是斜率，$y_j$是截距。看着总是怪怪的

很明显，我们可以维护一个凸包，表示在第$i$位最优的一条线段。
考虑李超线段树，我们很容易就可以将这些线段全部挂到树上，找$\frac{a_i}{b_i}$上的最优值。
可以先将所有的$\frac{a_i}{b_i}$求出来，我们要求的只有这些位置上的最优线段。
将它们排序后，作李超线段树的下标。

时间复杂度$O\left(nlogn\right)$。

源码

注意输入的$a,b,r$不是整数。

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<map>
using namespace std;
#define MAXN 100005
#define lowbit(x) (x&-x)
#define reg register
typedef long long LL;
typedef unsigned long long uLL;
typedef pair<int,int> pii;
const LL INF=0x7f7f7f7f7f7f;
const int mo=1e9+7;
const double PI=acos(-1.0);
template<typename _T>
_T Fabs(_T x){return x<0?-x:x;}
template<typename _T>
void read(_T &x){
	_T f=1;x=0;char s=getchar();
	while(s>'9'||s<'0'){if(s=='-')f=-1;s=getchar();}
	while('0'<=s&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=getchar();}
	x*=f;
}
int n,s,tot;double f[MAXN],b[MAXN];
struct ming{double a,b,r,x,y,k;}a[MAXN];
struct line{
	double k,b;
	double ask(const double x){return k*x+b;}
}tr[MAXN<<2];
void insert(int rt,int l,int r,line aw){
	if(l>r)return ;const int mid=l+r>>1;
	if(tr[rt].ask(b[mid])<aw.ask(b[mid]))swap(tr[rt],aw);
	if(tr[rt].ask(b[l])<aw.ask(b[l]))insert(rt<<1,l,mid,aw);
	if(tr[rt].ask(b[r])<aw.ask(b[r]))insert(rt<<1|1,mid,r,aw);
}
double query(int rt,int l,int r,const double ai){
	if(l>r)return 0;const int mid=l+r>>1;double res=0;
	res=tr[rt].ask(ai);if(l==r)return res;
	if(ai<b[mid])res=max(res,query(rt<<1,l,mid,ai));
	if(ai>b[mid])res=max(res,query(rt<<1|1,mid+1,r,ai));
	return res;
}
signed main(){
	read(n);read(s);f[1]=s;
	for(int i=1;i<=n;i++)
		scanf("%lf %lf %lf",&a[i].a,&a[i].b,&a[i].r),
		b[++tot]=a[i].k=a[i].a/a[i].b;
	sort(b+1,b+tot+1);
	a[1].x=f[1]*a[1].r/(a[1].a*a[1].r+a[1].b);
	a[1].y=f[1]/(a[1].a*a[1].r+a[1].b);
	insert(1,1,n,(line){a[1].x,a[1].y});
	for(int i=2;i<=n;i++){
		f[i]=max(f[i-1],a[i].b*query(1,1,n,a[i].k));
		a[i].x=f[i]*a[i].r/(a[i].a*a[i].r+a[i].b);
		a[i].y=f[i]/(a[i].a*a[i].r+a[i].b);
		insert(1,1,n,(line){a[i].x,a[i].y});
	}
	printf("%.3f",f[n]);
	return 0;
}

谢谢！！！