BZOJ1492：[NOI2007]货币兑换 （CDQ分治+斜率优化DP | splay动态维护凸包）

最新推荐文章于 2019-10-17 16:01:58 发布
最新推荐文章于 2019-10-17 16:01:58 发布
阅读量396 收藏 1
点赞数 1
分类专栏： cdq分治 Splay 文章标签： 动态凸包 cdq分治 斜率优化DP
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/kekxy/article/details/77854648
版权

BZOJ1492：[NOI2007]货币兑换

题目传送门

【问题描述】

小Y最近在一家金券交易所工作。该金券交易所只发行交易两种金券：A纪念券（以下简称A券）和B纪念券（以下简称B券）。每个持有金券的顾客都有一个自己的 帐户。金券的数目可以是一个实数。每天随着市场的起伏波动，两种金券都有自己当时的价值，即每一单位金券当天可以兑换的人民币数目。我们记录第K天中A券 和B券的价值分别为AK和BK（元/单位金券）。
为了方便顾客，金券交易所提供了一种非常方便的交易方式：比例交易法。
比例交易法分为两个方面：
a)卖出金券：顾客提供一个[0，100]内的实数OP作为卖出比例，其意义为：将OP%的A券和OP%的B券以当时的价值兑换为人民币；
b)买入金券：顾客支付IP元人民币，交易所将会兑换给用户总价值为IP的金券，并且，满足提供给顾客的A券和B券的比例在第K天恰好为RateK；
例如，假定接下来3天内的Ak、Bk、RateK的变化分别为：

时间akbkrk
第一天111
第二天122
第三天223

假定在第一天时，用户手中有100元人民币但是没有任何金券。
用户可以执行以下的操作：

时间用户操作人民币(元)A券的数量B券的数量
开户无10000
第一天买入100元05050
第二天卖出50%752525
第二天买入60元155540
第三天卖出100%20500

注意到，同一天内可以进行多次操作。
小Y是一个很有经济头脑的员工，通过较长时间的运作和行情测算，他已经知道了未来N天内的A券和B券的价值以及Rate。他还希望能够计算出来，如果开始时拥有S元钱，那么N天后最多能够获得多少元钱。
【输入格式】
第一行两个正整数N、S，分别表示小Y能预知的天数以及初始时拥有的钱数。接下来N行，第K行三个实数AK、BK、RateK，意义如题目中所述。
【输出格式】
只有一个实数MaxProfit，表示第N天的操作结束时能够获得的最大的金钱数目。答案保留3位小数。
【输入样例】
3 100
1 1 1
1 2 2
2 2 3
【输出样例】
225.000
【样例说明】

时间用户操作人民币(元)A券的数量B券的数量
开户无10000
第一天买入100元05050
第二天卖出100%15000
第二天买入150元07537.5
第三天卖出100%22500

【数据规模和约定】
输入第一行两个正整数N、S，分别表示小Y能预知的天数以及初始时拥有的钱数。接下来N行，第K行三个实数AK、BK、RateK，意义如题目中所述。对于100%的测试数据，满足：0＜AK≤10；0＜BK≤10；0＜RateK≤100；MaxProfit≤10^9。
【提示】
1.输入文件可能很大，请采用快速的读入方式。
2.必然存在一种最优的买卖方案满足：
每次买进操作使用完所有的人民币；
每次卖出操作卖出所有的金券。

【解题思路】

有点贪心的想法：一有亏损就不去碰，一有利益就去全占。这必然导致要不全部卖出，要不全部买入。

①CDQ分治+斜率优化DP

O(N^2)的DP很好想。F[n]表示到第n天，没有金券可以得到的最大价值。
X[n]表示第n天的钱买金券可以得到的A劵数目；
Y[n]表示第n天的钱买金券可以得到的B劵数目；
Y[n]=F[n]A[n]∗R[n]+B[n]

Y[n]=X[n]∗R[n]

