Codeforces.739E.Gosha is hunting(DP 带权二分)
\(Description\)
有\(n\)只精灵,两种精灵球(高级和低级),每种球能捕捉到第\(i\)只精灵的概率已知。求用\(A\)个低级球和\(B\)个高级球能捕捉到精灵数的最大期望。
\(n\leq10^5\)。
\(Solution\)
设\(f[i][a][b]\)表示前\(i\)只用了\(a\)个低级球,\(b\)个高级球的最大期望。转移时四种情况显然。复杂度\(\mathcal O(nAB)\)。
随着某种球可使用数的增多,f应是凸函数,即增长越来越慢。而且两种球都满足这个性质。
于是可以wqs二分套wqs二分了。。
没有个数限制的话,取个\(\max\),记一下个数就可以了。复杂度\(\mathcal O(nlog^2n)\)。
误差\(\leq 10^{-4}\),因为最后要\(*A/B\),所以\(eps\)应是\(10^{-8}\)...?
最后必须取\(r\)感觉不太懂。。
[Update] 19.2.11
二分的时候只要保证恰好取到\(k\)个就可以了,斜率具体是多少无所谓。而本题\(r\)才是正确的位置...
应该是这么理解吧...
总结:对于有着次数/段数之类的限制,可以使用带权二分来消掉这一限制,从而可以进行简单的快速DP。
//46ms 0KB
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
#define eps (1e-12)
const int N=2003;
int n,A,B,na[N],nb[N];
double pa[N],pb[N],Ans;
void Solve(double ca,double cb)
{
na[0]=nb[0]=0;
double las=0, now;
for(int i=1; i<=n; ++i, las=now)
{
now=las, na[i]=na[i-1], nb[i]=nb[i-1];
if(las+pa[i]-ca>now) now=las+pa[i]-ca, na[i]=na[i-1]+1;
if(las+pb[i]-cb>now) now=las+pb[i]-cb, nb[i]=nb[i-1]+1, na[i]=na[i-1];
if(las+pa[i]+pb[i]-pa[i]*pb[i]-ca-cb>now)//1-(1-pa)(1-pb)
now=las+pa[i]+pb[i]-pa[i]*pb[i]-ca-cb, na[i]=na[i-1]+1, nb[i]=nb[i-1]+1;
}
Ans=now;
}
int main()
{
scanf("%d%d%d",&n,&A,&B);
for(int i=1; i<=n; ++i) scanf("%lf",&pa[i]);
for(int i=1; i<=n; ++i) scanf("%lf",&pb[i]);
double l1=0,r1=1,mid1,l2,r2,mid2;//每个球0/1的权值就可以了啊
while(r1>=l1+eps)
{
mid1=(l1+r1)*0.5;
l2=0, r2=1;
while(r2>=l2+eps)
{
if(Solve(mid1,mid2=(l2+r2)*0.5),nb[n]>B) l2=mid2;
else r2=mid2;
}
if(Solve(mid1,r2),na[n]>A) l1=mid1;//最优可行的是r2?反正不是l2。。
else r1=mid1;
}
Solve(r1,r2);//最后Check一遍r。。
printf("%.5lf",Ans+A*r1+B*r2);
return 0;
}
Codeforces.739E.Gosha is hunting(DP 带权二分)的更多相关文章
- 洛谷 4383 [八省联考2018]林克卡特树lct——树形DP+带权二分
题目:https://www.luogu.org/problemnew/show/P4383 关于带权二分:https://www.cnblogs.com/flashhu/p/9480669.html ...
- 洛谷.4383.[八省联考2018]林克卡特树lct(树形DP 带权二分)
题目链接 \(Description\) 给定一棵边带权的树.求删掉K条边.再连上K条权为0的边后,新树的最大直径. \(n,K\leq3\times10^5\). \(Solution\) 题目可以 ...
- 6.13校内互测 (DP 带权二分 斜率优化)
丘中有麻plant 改自这儿,by ZBQ. 还有隐藏的一页不放了.. 直接走下去的话,如果开始时间确定那么到每个点的时间确定,把time减去dis就可以去掉路程的影响了. 这样对于减去d后的t,如果 ...
- codeforces 739E - Gosha is hunting
这道题有三种做法,感受一下: 感觉到了歪果仁费尽脑汁想出来的神仙贪心脑洞题被中国人套路算法踩爆的凄凉...(我的名字是不是暴露了我的真实实力) ============================ ...
- DP的各种优化(动态规划,决策单调性,斜率优化,带权二分,单调栈,单调队列)
前缀和优化 当DP过程中需要反复从一个求和式转移的话,可以先把它预处理一下.运算一般都要满足可减性. 比较naive就不展开了. 题目 [Todo]洛谷P2513 [HAOI2009]逆序对数列 [D ...
- CF 739E Gosha is Hunting
有 $n$ 个 Pokemon,你有 $A$ 个一类精灵球,$B$ 个二类精灵球 分别给出每个 Pokemon 被这两类精灵球捕捉的概率 求抓到 Pokemon 的最优期望个数 $n\leq 2000 ...
- P4383 [八省联考2018]林克卡特树lct 树形DP+凸优化/带权二分
$ \color{#0066ff}{ 题目描述 }$ 小L 最近沉迷于塞尔达传说:荒野之息(The Legend of Zelda: Breath of The Wild)无法自拔,他尤其喜欢游戏中的 ...
- 洛谷P2619 [国家集训队2]Tree I(带权二分,Kruscal,归并排序)
洛谷题目传送门 给一个比较有逼格的名词--WQS二分/带权二分/DP凸优化(当然这题不是DP). 用来解决一种特定类型的问题: 有\(n\)个物品,选择每一个都会有相应的权值,需要求出强制选\(nee ...
- BZOJ_4609_[Wf2016]Branch Assignment_决策单调性+带权二分
BZOJ_4609_[Wf2016]Branch Assignment_决策单调性+带权二分 Description 要完成一个由s个子项目组成的项目,给b(b>=s)个部门分配,从而把b个部门 ...
随机推荐
- 树形dp&&树的重心(D - Godfather POJ - 3107)
题目链接:https://cn.vjudge.net/contest/277955#problem/D 题目大意:求树的重心(树的重心指的是树上的某一个点,删掉之后形成的多棵树中节点数最大值最小). ...
- PHP返回Json数据函数封装
/** * 返回Json数据 * @param int $code * @param string $message * @param array $data * @return string */ ...
- Unity 3(二):Unity在AOP方面的应用
本文关注以下方面(环境为VS2012..Net Framework 4.5以及Unity 3): AOP简介: Interception using Unity示例 配置文件示例 一.AOP简介 AO ...
- Strusts2笔记6--拦截器
拦截器: Struts2的大多数核心功能都是通过拦截器实现的.拦截器之所以称之为“拦截器”,是因为它可以在执行Action之前或之后拦截下用户请求,执行一些操作,以增强Action方法的功能. Str ...
- web项目启动首页能访问接口报404
这个问题如果控制台没有报错,然后看一下日志,看看spring容器是否加载: 如果是一直卡在这里了,多半是resources文件下的配置有问题,或者是resources目录不是源文件,工具无法识别 我遇 ...
- 【HASPDOG】hasp_update参数f和i区别
[root@BICServer-TX shared]# ./hasp_update This is a simple demo program for the Sentinel Update and ...
- Flask:初次使用Flask-SQLAlchemy读取SQLite3
Windows 10家庭中文版,Python 3.6.4,Flask 1.0.2,Eclipse Oxygen.1a Release (4.7.1a),PyDev 6.3.2 SQLAlchemy是一 ...
- 【前端开发】localStorage的用法
localStorage.setItem("name","value") //存储name的值 var type = localStorage.getItem ...
- Python爬虫学习1: Requests模块的使用
Requests函数库是学习Python爬虫必备之一, 能够帮助我们方便地爬取. Requests: 让HTTP服务人类. 本文主要参考了其官方文档. Requests具有完备的中英文文档, 能完全满 ...
- No.6 selenium学习之路之下拉框Select
HTML中,标签显示为select,有option下拉属性的为Select弹框 1.Xpath定位 Xpath语法,顺序是从1开始,编程语言中是0开始