1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 | #include <iostream>#include <cstdio>#include <algorithm>#include <cstring>#include <cmath>using namespace std;const int N=100010;double a[N],b[N],ra[N],v[N];struct Data { int t; double k,b,x;} g[N],p[N],q[N];inline bool cmp(const Data &x,const Data &y) { return x.x<y.x;}inline void solve(int l,int r) { if(l==r) return; int mid=(l+r)>>1,p1=l,p2=mid+1,st=1,ed=0; double Mx=0; for(int i=l;i<=r;i++) { if(g[i].t<=mid) p[p1++]=g[i]; else p[p2++]=g[i]; } for(int i=l;i<=r;i++) g[i]=p[i]; solve(l,mid); for(int i=l;i<=mid;i++) { while(st<ed&&(g[i].k-q[ed].k)*(q[ed-1].b-g[i].b)>=(g[i].k-q[ed-1].k)*(q[ed].b-g[i].b)) ed--; q[++ed]=g[i],Mx=max(Mx,v[g[i].t]); } sort(g+mid+1,g+r+1,cmp); for(int i=mid+1;i<=r;i++) { while(st<ed&&q[st].k*g[i].x+q[st].b<=q[st+1].k*g[i].x+q[st+1].b) st++; v[g[i].t]=max(v[g[i].t],max(Mx,a[g[i].t]*(q[st].k*g[i].x+q[st].b))); g[i].k=v[g[i].t]/(a[g[i].t]*ra[g[i].t]+b[g[i].t]),g[i].b=g[i].k*ra[g[i].t]; } solve(mid+1,r);p1=l,p2=mid+1; for(int i=l;i<=r;i++) { if(p2>r||(p1<=mid&&g[p1].k<=g[p2].k)) p[i]=g[p1++]; else p[i]=g[p2++]; } for(int i=l;i<=r;i++) g[i]=p[i];}int main() { int n; cin>>n>>v[0]; for(int i=1;i<=n;i++) { scanf("%lf%lf%lf",&a[i],&b[i],&ra[i]); v[i]=v[0],g[i].t=i; g[i].x=b[i]/a[i],g[i].k=v[i]/(a[i]*ra[i]+b[i]),g[i].b=g[i].k*ra[i]; } solve(1,n); printf("%.3f\n",v[n]); return 0;} |
谁能帮我解读一下这段神奇的代码？？
算法：CDQ+斜率优化
0 分
0 人评分
C语言网提供由在职研发工程师或ACM蓝桥杯竞赛优秀选手录制的视频教程，并配有习题和答疑，点击了解：
一点编程也不会写的：零基础C语言学练课程
解决困扰你多年的C语言疑难杂症特性的C语言进阶课程
从零到写出一个爬虫的Python编程课程
只会语法写不出代码？手把手带你写100个编程真题的编程百练课程
信息学奥赛或C++选手的 必学C++课程
蓝桥杯ACM、信息学奥赛的必学课程：算法竞赛课入门课程
手把手讲解近五年真题的蓝桥杯辅导课程
发表评论 取消回复