bzoj-1492 货币兑换Cash (1)——平衡树维护凸包

原创 2015年07月11日 09:37:45

题意:

有n天和m的初始金钱,用来购买AB两种纪念券;

n天里每天都有AB的价格,每天可以进行这样的操作;

1.卖出手中x%的纪念券(AB分别都卖出x%);

2.用x的金钱买入纪念券,买入AB券的比例在第i天为Rate i;

求n天过去之后所获得的最大收益;

金钱和券数均为实数;

n<=100 000;


题解:

首先,虽然题中的买入和卖出都是任意数量的,但是同样的纪念券,分几天卖出得到的收 益,一定小于等于直接在一天卖出的收益;

同样,分几天买入也是不如一天花所有钱买入的;

令:

f[i]为第i天的最大收益;

X[i]为第i天将所有钱数都买入得到的A券数;

Y[i]为第i天将所有钱数都买入得到的B券数;

显然X,Y都可以由f[i]得来;

那么转移方程就是:

f[i]=max(f[i-1],A[i] * X[j] + B[i] * Y[j]);(1<=j< i);

这样转移是O(n^2)的,所以要对后面枚举j的部分优化;

将f[i]=A[i] * X[j] + B[i] * Y[j]整理;

得到Y[i]=(-A[i]/B[i]) * X[i] + f[i]/B[i];

这是一个直线方程的形式,而对于固定的i,直线斜率不变,而要最大化截距;

倘若将所有的[1,i-1]的点计算出(x,y)放在坐标系上;

找到最优值相当于在这个上凸包上找到某个点,使截距最大;

那么这个点左面的斜率一定大于当前i的斜率,右面的斜率一定小于当前i的斜率;

所以这其实就是一个斜率优化的形式;

通常我们做斜率优化都是找到不符合要求的点直接干掉就好的;

因为下一个i的斜率不是递增就是递减;

但是这个-A[i]/B[i]不单调,所以不能O(n)的维护队列处理凸包;

10^5的复杂度是支持O(nlogn)的,所以为了维护凸包可以选择一些数据结构;

那么就维护一个Splay,每个结点都在凸包上,中序遍历就是按x递增同时也按斜率递减的序列;

如果把Splay看做logn,那么每个点最多进出凸包一次,二分查询斜率共n次;

复杂度O(nlogn)还是听起来很好的;

但是写起来一点也不好玩!

总之维护凸包时对Splay中点较少时的讨论很烦。。。

照着对拍调数据改改也就过了;

代码3k+,时间960ms,这个跑的感觉也是挺快的了;

下一篇写写更神的CDQ分治,毕竟数据结构对代码能力要求颇高;

7/17

我被D了。。。斜率打成了shope,恩应该是slope无误;


代码:


#include<math.h>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 110000
#define which(x)	(tr[tr[x].fa].ch[1]==x)
const double INF = 1e100;
const double EPS = 1e-8;
using namespace std;
struct Point
{
	double x, y, s1, s2;
	int fa, ch[2];
}tr[N];
int root, tot;
double f[N], A[N], B[N], R[N], X[N], Y[N];
void get_slope(int a, int b)
{
	if (!a)			tr[b].s1 = INF;
	else if (!b)		tr[a].s2 = -INF;
	else
	{
		if (fabs(tr[a].x - tr[b].x)<EPS)
			tr[a].s2 = tr[b].s1 = (tr[a].y<tr[b].y ? INF : -INF);
		else
			tr[a].s2 = tr[b].s1 = (tr[a].y - tr[b].y) / (tr[a].x - tr[b].x);
	}
}
void Rotate(int x)
{
	int f = tr[x].fa;
	if (!f)	return;
	bool k = which(x);
	tr[f].ch[k] = tr[x].ch[!k];
	tr[x].ch[!k] = f;
	tr[tr[f].fa].ch[which(f)] = x;
	tr[x].fa = tr[f].fa;
	tr[tr[f].ch[k]].fa = f;
	tr[f].fa = x;
}
void Splay(int x, int g)
{
	if (!x)	return;
	while (tr[x].fa != g)
	{
		int f = tr[x].fa;
		if (tr[f].fa == g)
		{
			Rotate(x);
			break;
		}
		if (which(x) ^ which(f))
			Rotate(x);
		else
			Rotate(f);
		Rotate(x);
	}
	if (!g)	root = x;
}
int Pre(int x)
{
	if (!x)	return 0;
	int p = tr[x].ch[0];
	if (!p)	return 0;
	while (tr[p].ch[1])
		p = tr[p].ch[1];
	return p;
}
int Sub(int x)
{
	if (!x)	return 0;
	int p = tr[x].ch[1];
	if (!p)	return 0;
	while (tr[p].ch[0])
		p = tr[p].ch[0];
	return p;
}
int find(int p, double x)
{
	if (!p)	return 0;
	if (x<tr[p].x)
		return find(tr[p].ch[0], x);
	else
	{
		int t = find(tr[p].ch[1], x);
		return tr[p].x>tr[t].x ? p : t;
	}
}
void Insert(double X, double Y, int no)
{
	int x = find(root, X), y = 0;
	if (!x)
	{
		x = root;
		while (tr[x].ch[0])
			x = tr[x].ch[0];
		Splay(x, 0);
		y = x, x = 0;
	}
	else
	{
		Splay(x, 0);
		Splay(y = Sub(x), x);
	}
	tr[no].x = X, tr[no].y = Y;
	if (y)	tr[no].fa = y, tr[y].ch[0] = no;
	else	tr[no].fa = x, tr[x].ch[1] = no;
	get_slope(x, no);
	get_slope(no, y);
	if (tr[no].s1 <= tr[no].s2)
	{
		tr[y].ch[0] = 0;
		get_slope(x, y);
		return;
	}
	Rotate(no), Rotate(no);
	root = no;
	x = tr[no].ch[0];
	while (tr[x].s1 <= tr[x].s2&&x)
	{
		y = Pre(x);
		Splay(y, x);
		tr[y].fa = no;
		tr[no].ch[0] = y;
		get_slope(y, no);
		x = y;
	}
	x = tr[no].ch[1];
	while (tr[x].s1 <= tr[x].s2&&x)
	{
		y = Sub(x);
		Splay(y, x);
		tr[y].fa = no;
		tr[no].ch[1] = y;
		get_slope(no, y);
		x = y;
	}
}
int query(double S)
{
	int p = root;
	while (S>tr[p].s1 || S<tr[p].s2)
	{
		if (S>tr[p].s1)	p = tr[p].ch[0];
		else			p = tr[p].ch[1];
	}
	return p;
}
int main()
{
	int n, i, j, k;
	scanf("%d%lf", &n, &f[1]);
	for (i = 1; i <= n; i++)
		scanf("%lf%lf%lf", A + i, B + i, R + i);
	tr[0].x = tr[0].y = -INF;
	Y[1] = f[1] / (A[1] * R[1] + B[1]);
	X[1] = R[1] * Y[1];
	Insert(X[1], Y[1], 1);
	for (i = 2; i <= n; i++)
	{
		j = query(-A[i] / B[i]);
		f[i] = max(f[i - 1], A[i] * X[j] + B[i] * Y[j]);
		Y[i] = f[i] / (A[i] * R[i] + B[i]);
		X[i] = R[i] * Y[i];
		Insert(X[i], Y[i], i);
	}
	printf("%.3lf", f[n]);
	return 0;
}


