斜率优化(CDQ分治，Splay平衡树)：BZOJ 1492: [NOI2007]货币兑换Cash

Description

Input

第一行两个正整数N、S，分别表示小Y 能预知的天数以及初始时拥有的钱数。接下来N 行，第K 行三个实数AK、BK、RateK，意义如题目中所述

Output

只有一个实数MaxProfit，表示第N 天的操作结束时能够获得的最大的金钱数目。答案保留3 位小数。

Sample Input

3 100
1 1 1
1 2 2
2 2 3

Sample Output

225.000

HINT

测试数据设计使得精度误差不会超过10-7。
对于40%的测试数据，满足N ≤ 10；
对于60%的测试数据，满足N ≤ 1 000；
对于100%的测试数据，满足N ≤ 100 000；

 1 #include <iostream>
 2 #include <cstring>
 3 #include <cstdio>
 4 #include <algorithm>
 5 using namespace std;
 6 const double eps=1e-8;
 7 const int maxn=100010;
 8 struct Node{
 9     double a,b,r,x,y,k;
10     int id;
11 }d[maxn],t[maxn];
12 double f[maxn];
13 int st[maxn],cnt,n;
14 double Abs(double a){return a>0?a:-a;}
15 double K(int a,int b){
16     if(Abs(d[a].x-d[b].x)<eps)return 1e20;
17     return (d[a].y-d[b].y)/(d[a].x-d[b].x);
18 }
19 void Solve(int l,int r){
20     if(l==r){
21         f[l]=max(f[l],f[l-1]);
22         d[l].y=f[l]/(d[l].r*d[l].a+d[l].b);
23         d[l].x=d[l].y*d[l].r;
24         return;
25     }
26     int mid=(l+r)>>1,t1=l,t2=mid+1;
27     for(int i=l;i<=r;i++){
28         if(d[i].id<=mid)
29             t[t1++]=d[i];
30         else
31             t[t2++]=d[i];
32     }
33     for(int i=l;i<=r;i++)d[i]=t[i];
34     Solve(l,mid);
35     cnt=0;
36     for(int i=l;i<=mid;i++){
37         while(cnt>1&&K(i,st[cnt])-K(st[cnt],st[cnt-1])>=eps)cnt--;
38         st[++cnt]=i;
39     }
40     int fir=1;
41     for(int i=mid+1;i<=r;i++){
42         while(fir<cnt&&K(st[fir+1],st[fir])-d[i].k>=eps)fir++;
43         f[d[i].id]=max(f[d[i].id],d[st[fir]].x*d[i].a+d[st[fir]].y*d[i].b);
44     }
45     Solve(mid+1,r);
46     t1=l;t2=mid+1;
47     for(int i=l;i<=r;i++){
48         if(t2==r+1||(d[t2].x-d[t1].x>=eps||Abs(d[t2].x-d[t1].x)<eps&&d[t2].y-d[t1].y>=eps)&&t1<=mid)
49             t[i]=d[t1++];
50         else
51             t[i]=d[t2++];
52     }
53     for(int i=l;i<=r;i++)
54         d[i]=t[i];
55 }
56
57
58 bool cmp(Node a,Node b){
59     return a.k>b.k;
60 }
61 int main(){
62     scanf("%d%lf",&n,&f[1]);
63     for(int i=1;i<=n;i++){
64         scanf("%lf%lf%lf",&d[i].a,&d[i].b,&d[i].r);
65         d[i].k=-d[i].a/d[i].b;
66         d[i].id=i;
67     }
68     sort(d+1,d+n+1,cmp);
69     Solve(1,n);
70     printf("%.3lf\n",f[n]);
71     return 0;
72 }