F[n]=Max(X[i]∗A[n]+Y[i]∗B[n],F[n−1])
假如X[j]＜X[k]，并且k比j优。那么
X[k]∗A[n]+Y[k]∗B[n]>X[j]∗A[n]+Y[j]∗B[n]
转化一下式子：
Y[k]−Y[j]X[k]−X[j]>−A[n]B[n]
slop(j,k)=Y[k]−Y[j]X[k]−X[j]
如果k比j优，则slop(j,k)>-A[n]/B[n]
单调队列加入F[i]时，当slop(q[qe-1],q[qe])＜slop(q[qe],i)，那么q[qe]是没必要存的。
假如-A[n]/B[n]是从大到小的。
∵若slop(q[qe],i)>-A[n]/B[n]，则i比q[qe]优，剔除q[qe]。
若slop(q[qe],i)＜-A[n]/B[n]，则q[qe]比i优，但slop(q[qe-1],q[qe])＜-A[n]/B[n]，q[qe-1]比q[qe]优，剔除q[qe]。
但这不能直接用斜率优化。X[i]和-A[i]/B[i]不一定同时有序。
CDQ分治可以很好的解决这一问题。
CDQ(L,R)表示处理id为L到R的F[]。
先按-A[i]/B[i]从大到小排序。
CDQ(1,n)
1、处理id在L~Mid的F[i]，CDQ(L,Mid)
2、用id在L~Mid的F[i]，更新id在Mid+1~R的F[i]。
有可能对于F[i]，最优解在F[j]处取得（Mid+1≤i≤N，1≤j≤Mid）
（因为L~Mid的X[i]是有序的，Mid+1~R的-A[i]/B[i]是有序的，满足斜率优化的条件）
3、处理id在Mid+1~R的F[i]，CDQ(Mid+1,R)。
有可能对于F[i]，最优解在F[j]处取得（Mid+1≤i≤N，Mid+1≤j＜i）
4、对id在L~R，按X[i]从小到大排序（类似归并）。
【代码】

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>

using namespace std;

const int N=100010;
const double inf=1e20;
const double eps=1e-8;
double f[N];
struct data
{
    double x,y,a,b,k,r;
    int id;
} p[N],t[N];
int n,qe,q[N];

double slop(int a,int b)
{
    if(!b) return -inf;
    if(fabs(p[a].x-p[b].x)<eps) return inf;
    return ((p[b].y-p[a].y)/(p[b].x-p[a].x));
}

bool cmp(data A,data B) { return (A.k>B.k); }

void CDQ(int L,int R)
{
    if(L==R)
    {
        f[L]=max(f[L],f[L-1]);
        p[L].y=f[L]/(p[L].a*p[L].r+p[L].b);
        p[L].x=p[L].y*p[L].r;
        return;
    }
    int Mid=(L+R)>>1;
    int l1=L,l2=Mid+1; int h=1;
    for(int i=L;i<=R;i++)
    if(p[i].id<=Mid) t[l1++]=p[i]; else t[l2++]=p[i];
    for(int i=L;i<=R;i++) p[i]=t[i];
    CDQ(L,Mid);
    qe=0;
    for(int i=L;i<=Mid;i++)
    {
        while(qe>1 && slop(q[qe-1],q[qe])<slop(q[qe],i)+eps) qe--;
        q[++qe]=i;
    }
    q[++qe]=0;
    for(int i=Mid+1;i<=R;i++)
    {
        while(h<qe && slop(q[h],q[h+1])+eps>p[i].k) h++;
        f[p[i].id]=max(f[p[i].id],p[q[h]].x*p[i].a+p[q[h]].y*p[i].b);
    }
    CDQ(Mid+1,R);
    l1=L; l2=Mid+1;
    for(int i=L;i<=R;i++)
    {
        if(((p[l1].x<p[l2].x || (fabs(p[l1].x-p[l2].x)<eps && p[l1].y<p[l2].y)) || l2>R) && l1<=Mid) t[i]=p[l1++];
        else t[i]=p[l2++];
    }
    for(int i=L;i<=R;i++) p[i]=t[i];
}

int main()
{
    freopen("cash.in","r",stdin);
    freopen("cash.out","w",stdout);
    scanf("%d%lf",&n,&f[0]);
    for(int i=1;i<=n;i++)
    {
        scanf("%lf%lf%lf",&p[i].a,&p[i].b,&p[i].r);
        p[i].id=i; p[i].k=-p[i].a/p[i].b;
    }
    sort(p+1,p+1+n,cmp);
    CDQ(1,n);
    printf("%.3lf\n",f[n]);
    return 0;
}

②splay动态维护凸包

