[关闭]
@ysner 2018-10-05T01:05:58.000000Z 字数 1638 阅读 2094

[CF739E]Gosha is hunting

DP 二分 贪心


题面

现在一共有只神奇宝贝。 你有个宝贝球和个超级球。
宝贝球抓到第只神奇宝贝的概率是 ​ ,超级球抓到的概率则是
不能往同一只神奇宝贝上使用超过一个同种的球,但是可以往同一只上既使用宝贝球又使用超级球(都抓到算一个)。
请合理分配每个球抓谁,使得你抓到神奇宝贝的总个数期望最大,并输出这个值。

解析

这题讲课时有一种贪心做法。
先按超级球捕捉概率排序。
二分最后一次用超级球的位置,则前面要不,要不;后面要不,要不
这样就可以拿个堆贪心地选取以最大化期望。

看到这题,裸是三维的。
表示到第个宝贝,用个宝贝球,个超级球的期望。
然而没分。

以下摘自巨佬yyb的博客
发现可以凸优化,对于其中一个球给它二分一个权值,表示每使用一次就需要额外花费掉这么多的权值,同时不再限制使用的个数。
然后忽略这一个限制,做,利用最优解使用的这种球的个数以及限制个数继续二分。
两维都可以这么做,复杂度

怎么说呢,通过二分强加权值来代替限制这个操作好像出现不止一次了。例如[国家集训队2]Tree I
知识学不学得会是一回事,会不会灵活运用又是一回事。。。

还有个细节,如果要嵌套二分,精度要求要翻倍。

  1. #include<iostream>
  2. #include<cmath>
  3. #include<cstdio>
  4. #include<cstdlib>
  5. #include<cstring>
  6. #include<algorithm>
  7. #include<set>
  8. #define ll long long
  9. #define re register
  10. #define il inline
  11. #define db double
  12. #define eps 1e-8
  13. #define fp(i,a,b) for(re int i=a;i<=b;++i)
  14. #define fq(i,a,b) for(re int i=a;i>=b;--i)
  15. using namespace std;
  16. const int N=5e3+100;
  17. int n,ta,tb;
  18. db f[N],fa[N],fb[N];
  19. struct dat{db a,b;}t[N];
  20. il void check(re db w1,re db w2)
  21. {
  22. fp(i,1,n)
  23. {
  24. f[i]=f[i-1];fa[i]=fa[i-1];fb[i]=fb[i-1];
  25. if(f[i-1]+t[i].a-w1>f[i]) f[i]=f[i-1]+t[i].a-w1,fa[i]=fa[i-1]+1,fb[i]=fb[i-1];
  26. if(f[i-1]+t[i].b-w2>f[i]) f[i]=f[i-1]+t[i].b-w2,fa[i]=fa[i-1],fb[i]=fb[i-1]+1;
  27. if(f[i-1]+t[i].a+t[i].b-t[i].a*t[i].b-w1-w2>f[i]) f[i]=f[i-1]+t[i].a+t[i].b-t[i].a*t[i].b-w1-w2,fa[i]=fa[i-1]+1,fb[i]=fb[i-1]+1;
  28. }
  29. }
  30. int main()
  31. {
  32. ios::sync_with_stdio(false);
  33. cin>>n>>ta>>tb;
  34. fp(i,1,n) cin>>t[i].a;
  35. fp(i,1,n) cin>>t[i].b;
  36. re db l1=0,r1=1,mid1,l2,r2,mid2;
  37. while(l1+eps<r1)
  38. {
  39. mid1=(l1+r1)/2;
  40. l2=0,r2=1;
  41. while(l2+eps<r2)
  42. {
  43. mid2=(l2+r2)/2;
  44. check(mid1,mid2);
  45. if(fb[n]>tb) l2=mid2;else r2=mid2;
  46. }
  47. check(mid1,r2);
  48. if(fa[n]>ta) l1=mid1;else r1=mid1;
  49. }
  50. check(r1,r2);
  51. printf("%.4lf\n",f[n]+r1*ta+r2*tb);
  52. return 0;
  53. }
添加新批注
在作者公开此批注前,只有你和作者可见。
回复批注