　　这里还有Splay。

  1 #include <iostream>
  2 #include <cstring>
  3 #include <cstdio>
  4 using namespace std;
  5 const double eps=1e-8;
  6 const double INF=1e20;
  7 const int maxn=100010;
  8 int ch[maxn][2],fa[maxn],rt,tot,n;
  9 double X[maxn],Y[maxn],lk[maxn],rk[maxn],ans;
 10 double fabs(double x){return (x>0)?x:-x;}
 11 void Rotate(int x){
 12     int y=fa[x],g=fa[y],c=ch[y][1]==x;
 13     ch[y][c]=ch[x][c^1];fa[ch[y][c]]=y;
 14     ch[x][c^1]=y;fa[y]=x;fa[x]=g;
 15     if(g)ch[g][ch[g][1]==y]=x;
 16 }
 17
 18 void Splay(int x,int g=0){
 19     for(int y;(y=fa[x])!=g;Rotate(x))
 20         if(fa[y]!=g)Rotate((ch[fa[y]][1]==y)==(ch[y][1]==x)?y:x);
 21     if(!g)rt=x;
 22 }
 23
 24 double Get_K(int j,int k){
 25     if(fabs(X[j]-X[k])<=eps)return INF;
 26     else return (Y[j]-Y[k])/(X[j]-X[k]);
 27 }
 28
 29 int Get_Prev(){
 30     int p=ch[rt][0],ret=p;
 31     while(p){
 32         if(Get_K(rt,p)+eps>=lk[p])p=ch[p][0];
 33         else ret=p,p=ch[p][1];
 34     }
 35     return ret;
 36 }
 37
 38 int Get_Succ(){
 39     int p=ch[rt][1],ret=p;
 40     while(p){
 41         if(Get_K(p,rt)<=rk[p]+eps)p=ch[p][1];
 42         else ret=p,p=ch[p][0];
 43     }
 44     return ret;
 45 }
 46
 47 void Insert(int &r,int pre,int p){
 48     if(r==0){r=p;fa[p]=pre;return;}
 49     if(X[p]<=X[r]+eps)Insert(ch[r][0],r,p);
 50     else Insert(ch[r][1],r,p);
 51 }
 52
 53 void Update(int p){
 54     Splay(p);
 55     if (ch[p][0]){
 56         int l=Get_Prev();
 57         Splay(l,p);ch[l][1]=0;
 58         lk[p]=rk[l]=Get_K(p,l);
 59     }
 60     else lk[p]=INF;
 61     if (ch[p][1]){
 62         int r=Get_Succ();
 63         Splay(r,p); ch[r][0]=0;
 64         rk[p]=lk[r]=Get_K(r,p);
 65     }
 66     else rk[p]=-INF;
 67     if (lk[p]<=rk[p]+eps){
 68         rt=ch[p][0]; ch[rt][1]=ch[p][1]; fa[ch[p][1]]=rt; fa[rt]=0;
 69         rk[rt]=lk[ch[p][1]]=Get_K(ch[rt][1],rt);
 70     }
 71 }
 72
 73 int Get_Pos(double k){
 74     int p=rt;
 75     while(p){
 76         if(lk[p]+eps>=k&&k+eps>=rk[p])break;
 77         if(lk[p]<k+eps)p=ch[p][0];
 78         else p=ch[p][1];
 79     }
 80     return p;
 81 }
 82
 83 double Get_Ans(double a,double b){
 84     int p=Get_Pos(-b/a);
 85     return a*Y[p]+b*X[p];
 86 }
 87
 88 int main(){
 89 #ifndef ONLINE_JUDGE
 90     freopen("cash.in","r",stdin);
 91     freopen("cash.out","w",stdout);
 92 #endif
 93     double a,b,rate;
 94     scanf("%d%lf",&n,&ans);
 95     for(int i=1;i<=n;i++){
 96         scanf("%lf%lf%lf",&a,&b,&rate);
 97         if(i!=1)ans=max(ans,Get_Ans(a,b));
 98         X[i]=ans/(rate*a+b);
 99         Y[i]=X[i]*rate;
100         Insert(rt,0,i);
101         Update(i);
102     }
103     printf("%.3f\n",ans);
104     return 0;
105 }

转载于:https://www.cnblogs.com/TenderRun/p/5307998.html

斜率优化(CDQ分治，Splay平衡树)：BZOJ 1492: [NOI2007]货币兑换Cash相关推荐

BZOJ 1492: [NOI2007]货币兑换Cash [CDQ分治斜率优化DP]
传送门题意:不想写... 扔链接就跑好吧我回来了首先发现每次兑换一定是全部兑换,因为你兑换说明有利可图,是为了后面的某一天两种卷的汇率差别明显而兑换那么一定拿全利啊,一定比多天的组合好 $f[ ...
bzoj 1492: [NOI2007]货币兑换Cash
Description 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和 B纪念券(以下简称B券).每个持有金券的顾客都有一个自己的帐户.金券的数目可以是一个 ...
【BZOJ1492】[NOI2007]货币兑换Cash 斜率优化+cdq分治
[BZOJ10492][NOI2007]货币兑换Cash Description 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和 B纪念券(以下简称B券).每 ...
bzoj 2149 拆迁队斜率优化+cdq分治
题面题目传送门解法从来没写过这样的-- 第一问非常简单,能够从 j j j转移到 i i i的条件显然为 a [ i ] − a [ j ] ≥ i − j a[i]-a[j]≥i-j a[i] ...
【bzoj 1492】【codevs 1797】 [NOI2007]货币兑换Cash (dp+cdq分治)
1492: [NOI2007]货币兑换Cash Time Limit: 5 Sec Memory Limit: 64 MB Submit: 3803 Solved: 1604 [Submit][S ...
[BZOJ1492] [NOI2007]货币兑换Cash 斜率优化+cdq/平衡树维护凸包
1492: [NOI2007]货币兑换Cash Time Limit: 5 Sec Memory Limit: 64 MB Submit: 5907 Solved: 2377 [Submit][S ...
bzoj千题计划237：bzoj1492: [NOI2007]货币兑换Cash
http://www.lydsy.com/JudgeOnline/problem.php?id=1492 dp[i] 表示第i天卖完的最大收益朴素的dp: 枚举从哪一天买来的在第i天卖掉,或者是不 ...
bzoj1492 [NOI2007]货币兑换Cash (斜率DP+cdq分治)
题意:到处都找得到. 我没看错的话当年考试的时候的题面里头,是提示了买卖一定是全部买入和卖出的.这样一来就好办了.cdq的论文里面那个F并不是她所说的那样,而是就是那个最优值.方程转移的时候实际上是枚 ...
YbtOJ#493-最大分数【斜率优化dp,分治】
正题题目链接:http://172.17.55.160/contest/117/problem/1 题目大意 nnn个数的一个序列,给其中的一些数打上标记. 一个标记方案的贡献为s1s_1s1表示 ...