F[n]=X[i]∗A[n]+Y[i]∗B[n]
转化一下式子：
Y[i]=−A[n]B[n]X[i]+F[n]B[n]
其实就是使一次函数的截距最大（B[N]为定值，不用理）。
这里写图片描述
F[n]B[n]=Y[i]+A[n]B[n]X[i]
相当于一条斜率为-A[n]/B[n]的直线从上往下扫，碰到的第一个点(Xj,Yj)就是最优解（j＜n）
这里写图片描述
只需要维护一个凸包就行了，最优解只可能在凸包上，凸包内的点不会影响答案。
记得特判……
这里写图片描述
【代码】

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>

using namespace std;

const double eps=1e-8;
const double inf=1e9;
const int N=100010;
int n,root;
double f[N],A[N],B[N],r[N],X[N],Y[N];
double lk[N],rk[N];
int fa[N],ch[N][2],D[N];

int son(int x) { return (ch[fa[x]][1]==x); }

void rotate(int x)
{
    int y=fa[x],z=son(x);
    fa[x]=fa[y];
    if(fa[y]) ch[fa[y]][son(y)]=x;
    if(ch[x][1-z]) fa[ch[x][1-z]]=y;
    fa[y]=x; ch[y][z]=ch[x][1-z]; ch[x][1-z]=y;
}

void splay(int x,int y)
{
    for(;fa[x]!=y;)
    {
        if(fa[fa[x]]!=y)
        (son(x)==son(fa[x]))?rotate(fa[x]):rotate(x);
        rotate(x);
    }
    if(!y) root=x;
}

double getk(int i,int j)
{
    if(fabs(X[j]-X[i])<eps) return -inf;
    return((Y[j]-Y[i])/(X[j]-X[i]));
}

int pre(int x)
{
    int y=ch[x][0],z=x;
    for(;y;) if(getk(y,x)<=lk[y]+eps) z=y,y=ch[y][1]; else y=ch[y][0];
    return z;
}

int suc(int x)
{
    int y=ch[x][1],z=x;
    for(;y;) if(getk(x,y)+eps>=rk[y]) z=y,y=ch[y][0]; else y=ch[y][1];
    return z;
}

void inse(int &v,int f,int id)
{
    if(!v) { v=id; fa[v]=f; splay(v,0); return; }
    if(X[id]<=X[v]+eps) inse(ch[v][0],v,id); else inse(ch[v][1],v,id);
}

int find(int v,double k)
{
    if(!v) return 0;
    if(lk[v]+eps>=k && k+eps>=rk[v]) return v;
    if(k>lk[v]) return find(ch[v][0],k); else return find(ch[v][1],k);
}

void updata(int x)
{
    splay(x,0);
    if(ch[x][0])
    {
        int lef=pre(x); splay(lef,x); ch[lef][1]=0;
        lk[x]=rk[lef]=getk(lef,x);
    } else lk[x]=inf;
    if(ch[x][1])
    {
        int rig=suc(x); splay(rig,x); ch[rig][0]=0;
        lk[rig]=rk[x]=getk(x,rig);
    } else rk[x]=-inf;
    if(lk[x]<=rk[x]+eps)
    {
        root=ch[x][0]; ch[root][1]=ch[x][1];
        fa[ch[x][1]]=root; fa[x]=0;
        rk[root]=lk[ch[root][1]]=getk(root,ch[root][1]);
    }
}

