题目链接:

http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1010&cid=459

找出公式,公式有实际意义,某种情形当重复做n次实验时会出现一次,即出现的概率为1/n,现在要想出现这种情形,平均要做多少次实验,显然平均要做n次。

说一个具体的,比如掷色子,有6个点,6个点随机等概率出现。掷一次色子出现1的概率为1/6,现在想掷出1来,平均要掷色子多少次,即次数的数学期望是多少。

可以证明:

设掷i次色子才出现1的概率为p[i],则有

p[1] = 1/6;

p[2] = (5/6) *(1/6);第一次没有出现

p[3] = (5/6)^2*(1/6)

`````

p[n] = (5/6)^n-1*(1/6);前n-1次没有出现

``````

所以次数的数学期望为sum(k*p[k]),(k>=1,k为自然数,可以取到无穷)  ,据观察,有p[i] = p[i-1]*(5/6);

式子1:s = 1*p[1] + 2*p[2] + ``` +n*p[n] +````

式子2:(5/6)*s =     1*p[1] + ````(n-1)*p[n]+````(式子2为式子1的左右两边同时乘以5/6得到)

用式子1-式子2得

(1/6)*s = 1/6 + (1/6)*(5/6) + (1/6)*(5/6)^2````+(1/6)*(5/6)^n+```

又等比数列的公式得(1/6)*s =( 1/6(1 -(5/6)^n)/(1-5/6).由于n为正无穷,所有(5/6)^n =0;所以s = 6.

这个题的公式为:

情形1,出现连续n个相同后停止掷色子的次数的期望 = 1+m+m^2 + ``` + m^(n-1);

情形2,出现连续n个不同后停止掷色子的次数的期望 = 1+m/(m-1) + m^2/((m-1)*(m-2))+```+m^(n-1)/((m-1)*(m-2)*```*(m-n+1));

可以这样解释,情形1,先随便掷色子一次,后面要掷出和这个相同的期望数为m,然后仍要掷出相同的期望为m^2```

类推,后面要掷出和这个色子同色的概率为1/m,所以次数的期望为m。第三次掷色子要和前两次相同的概率为1/(m*m),所以期望为m*m```

其余就自己推吧···(未必对···)

另外一种方法:概率dp,摘自杭电解题报告

设dp[i]表示当前在 已经投掷出 i个 不相同/相同 这个状态时期望还需要投掷多少次,然后dp[i] 有如下等式:

相同:

//dp[0] = 1 + dp[1]

//dp[1] = 1 + ((m-1)dp[1] + dp[2]) / m

//dp[i] = 1 + ((m-1)dp[1] + dp[i+1]) / m

//...

//dp[n] = 0;

不相同:

//dp[0] = 1 + dp[1]

//dp[1] = 1 + (dp[1] + (m-1) dp[2]) / m

//dp[2] = 1 + (dp[1] + dp[2] + (m-2) dp[3]) / m

//dp[i] = 1 + (dp[1] + dp[2] + ... dp[i] + (m-i)dp[i+1]) / m

//...

//dp[n] = 0;

。。于是可以高斯消元。。对于第一问。。我们发现就是相当于 Typing Monkey 问题中字符串是 AAAA..AA 这一特殊情况。。解得递推式:

dp[n] = 0

dp[n-1] = dp[n] * m + 1

。。。

解开后等于等比数列求和。

(也可以直接得到这个公式。。。因为在当前状态只有 m/1 的概率可以进入下一状态,否则要重新来过。。而这一步会另总的步数 + 1。)

对于第二问。。

现在设s[i]=sigma{dp[i], 1..i},对s[i] 列方程
每个方程是关于三个相邻的s[i] 的,然后就可以线性时间解出来了。

也可以设 d[i] = dp[i] - dp[i+1].

可以得到 d[i] =  m * d[i-1]  / (m-i)

然后就是解一元一次方程...

参见:

http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=13614

http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=21631

http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=24402

http://www.cnblogs.com/ch3656468/archive/2011/05/04/2036332.html

上面的两种递推式肯定能推出上面的公式来,所以我直接用公式写的

贴代码:

 #include <cstdio>
int main()
{
// freopen("in.c","r",stdin);
int t;
while(scanf("%d",&t) != EOF)
{
while(t--)
{
double flag,m,n;
scanf("%lf%lf%lf",&flag,&m,&n);
double ans=,tmp =;
if(flag == )
{
for(int i=; i<n; ++i)
{
tmp *= m;
ans += tmp; }
}
else
{
for(int i=; i<n; ++i)
{
tmp *= (m/(m-i));
ans += tmp;
}
}
printf("%lf\n",ans);
}
}
return ;
}

