stone2 [期望]
\(\mathcal{Description}\)
有 \(n\) 堆石子,依次编号为 \(1, 2,\ldots , n\),其中第 \(i\) 堆有 \(a_i\) 颗石子
你每次在仍然有石子的石子堆中等概率随机选择一堆石子,并取走其中一颗石子
求第 \(1\) 堆石子被取走的时间的期望
\(n\leq 5\times 10^5,a_i\leq 5\times 10^5\)
\(\mathcal{Solution}\)
这题其实也不难,然而也不是考虑\(DP\),和stone一样
问题要求的实际就是在第\(1\)堆石子被取完之前,总共有多少个石子被拿走了
显然\(a_1\)被拿完了,再考虑期望的线性性,由于你取走其他堆石子对当前堆没有影响,可以单独考虑每一堆石子被取了多少个
先考虑\(a_i\)没有全部被取完
假设第\(i\)堆石子被取走了\(x\left(0\leq x\leq a_i-1\right)\)颗石子,因为第\(i\)堆石子没被取完,而没有继续被取走肯定是因为第\(1\)堆石子被取完了
设\(1\)为从第一堆石子中取一颗石子,\(0\)为从第二堆石子中取一颗石子
那么方案数就是有\(a_1\)个\(1\)和\(x\)个\(0\)的长度为\(a_1+x\)的且最后一位是\(1\)的二进制串的个数
则有\(\begin{pmatrix}a_1+x-1 \\ x\end{pmatrix}\)种方案数,总方案数为\(2^{a_1+x}\),那么概率就是\(p_x=\frac{\begin{pmatrix}a_1+x-1\\ x\end{pmatrix}}{2^{a_1+x}}\)
若\(x=a_i\),这个看起来没原来那么好算了,想到这两种情况的概率加起来应该等于\(1\),所以这种情况的概率就是\(q=1-\sum\limits_{j=0}^{a_i-1}p_j\)
则我们得到\(E_i=\sum\limits_{j=0}^{a_i-1}j*p_j+qa_i\)
总期望就是\(E=\sum\limits_{i=2}^nE_i\),写复杂点就是
\(E=\left(\sum\limits_{i=2}^n\left(\sum\limits_{j=0}^{a_i-1}j*\frac{\begin{pmatrix}a_1+j-1 \\ j\end{pmatrix}}{2^{a_1+j}}\right)+a_i\left(1-\sum\limits_{j=0}^{a_i-1}\frac{\begin{pmatrix}a_1+j-1 \\ j\end{pmatrix}}{2^{a_1+j}}\right)\right)+a_1\)
这个东西怎么维护呢,直接考虑\(a_i\)变成\(a_i+1\)的情况,我们考虑里面那个\(\sum\)的变化,实际上只有枚举上界增大\(1\)变成了\(a_i\)我们直接对\(a_i\)的所有取值的答案都预处理出来就行了
\(\mathcal{Code}\)
/*******************************
Author:Morning_Glory
LANG:C++
Created Time:2019年11月08日 星期五 16时20分30秒
*******************************/
#include <cstdio>
#include <fstream>
using namespace std;
const int maxn = 1000006;
const int lim = 1000000;
const int h = 500000;
const int mod = 323232323;
//cin为我打的快读板子,详细内容可去看以前的代码,总是在这里写感觉有点影响阅读
int n,p,ans,a;
int fac[maxn],ifac[maxn],inv[maxn],mi[maxn],f[maxn],g[maxn];
int C (int n,int m){ return 1ll*fac[n]*ifac[n-m]%mod*ifac[m]%mod;}
int main()
{
fac[0]=ifac[0]=mi[0]=inv[1]=1;
for (int i=2;i<=lim;++i) inv[i]=(mod-1ll*mod/i*inv[mod%i]%mod);
for (int i=1;i<=lim;++i) fac[i]=1ll*fac[i-1]*i%mod,ifac[i]=1ll*ifac[i-1]*inv[i]%mod,mi[i]=1ll*mi[i-1]*inv[2]%mod;
cin>>n>>a;
ans=a;
//枚举上界为i
g[0]=mi[a];
for (int i=1;i<=h;++i){
int p=1ll*C(a+i-1,i)*mi[a+i]%mod;
f[i]=(f[i-1]+1ll*i*p%mod)%mod;
g[i]=(g[i-1]+p)%mod;
}
//石子数为i的答案
for (int i=h;i>=1;--i) f[i]=(f[i-1]+1ll*i*(mod+1-g[i-1])%mod)%mod;
for (int i=2;i<=n;++i) cin>>a,ans=(ans+f[a])%mod;
printf("%d\n",ans);
return 0;
}
如有哪里讲得不是很明白或是有错误,欢迎指正
如您喜欢的话不妨点个赞收藏一下吧
stone2 [期望]的更多相关文章
- 【BZOJ-3143】游走 高斯消元 + 概率期望
3143: [Hnoi2013]游走 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2264 Solved: 987[Submit][Status] ...
- bzoj1415[NOI2005]聪聪和可可-期望的线性性
这道题之前我写过一个巨逗比的写法(传送门:http://www.cnblogs.com/liu-runda/p/6220381.html) 当时的原因是这道题可以抽象出和"绿豆蛙的归宿&qu ...
- hdu 4481 Time travel(高斯求期望)(转)
(转)http://blog.csdn.net/u013081425/article/details/39240021 http://acm.hdu.edu.cn/showproblem.php?pi ...
- 【BZOJ3036】绿豆蛙的归宿 概率与期望
最水的概率期望,推荐算法合集之<浅析竞赛中一类数学期望问题的解决方法> #include <iostream> #include <cstdio> using na ...
- UVA&&POJ离散概率与数学期望入门练习[4]
POJ3869 Headshot 题意:给出左轮手枪的子弹序列,打了一枪没子弹,要使下一枪也没子弹概率最大应该rotate还是shoot 条件概率,|00|/(|00|+|01|)和|0|/n谁大的问 ...
- 【BZOJ-1426】收集邮票 概率与期望DP
1426: 收集邮票 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 261 Solved: 209[Submit][Status][Discuss] ...
- 【BZOJ-1419】Red is good 概率期望DP
1419: Red is good Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 660 Solved: 257[Submit][Status][Di ...
- [NOIP2016]换教室 D1 T3 Floyed+期望DP
[NOIP2016]换教室 D1 T3 Description 对于刚上大学的牛牛来说, 他面临的第一个问题是如何根据实际情况中情合适的课程. 在可以选择的课程中,有2n节课程安排在n个时间段上.在第 ...
- 关于在left join的on子句中限制左边表的取值时出现非期望的结果
使用的SQL大概是这样的: select * from A left join B on A.id=B.id and A.id>10; --错误的使用 我们期望的结果集应该是 A中的id> ...
随机推荐
- js中实现函数防抖跟函数节流
最近刚接触两个新概念函数防抖与函数节流,虽然这些内容网上可以搜到很多,大家都有自己的一套的理解方式,都写得很好, 而自己则想在理解的基础上自己把代码实现一遍,加深印象. 一.函数防抖 假如我们有这样的 ...
- nginx 反向代理之 proxy_cache
proxy_cache将从C上获取到的数据根据预设规则存放到B上(内存+磁盘)留着备用,A请求B时,B会把缓存的这些数据直接给A,而不需要再去向C去获取. proxy_cache相关功能生效的前提是, ...
- Exploiting ConvNet Diversity for Flooding Identification
语义分割洪水区域. 空洞卷积和反卷积组合,结果再用svm学习如何组合,能获得更好的效果. 直接对不同网络的结果进行投票会得到更差的结果. 消融研究(Ablation Study):类似控制变量法,就对 ...
- 使用adb连接Mumu模拟器
1)下载Mumu模拟器 2)运行Mumu模拟器 3)找到mumu安装目录下的MuMu\emulator\nemu\vmonitor\bin目录 4)在当前目录打开cmd,执行 adb connect ...
- Nginx服务器作反向代理实现内部局域网的url转发配置
情景 由于公司内网有多台服务器的http服务要映射到公司外网静态IP,如果用路由的端口映射来做,就只能一台内网服务器的80端口映射到外网80端口,其他服务器的80端口只能映射到外网的非80端口.非80 ...
- <每日 1 OJ> -内存文件系统
蛮有意思的,主要考查链表和树结构的知识. 代码如下: ************************************************************************* ...
- 利用python实现汉字转拼音
安装:pip install pypinyin import pypinyin # 不带声调的(style=pypinyin.NORMAL) def pinyin(word): s = '' for ...
- virtualenv, conda, pip分别是什么
自己一直使用virtualenv,但是发现很多工具或框架都是以来conda,于是就网上搜了下二者的区别,感觉这篇文章讲的比较清楚:https://blog.csdn.net/zhouchen1998/ ...
- vue---props进行双向数据绑定报错
在使用vue进行组件开发的时候,遇到一个问题,父组件传递到子组件里面的值,如果在子组件里面进行改变 传递过来的"值",会报错: [Vue warn]: Avoid mutating ...
- ThinkPhp5 数据迁移(think-migration)
ThinkPhp5 数据迁移(think-migration) 在thinkphp5中提供了数据迁移工具(think-migration),它是机遇phinx开发(文档地址:http://docs ...