BZOJ1492 [NOI2007]货币兑换

Address

• BZOJ1492 洛谷P4027

Solution

• 因为傻逼错误调了一下午。
• 有一个贪心策略：每次买入或卖出时都是全部买入或全部卖出。
• 于是分别设 $f_i,x_i,y_i$ 表示到第 $i$ 天时的最大收益、最多能购买的 $A$ 劵数量、最多能购买的 $B$ 劵数量，则转移为： f i = max ⁡ j = 1 i − 1 { f i − 1 , A i × x j + B i × y j } f_i = \max\limits^{i - 1}_{j = 1}\{f_{i - 1}, A_i \times x_j + B_i \times y_j\} fi​=j=1maxi−1​{fi−1​,Ai​×xj​+Bi​×yj​}$$y_i=\frac{f_i}{A_i\times Rat{e}_i+B_i}$$$$x_i=Rat{e}_i\times y_i$$
• 时间复杂度 $O(N^2)$，不能通过。
• $x_i,y_i$ 只由 $f_i$ 确定，考虑怎样优化 $f_i$ 的转移。
• 先不管转移中的 $f_{i-1}$ 和 $max$，则有：$$f_i=A_i\times x_j+B_i\times y_j$$
• 选择两个含 $j$ 的不同项，分别移到等式两边，并使右边含 $j$ 的项系数为1，把只含 $i$ 的项移到左边：$$-\frac{A_i}{B_i}\times x_j+\frac{f_i}{B_i}=y_j$$
• 因为要使 $f_i$ 最大，即为让一条斜率为 $-\frac{A_i}{B_i}$ 的直线经过点 $(x_j,y_j)$，使得它的截距（$\frac{f_i}{B_i}$）最大。
• 通过画图我们可以发现，最优解必然在点 $(x_j,y_j)$ 构成的上凸壳上（如果取最小就是下凸壳），并且最优解的 $j$ 即满足它与它在凸壳上的上一个点所在直线的斜率大于 $-\frac{A_i}{B_i}$，它与它在凸壳上的下一个点所在直线的斜率小于 $-\frac{A_i}{B_i}$。
• 然而这题 $x_j,y_j$ 都不递增，不能用单调队列或者二分来维护。
• $Splay$ 大法好！
• 按 $x_j$ 的大小在 $Splay$ 上排好顺序，记录 $l{k}_i,r{k}_i$分别表示与凸壳上一个点和下一个点所在直线的斜率。
• 查找最优解只要在 $Splay$ 上走，主要要解决的是如何动态插入一个点 $x$。
• 同样考虑二分，能作为 $x$ 在凸壳上的上一个点 $y$ 要满足这两点所在直线的斜率小于 $l{k}_y$，并且 $y$ 要取到最大，$x$ 在凸壳上的下一个点同理。
• 若插入后 $l{k}_x<r{k}_x$，点 $x$ 在凸壳内，直接把点 $x$ 删除。
• 注意精度问题，时间复杂度 $O(均摊n\log n)$。

Code

• 数据中似乎不用考虑点在凸壳内的情况。
• 两个地方把 x 写成 y 还有 90pts。
• 这数据到底有多水……

#include <iostream>
#include <cstdio>
#include <cctype>
#include <algorithm>
#include <cstring>
#include <cstdlib>

using namespace std;

const double eps = 1e-9;
const double feps = -1e-9;
const double inf = 1e15;
const double finf = -1e15;
const int N = 1e5 + 5;
double lk[N], rk[N], f[N], vx[N], vy[N];
int n, T, rt, lc[N], rc[N], fa[N];

inline double Max(double x, double y) {return x + eps >= y ? x : y;}

inline bool Which(int x) {return lc[fa[x]] == x;}

inline void Rotate(int x)
{
	int y = fa[x], z = fa[y],
		b = lc[y] == x ? rc[x] : lc[x];
	if (z) (lc[z] == y ? lc[z] : rc[z]) = x;
	fa[x] = z; fa[y] = x;
	if (b) fa[b] = y;
	if (lc[y] == x) rc[x] = y, lc[y] = b;
		else lc[x] = y, rc[y] = b;
}

inline void Splay(int x, int tar)
{
	while (fa[x] != tar)
	{
		if (fa[fa[x]] != tar)
			Which(fa[x]) == Which(x) ? Rotate(fa[x]) : Rotate(x);
		Rotate(x);
	}
	if (!tar) rt = x;
}

inline double Slope(int x, int y)
{
	if (vx[x] - vx[y] <= eps && vx[x] - vx[y] >= feps) 
		return finf;
	return (vy[y] - vy[x]) / (vx[y] - vx[x]);
}

inline int findLeft(int x, int y)
{
	int res = x;
	while (x)
	{
		if (lk[x] + eps >= Slope(x, y)) res = x, x = rc[x];
			else x = lc[x];
	}
	return res;
}

inline int findRight(int x, int y)
{
	int res = x;
	while (x)
	{
		if (rk[x] <= Slope(y, x) + eps) res = x, x = lc[x];
			else x = rc[x];
	}
	return res;
}
 
inline void Insert()
{
	int x = rt, y = 0, dir; ++T;
	while (x)
	{
		y = x;
		if (vx[T] <= vx[x] + eps) x = lc[x], dir = 0;
			else x = rc[x], dir = 1;
	}
	fa[x = T] = y;
	if (y) (dir ? rc[y] : lc[y]) = x;
	Splay(x, 0); 
	if (lc[x])
	{
		int z = findLeft(lc[x], x);
		Splay(z, x); fa[rc[z]] = rc[z] = 0; 
		rk[z] = lk[x] = Slope(z, x); 
	}
	
	if (rc[x])
	{
		int z = findRight(rc[x], x);
		Splay(z, x); fa[lc[z]] = lc[z] = 0;
		lk[z] = rk[x] = Slope(x, z);
	}
	
	if (lk[x] <= rk[x] + eps)
	{
		rt = lc[x]; rc[rt] = rc[x];
		fa[rc[x]] = rt; fa[rt] = 0;
		lc[x] = rc[x] = fa[x] = 0;
		rk[rt] = lk[rc[rt]] = Slope(rt, rc[rt]);
	}
}

inline int Query(double v)
{
	int x = rt, res = 0;
	while (x)
	{
		if (lk[x] + eps >= v) res = x, x = rc[x];
			else x = lc[x];
	}
	return res;
}

int main()
{
	scanf("%d%lf", &n, &f[0]); double a, b, r;
	for (int i = 1; i <= n; ++i)
		lk[i] = inf, rk[i] = finf;
	
	for (int i = 1; i <= n; ++i)
	{ 
		scanf("%lf%lf%lf", &a, &b, &r);
		int j = Query(-a / b);
		f[i] = Max(f[i - 1], a * vx[j] + b * vy[j]);
		vy[i] = f[i] / (a * r + b);
		vx[i] = vy[i] * r; Insert();
	}
	
	printf("%.3lf\n", f[n]);
}