[NOI2007]货币兑换

题目

传送门 to luogu

思路

都是利用「买则疯狂买，卖则疯狂卖」进行 $\mathtt{d}\mathtt{p}$ 的。而这种性质也很容易证明：如果第 $i$ 天买入、第 $j$ 天卖出是有利可图的，那么增加本金只会赚更多。

回到题目上。用 $f(i)$ 表示第 $i$ 天卖光金券后可以获得最多的钱。考虑一下第 $j$ 天买入、第 $i$ 天卖出，能赚多少。买入的数量是 $A$ 券 $\frac{r_j}{a_j{r}_j+b_j}\cdot money$ 与 $B$ 券 $\frac{1}{a_j{r}_j+b_j}\cdot money$ ，所以转移式
$$f(i)=\underset{j=1}{\overset{i-1}{\max }}\left[\frac{a_i{r}_j+b_i}{a_j{r}_j+b_j}\cdot f(j)\right]$$

记 $v_j=\frac{f(j)}{a_j{r}_j+b_j},p_j=r_j{v}_j$ 则
$$f(i)=b_i\cdot \underset{j=1}{\overset{i-1}{\max }}\left(p_j\cdot \frac{a_i}{b_i}+v_j\right)$$

这就是一个李超树板题了嘛。复杂度 $\mathcal{O}(n\log n)$ 。

代码

一开始我傻了，直接用 $\mathtt{d}\mathtt{o}\mathtt{u}\mathtt{b}\mathtt{l}\mathtt{e}$，结果精度爆炸。事实上，查询点 $\frac{a_i}{b_i}$ 是固定的，所以可以离散。

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long int_;
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
inline void writeint(int_ x){
	if(x > 9) writeint(x/10);
	putchar((x-x/10*10)^48);
}

struct Line{
	double k, b; int id;
	Line(){ k = b = 0; }
	double operator()(const double &x) const {
		return k*x+b; // definition
	}
};
bool low(const Line &A,const Line &B,const double &x){
	return (A.k-B.k)*x < B.b-A.b;
}

const int MaxN = 500005;
int son[MaxN][2], cntNode;
Line line[MaxN]; double all[MaxN];
void modify(Line x,int &o,int l,int r){
	if(!o) o = ++ cntNode; // new node
	int mid = (l+r)>>1; // all[mid]
	if(low(line[o],x,all[l])) swap(x,line[o]);
	if(low(x,line[o],all[r])) return ;
	if(low(x,line[o],all[mid]))
		modify(x,son[o][1],mid+1,r);
	else{
		modify(line[o],son[o][0],l,mid);
		line[o] = x; // a better ans
	}
}

double dp[MaxN], a[MaxN], b[MaxN], R[MaxN];
double query(int qid,int real,int o,int l,int r){
	double res = 0;
	for(int d; o; o=son[o][d]){
		res = max(res,
			(a[real]*R[line[o].id]
			+b[real])*dp[line[o].id]
			/(a[line[o].id]*R[line[o].id]
			+b[line[o].id])
		);
		if((d = (qid*2 > l+r)))
			l = (l+r)/2+1;
		else r = (l+r)>>1;
	}
	return res;
}

int main(){
	int n = readint(); dp[0] = readint();
	for(int i=1; i<=n; ++i){
		scanf("%lf %lf %lf",&a[i],&b[i],&R[i]);
		all[i] = a[i]/b[i]; // all needed X
	}
	sort(all+1,all+n+1); Line now;
	int tot = unique(all+1,all+n+1)-all;
	for(int i=1,rt=0; i<=n; ++i){
		int haxi = lower_bound(all+1,all+tot,a[i]/b[i])-all;
		dp[i] = max(dp[i-1],query(haxi,i,rt,1,tot-1));
		now.b = dp[i]/(a[i]*R[i]+b[i]); // v_j
		now.k = now.b*R[i], now.id = i;
		modify(now,rt,1,tot-1); // add i
	}
	printf("%.3f\n",dp[n]);
	return 0;
}