版权声明:

相关文章推荐

BZOJ 1492 货币兑换 Cash 平衡树维护凸包 CDQ分治

题目大意: 这题真是为难了我1B。。。做了三天,两种方法都试过了一遍,真是膜拜NOI的神犇们能在一个多小时内干掉这道题 首先这题一看就是斜率优化 连递推式都是P=A[i]*X[i]+B...

BZOJ 1492 Cash(CDQ分治论文题)

CDQ分治论文题 #include #include #include #include #include #include #include #include #include #include #...

bzoj 1492: [NOI2007]货币兑换Cash (dp+cdq分治)

1492: [NOI2007]货币兑换Cash Time Limit: 5 Sec  Memory Limit: 64 MB Submit: 3801  Solved: 1602 [Submit...

bzoj1492(cdq分治+凸包+优化dp)

我们会发现dp不满足单调性,但是可以构建凸包,其中点的表示方法很巧妙   要对一些问题转化成凸包敏感一点,尝试找出表示凸包的点,看看是否可以用斜率来表示答案。   一定要写eps,wa到死。。...

bzoj1492: [NOI2007]货币兑换Cash

货币兑换   问题描述      小 Y 最近在一家金券交易所工作。该金券交易所只发行交易两种金券:A 纪  念券(以下简称 A 券)和 B 纪念券(以下简称 B 券)。每个持有金券的顾客都有  ...

bzoj1492 货币兑换 CDQ分治

CDQ分治。。用来解决不满足斜率优化条件的dp。        首先,原题的Hint有一条就是要么全部买入,要么全部卖出。否则不最优。这是显然的。        那么令f[i]表示第i天的最大获利,f...

1492: [NOI2007]货币兑换Cash DP+斜率优化+splay维护凸包/CDQ分治

似乎两个月之前我写过一次这题写挂了? 《多年的心头大恨终于解决了系列》 我们可以发现,存在最优解满足每次将全部的钱买了券或每次将全部的券卖成钱,因为有利益我们要尽可能的去获得,有亏损就一点也不要碰...

[挖坑]待学的算法

我太弱啦= =一大波东西没学,有很多东西还是有必要学一下的。可并堆 动态树分治 树套树 可持久化线段树 可持久化trie 可持久化平衡树 虚树 基环树 分块 莫队算法 块状链表 ...

NOI2016滚粗记

Day-1由于飞机是明天一早,所以今天就要先去济南了。下午2:00出发,与之前不同的是这次是商务车直达。路上老师检查了一下笔试,似乎没什么问题。傍晚就到了,老师们请客一起吃饭,然而在机场旁太偏了并没什...

bzoj 2300: [HAOI2011]防线修建(splay动态维护凸包)

2300: [HAOI2011]防线修建 Time Limit: 10 Sec  Memory Limit: 512 MB Submit: 800  Solved: 440 [Submit][S...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)