int main()
{
    freopen("cash.in","r",stdin);
    freopen("cash.out","w",stdout);
    scanf("%d%lf",&n,&f[0]);
    for(int i=1;i<=n;i++) scanf("%lf%lf%lf",&A[i],&B[i],&r[i]);
    root=0;
    for(int i=1;i<=n;i++)
    {
        double yu=-A[i]/B[i];
        int now=find(root,yu);
        f[i]=max(f[i-1],X[now]*A[i]+Y[now]*B[i]);
        Y[i]=f[i]/(r[i]*A[i]+B[i]);
        X[i]=Y[i]*r[i];
        inse(root,0,i);
        updata(i);
    }
    printf("%.3lf\n",f[n]);
    return 0;
}
kekxy
关注
专栏目录
1492: [NOI2007]货币兑换Cash DP+斜率优化+splay维护凸包/CDQ分治
ws_fqk
02-29 993
似乎两个月之前我写过一次这题写挂了？ 《多年的心头大恨终于解决了系列》 我们可以发现，存在最优解满足每次将全部的钱买了券或每次将全部的券卖成钱，因为有利益我们要尽可能的去获得，有亏损就一点也不要碰。 我们用fif_i表示到ii天时的最大收益，令xi,yix_i,y_i分别表示第ii天购买AA券和BB券的数量。则如果将第ii天的收益全部买成券，那么满足 fi=Ai∗xi+Bi∗yif_i=A_
BZOJ3963: [WF2011]MachineWorks 【CDQ+斜率优化DP】*
Dream_maker的博客
09-02 243
BZOJ3963: [WF2011]MachineWorks Description 你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理，公司使用更加先进的机械设备生产先进的机器。原来的那一台生产机器已经坏了，所以你要去为公司买一台新的生产机器。你的任务是在转型期内尽可能得到更大的收益。在这段时间内，你要买卖机器，并且当机器被ACM公司拥...
[NOI2007]货币兑换Cash（DP+动态凸包）
weixin_30483697的博客
03-20 80
第一次打动态凸包维护dp，感觉学到了超级多的东西。 首先，set是如此的好用！！！可以通过控制一个flag来实现两种查询，维护凸包和查找斜率k 不过就是重载运算符和一些细节方面有些恶心，90行解决 后面还有一个cdq分治，找时间学下，看下能不能处理一大类恶心的问题 github还是不会用，找时间搞下吧 CODE： 1 #include<cstdio> 2 #inc...
[DP 斜率优化 CDQ分治||动态维护凸包] BZOJ 1492 [NOI2007]货币兑换Cash
雯舞
08-08 683
打了个set维护凸包 cdq等待填坑 #include #include #include #include #include #include #define dprintf(...) fprintf(stderr,__VA_ARGS__) using namespace std; //typedef long double ld; typedef double ld; inline
bzoj 1492: [NOI2007]货币兑换Cash
lqybzx的专栏
04-13 1036
Description Input 第一行两个正整数N、S，分别表示小Y 能预知的天数以及初始时拥有的钱数。 接下来N 行，第K 行三个实数AK、BK、RateK，意义如题目中所述 Output 只有一个实数MaxProfit，表示第N 天的操作结束时能够获得的最大的金钱 数目。答案保留3 位小数。 Sample Input 3 100 1 1 1 1 2 2
【bzoj1492】[NOI2007]货币兑换Cash 斜率优化+set+凸包
Oxer的专栏
03-16 764
斜率优化boss题？ 首先这道题有一个性质，一定是在一天买进一天卖出，而且全都是倾巢买入卖出 先写方程 f[i]表示到前i天最多获得多少元 f[i]=max{f[j]/(a[j]*rate[j]+b[j])*rate[j]*a[i]+f[j]/(a[j]*rate[j]+b[j])*b[i]} (1 f[i]=max{x[j]*a[i]+y[j]*b[i]} =(x[j]*a[i]/
【BZOJ3963】【ACM-WF2011】MachineWorks（CDQ分治+斜率优化）
(■-■)
03-17 1275
Description你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理，公司使用更加先进的机械设备生产先进的机器。原来的那一台生产机器已经坏了，所以你要去为公司买一台新的生产机器。你的任务是在转型期内尽可能得到更大的收益。在这段时间内，你要买卖机器，并且当机器被ACM公司拥有的时候，操控这些机器以获取利润。因为空间的限制，ACM公司在任何时候都只能最
bzoj2726: [SDOI2012]任务安排 斜率优化Dp+cdq分治
maxtir的博客
07-16 327
bzoj2726: [SDOI2012]任务安排 Description 机器上有N个需要处理的任务，它们构成了一个序列。这些任务被标号为1到N，因此序列的排列为1,2,3…N。这N个任务被分成若干批，每批包含相邻的若干任务。从时刻0开始，这些任务被分批加工，第i个任务单独完成所需的时间是Ti。在每批任务开始前，机器需要启动时间S，而完成这批任务所需的时间是各个任务需要时间的总和。注意，同...
bzoj3451: Tyvj1953 Normal（点分治+ntt）
苟为蒟蒻又何妨
08-22 176
传送门 利用期望的线性性。 考虑到iii在jjj子树中的概率为1dist(i,j)\frac1{dist(i,j)}dist(i,j)1​，于是问题变成了统计两点间距离的出现次数。 代码： #include<bits/stdc++.h> #define ri register int using namespace std; const int rlen=1<<18|1; ...
BZOJ 4311: 向量（线段树分治+凸包+三分）
Ab.Ever
11-05 573
Description你要维护一个向量集合，支持以下操作： 1.插入一个向量(x,y) 2.删除插入的第i个向量 3.查询当前集合与(x,y)点积的最大值是多少。如果当前是空集输出0Input第一行输入一个整数n，表示操作个数 接下来n行，每行先是一个整数t表示类型，如果t=1，输入向量 (x,y)；如果t=2，输入id表示删除第id个向量；否则输入(x,y)，查询 与向量(x,y)点积
bzoj 2300 [HAOI2011]防线修建 splay维护凸包
make_it_for_good的博客
12-28 405
正常的splay维护凸包不支持删点，把操作倒过来变成加点。 HAOI真心良心呀，边界都给好了。。。#include <bits/stdc++.h> using namespace std; #define N 110000 #define M 210000 #define which(x) (ch[fa[x]][1]==x) int d,n,Q,root; int X[N],Y[N],vis[N]
bzoj 2300: [HAOI2011]防线修建（splay动态维护凸包）
clover_hxy的博客
09-10 544
2300: [HAOI2011]防线修建 Time Limit: 10 Sec  Memory Limit: 512 MB Submit: 800  Solved: 440 [Submit][Status][Discuss] Description 近来A国和B国的矛盾激化，为了预防不测，A国准备修建一条长长的防线，当然修建防线的话，肯定要把需要保护的城市修在防线内部了。可是A国上层
bzoj3963 [WF2011]MachineWorks（斜率优化 & splay维护动态凸包）
Bfk的博客
01-15 480
bzoj3963 [WF2011]MachineWorks 原题地址：http://www.lydsy.com/JudgeOnline/problem.php?id=3963 题意： 多组数据。 你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理，公司使用更加先进的机械设备生产先进的机器。原来的那一台生产机器已经坏了，所以你要去为公司买一台新
[BZOJ1492][NOI2007]货币兑换 斜率优化DP+splay维护凸壳
DOFYPXY的博客
04-26 727
写着题写的想死。。。来扯扯写着题的过程。。。（我菜啊） 早晨总结了一下昨天写的几道决策单调性和斜率优化dp，然后突然看到这题，准备来写，此时8:30。 看了popoQQQ的题解，感觉这题代码有点难度，不过还是要挑战一下的，然后大概理了理思路和细节，就开始码了，此时9:00。 10:40码完，调了几个splayRE的地方，11:00过了样例（样例很水啊），交一发，TLE（后来分析应该是死循环了）
BZOJ 3963: [WF2011]MachineWorks 斜率优化 + splay动态维护凸包
EM LGH
07-08 146
Code: #include<bits/stdc++.h> #define setIO(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout) #define maxn 300000 #define inf 10000000000000009 #define ll long long ...
3963: [WF2011]MachineWorks
CRZbulabula的博客
11-10 475
3963: [WF2011]MachineWorks Time Limit: 30 Sec  Memory Limit: 256 MB Submit: 93  Solved: 28 [Submit][Status][Discuss] Description 你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理，公司使用更加先进的机械设备生
2300: [HAOI2011]防线修建 set维护凸包
ws_fqk
02-23 490
从前写货币兑换的时候，写的splay维护凸包，没调出来。 后来写防线修建，我用了set维护凸包，A啦。 从前写永无乡的时候，写的splay启发式合并，没调出来。 后来写dispatching，我用了trep启发式合并，A啦。 我再也不相信splay了T T。以上纯属一个蒟蒻的扯淡。#include<iostream> #include<cstdio> #include<set> #inclu
bzoj-1492 货币兑换Cash (1)——平衡树维护凸包
140142
07-11 1985
题意： 有n天和m的初始金钱，用来购买AB两种纪念券； n天里每天都有AB的价格，每天可以进行这样的操作； 1.卖出手中x%的纪念券(AB分别都卖出x%)； 2.用x的金钱买入纪念券，买入AB券的比例在第i天为Rate i； 求n天过去之后所获得的最大收益； 金钱和券数均为实数； n 题解： 首先，虽然题中的买入和卖出都是任意数量的，但是同样的纪念券，分几天卖出得到的收
「NOI 2007」货币兑换
forever_dreams的博客
10-17 195
算法标签：斜率优化dp+cdq维护凸包
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值