码迷,mamicode.com
首页 > 其他好文 > 详细

cdq分治入门--BZOJ1492: [NOI2007]货币兑换Cash

时间:2017-11-23 08:15:37      阅读:128      评论:0      收藏:0      [点我收藏+]

标签:name   return   ==   one   void   std   amp   cdq分治   转化   

n<=100000天,一开始有s块钱,每天股票A价格ai,B价格bi,每天可以做的事情:卖出股票;按A:B=RTi的比例买入股票。问最后的最大收益。股票可以为浮点数,答案保留三位。

脚指头想想就知道是:某一天全部买进来,某一天全部卖出去,没有说买一半卖一半的。

那就可以dp了,f(i)表示前i天最大收益,技术分享图片,其中Xi表示用f(i)块钱在第i天能买多少A券,Yi表示f(i)块前第i天买多少B券,可以自己算,n方过不了。

现要找max(Ai*Xj+Bi*Yj),考虑两个状态j,k,j比k优时技术分享图片

整理得技术分享图片,前提是Xj-Xk

那就平衡树维护一下一个递增的(Yj-Yk)/(Xj-Xk),即维护一个凸包即可,难写,略。

cdq就是这样把一个在线的东西强行转化成离线。

solve(l,r)表示把这个区间里的f算完,solve(l,mid)之后,用(l,mid)的状态更新(mid+1,r)的状态,然后solve(mid+1,r),这就是一个分治。

为了使这个更新过程顺利完成,在solve(l,mid)时需要找到这个凸包,可以通过维护Xi的单调,然后直接一个栈保存单调的斜率即可;mid+1到r这一段,也需要保证-Ai/Bi的单调,这个可以预处理出来。现在就是一个离线问题,一边凸包单增被离线实现了,一边-Ai/Bi预处理排序好了,那就可以两个指针直接扫一遍更新了。

trick!!(Xj-Xk)可能等于0。。。。。。

技术分享图片
 1 #include<stdio.h>
 2 #include<string.h>
 3 #include<stdlib.h>
 4 #include<algorithm>
 5 #include<math.h>
 6 //#include<iostream>
 7 using namespace std;
 8 
 9 int n,s;
10 #define maxn 200011
11 double a[maxn],b[maxn],rt[maxn];
12 
13 int num[20][maxn];
14 void mergesort(int L,int R,int cur)
15 {
16     if (L==R)
17     {
18         num[cur][L]=L;
19         return;
20     }
21     const int mid=(L+R)>>1;
22     mergesort(L,mid,cur+1),mergesort(mid+1,R,cur+1);
23     int i=L,j=mid+1,k=i;
24     while (i<=mid && j<=R)
25     {
26         if (-a[num[cur+1][i]]/b[num[cur+1][i]]>-a[num[cur+1][j]]/b[num[cur+1][j]])
27             num[cur][k++]=num[cur+1][i++];
28         else num[cur][k++]=num[cur+1][j++];
29     }
30     while (i<=mid) num[cur][k++]=num[cur+1][i++];
31     while (j<=R) num[cur][k++]=num[cur+1][j++];
32 }
33 
34 double f[maxn],xx[maxn],yy[maxn];int sta[maxn],top,numf[20][maxn];
35 double calc(int i,int j) {return fabs(xx[i]-xx[j])>1e-9?(yy[i]-yy[j])/(xx[i]-xx[j]):1e300;}
36 //i>j
37 void solve(int L,int R,int cur)
38 {
39     if (L==R)
40     {
41         numf[cur][L]=L;
42         f[L]=max(f[L],f[L-1]);
43         yy[L]=f[L]/(rt[L]*a[L]+b[L]);
44         xx[L]=f[L]*rt[L]/(rt[L]*a[L]+b[L]);
45         return;
46     }
47     const int mid=(L+R)>>1;
48     solve(L,mid,cur+1);
49     top=0;
50     for (int i=L;i<=mid;i++)
51     {
52         const int id=numf[cur+1][i];
53         while (top>1 && calc(id,sta[top])-calc(sta[top],sta[top-1])>-1e-9) top--;
54         sta[++top]=id;
55     }
56     for (int i=mid+1,j=1;i<=R;i++)
57     {
58         const int &id=num[cur+1][i];
59         while (j<top && -a[id]/b[id]-calc(sta[j+1],sta[j])<1e-9) j++;
60         f[id]=max(f[id],a[id]*xx[sta[j]]+b[id]*yy[sta[j]]);
61     }
62     solve(mid+1,R,cur+1);
63     int i=L,j=mid+1,k=i;
64     while (i<=mid && j<=R)
65     {
66         if (xx[numf[cur+1][i]]<xx[numf[cur+1][j]]) numf[cur][k++]=numf[cur+1][i++];
67         else numf[cur][k++]=numf[cur+1][j++];
68     }
69     while (i<=mid) numf[cur][k++]=numf[cur+1][i++];
70     while (j<=R) numf[cur][k++]=numf[cur+1][j++];
71 }
72 
73 int main()
74 {
75     scanf("%d%d",&n,&s);
76     for (int i=1;i<=n;i++) scanf("%lf%lf%lf",&a[i],&b[i],&rt[i]);
77     mergesort(1,n,0);
78 //    for (int i=0;i<=2;i++){
79 //        for (int j=1;j<=n;j++)
80 //            cout<<(-a[num[i][j]]/b[num[i][j]])<<‘ ‘;cout<<endl;}
81     f[0]=s;for (int i=1;i<=n;i++) f[i]=0;
82     solve(1,n,0);
83 //    for (int i=1;i<=n;i++) cout<<f[i]<<‘ ‘;cout<<endl;
84     printf("%.3f\n",f[n]);
85     return 0;
86 }
View Code

 

cdq分治入门--BZOJ1492: [NOI2007]货币兑换Cash

标签:name   return   ==   one   void   std   amp   cdq分治   转化   

原文地址:http://www.cnblogs.com/Blue233333/p/7881365.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!