dice 概率论 概率DP的更多相关文章

  1. hdu 4625 Dice(概率DP)

    Dice Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submi ...

  2. hdu 4586 Play the Dice(概率dp)

    Problem Description There is a dice with n sides, which are numbered from 1,2,...,n and have the equ ...

  3. HihoCoder - 1339 Dice Possibility(概率dp)

    题意:求用N(1<=N<=100)个骰子掷出M(1<=M<=600)的概率 分析:直接求概率可能出现6^100次方,会爆精度.可以用一个数组dp[i][j]记录用i个骰子掷出j ...

  4. HihoCoder1339 Dice Possibility(概率DP+母函数)

    时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 What is possibility of rolling N dice and the sum of the numb ...

  5. HDU 4652 Dice (概率DP)

    版权声明:欢迎关注我的博客,本文为博主[炒饭君]原创文章,未经博主同意不得转载 https://blog.csdn.net/a1061747415/article/details/36685493 D ...

  6. Dice (III) 概率dp

    #include <cstdio> #include <iostream> #include <cstring> #include <algorithm> ...

  7. Throwing Dice(概率dp)

    C - Throwing Dice Time Limit:2000MS     Memory Limit:32768KB     64bit IO Format:%lld & %llu Lig ...

  8. HDU 4599 Dice (概率DP+数学+快速幂)

    题意:给定三个表达式,问你求出最小的m1,m2,满足G(m1) >= F(n), G(m2) >= G(n). 析:这个题是一个概率DP,但是并没有那么简单,运算过程很麻烦. 先分析F(n ...

  9. SPOJ Favorite Dice(概率dp)

    题意: 一个骰子,n个面,摇到每一个面的概率都一样.问你把每一个面都摇到至少一次需要摇多少次,求摇的期望次数 题解: dp[i]:已经摇到i个面,还需要摇多少次才能摇到n个面的摇骰子的期望次数 因为我 ...

随机推荐

  1. 96D - Police Stations

    96D - Police Stations 思路:bfs,从所有的警察局开始bfs,因为bfs的深度一样,而且题目给的树保证满足条件,所以不用考虑深度. 如果搜索到一个点a,他的下一个点b已经被搜索过 ...

  2. Python爬虫Urllib库的高级用法

    Python爬虫Urllib库的高级用法 设置Headers 有些网站不会同意程序直接用上面的方式进行访问,如果识别有问题,那么站点根本不会响应,所以为了完全模拟浏览器的工作,我们需要设置一些Head ...

  3. npm 安装文件 运行报错 %1 is not a valid Win32 application

    安装了那个模板出了错报这样的错误 “%1 is not a valid Win32 application” 你就除那个模板新安装. 如下例: 运行 npm install -g @angular/c ...

  4. Windows Update 自动更新 设定 被锁(变灰)

    估计是McAfee自动更改掉的. 真TM烦人. 方法 1 不过找到了回复方法了: http://www.askvg.com/how-to-change-windows-update-settings- ...

  5. eclipse 与 tomcat 的那些路径

    我们用mvn创建了一个web工程,同时希望在eclipse里调试开发.mvn有mvn的路径要求,eclispe有eclipse的默认路径,怎么整合二者? 首先介绍一下eclipse的默认路径. 重点在 ...

  6. React 介绍

    ttps://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Function/bind The sm ...

  7. zzuli1728(数学期望,组合数)

    1728: 社交网络 Time Limit: 2 Sec  Memory Limit: 128 MBSubmit: 232  Solved: 64 SubmitStatusWeb Board Desc ...

  8. Oracle to_char()和to_date()函数的用法

    to_char()函数是我们经常使用的函数,下面就为您详细介绍Oracle to_date()函数的用法 1.to_char()函数分析 1)SQL中不区分大小写,MM和mm被认为是相同的格式代码 先 ...

  9. shiro学习笔记-Subject#login(token)实现过程

    本博文所有的代码均为shiro官网(http://shiro.apache.org/)中shiro 1.3.2版本中的源码. 追踪Subject的login(AuthenticationToken t ...

  10. Nginx在windows上安装 及 Nginx的配置及优化

    https://www.cnblogs.com/Chiler/p/8027167.html http://www.runoob.com/linux/nginx-install-setup.html 前 ...