题目

太郎和一只免子正在玩一个掷骰子游戏。有一个有N个格子的长条棋盘,太郎和兔子轮流掷一个有M面的骰子,骰子M面分别是1到M的数字.且掷到任意一面的概率是相同的.掷到几.就往前走几步.当谁走到第N格时,谁就获胜了。游戏中还有一个规则“反弹”.就是当一位选手要走到第N格外时.他就会后退(就像飞行棋进营一样)。

假设现在一位追手在A格.当他掷出B时:

1.A+B<N,走到第A+B.络,

2.A+B=N,走到第N格,获胜。

3.A+B≥N,走到第(N-(A+B-N)格

现在太郎和兔子分别在第x和y格.接下来是太郎掷骰子,太郎想知道他赢得比赛的概率就多少。

分析

设\(f_{i,j}\)表示太郎在i,兔子在j,太郎的胜率。我们从后往前转移。

我们分四种情况:

1、i+m<=n and j+m<=n
2、i+m>n and j+m<=n
3、i+m<=n and j+m>n
4、i+m>n and j+m>n

why?

因为发现,当\(i+m>n\)时,i怎么跳总是\(i+m>n\),那么就可以把它们当做同一种状态。j也一样。

情况一:i+m<=n and j+m<=n

因为i走到k的概率为\(\dfrac{1}{m}\),j走到l的概率也为\(\dfrac{1}{m}\),

那么状态(i,j)的胜率就是状态(k,l)胜率的总和。

\[f_{i,j}=\dfrac{1}{m^2}\sum_{k=i+1}^{i+m+1}\sum_{l=j+1}^{j+m+1}f_{k,l}
\]

情况二:i+m>n and j+m<=n

\[f_{i,j}=\dfrac{(m-1)\sum_{l=j+1}^{j+m+1}f_{i,l}}{m^2}+\dfrac{1}{m}
\]

显然i有\(\dfrac{1}{m}\)的概率到达终点,而i有m种可能,又那么既然已经算了到达终点的概率,那么就不用在计算,所以乘以(m-1)。

情况三:i+m<=n and j+m>n

\[f_{i,j}=\dfrac{(m-1)\sum_{k=i+1}^{i+m+1}f_{k,j}}{m^2}+\dfrac{x(如果i=n-m,x=1,否则x=0)}{m}
\]

同样j有m种可能,但不能让他到达终点,那么就不用在计算,所以乘以(m-1)。

而当i在n-m这个位置时,i也有\(\dfrac{1}{m}\)的概率到达终点,有可能出现状态(n,n),由于太郎是先手,所以算太郎赢,而前面有减去了j到达终点的情况,所以加上去。

情况四:i+m>n and j+m>n

要求i赢,所以

当i第一回合就走到了n,概率为\(\dfrac{1}{m}\),

当i第一回合没有走到n,而j也不能走到n,在第二回合i走到了n概率为\(\dfrac{1}{m}(1-\dfrac{1}{m})^2\)

如果在第二回合i还是没有走到n,而j还是不能走到n,在第三回合i走到了n概率为\(\dfrac{1}{m}(1-\dfrac{1}{m})^4\)

如此类推,

\[f_{i,j}=limit_{n\to\infty}\dfrac{1}{m}+\dfrac{1}{m}(1-\dfrac{1}{m})^2+\dfrac{1}{m}(1-\dfrac{1}{m})^4+...+\dfrac{1}{m}(1-\dfrac{1}{m})^{2n}
\]

\[=\dfrac{1}{m}(1+(1-\dfrac{1}{m})^2+(1-\dfrac{1}{m})^4+...)+(1-\dfrac{1}{m})^{2n}
\]

等比数列求和

\[=\dfrac{1}{m}\dfrac{1*[1-(1-\dfrac{1}{m})^{2n}]}{1-(1-\dfrac{1}{m})^2}
\]

因为数列的公比小于1,\([1-(1-\dfrac{1}{m})^{2n}]\)无限趋近于1,所以

\[=\dfrac{1}{m}\dfrac{1}{1-(1-\dfrac{1}{m})^2}
\]

解得

\[f_{i,j}=\dfrac{m}{2m-1}
\]

但是这样是\(O(n^4)\)的,用矩阵后缀和优化,变成\(O(n^2)\)。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
const int mo=1000000007;
const int N=2005;
using namespace std;
double f[N][N],n,m,x,y;
double val(int x1,int y1,int x2,int y2)
{
return f[x1][y1]-f[x1][y2]-f[x2][y1]+f[x2][y2];
}
int main()
{
scanf("%lf%lf%lf%lf",&n,&m,&x,&y);
for(int i=n;i>=1;i--)
for(int j=n;j>=1;j--)
{
f[i][j]=f[i+1][j]+f[i][j+1]-f[i+1][j+1];
if(i==n)
{
f[i][j]++;
continue;
}
if(j==n) continue;
if(i+m>n && j+m>n)
f[i][j]+=m/(2*m-1);
else
if(i+m>n)
f[i][j]+=val(i,j+1,i+1,j+m+1)*(m-1)/m/m+1/m;
else
if(j+m>n)
f[i][j]+=(m-1)*val(i+1,j,i+1+m,j+1)/m/m+(i==n-m)/m/m;
else
f[i][j]+=val(i+1,j+1,i+m+1,j+m+1)/(m*m); }
printf("%.6lf",val(x,y,x+1,y+1));
}

【NOIP2012模拟10.31】掷骰子的更多相关文章

  1. 【JZOJ4835】【GDOI2017模拟10.31】量化交易

    题目描述 数据范围 解法 贪心: 从左往右枚举,设枚举到元素为x,并维护一个堆: 设此时堆顶元素为y, 如果x大于y,那么x可以与y产生差价,立即将差价贡献给答案. 如果y之前已经和其他元素z产生过差 ...

  2. 【NOIP2012模拟10.25】旅行

    题目 给定一个n行m列的字符矩阵,'.'代表空地,'X'代表障碍.移动的规则是:每秒钟以上下左右四个方向之一移动一格,不能进入障碍. 计算:在空地中随机选择起点和终点(可以重合,此时最短耗时为0),从 ...

  3. 【NOIP2012模拟10.25】剪草

    题目 有N棵小草,编号0至N-1.奶牛Bessie不喜欢小草,所以Bessie要用剪刀剪草,目标是使得这N棵小草的高度总和不超过H.在第0时刻,第i棵小草的高度是h[i],接下来的每个整数时刻,会依次 ...

  4. 【NOIP2012模拟10.25】单元格

    题目 在一个R行C列的表格里,我们要选出3个不同的单元格.但要满足如下的两个条件: (1)选中的任意两个单元格都不在同一行. (2)选中的任意两个单元格都不在同一列. 假设我们选中的单元格分别是:A, ...

  5. LeetCode 1223. 掷骰子模拟 Dice Roll Simulation - Java - DP

    题目链接:1223. 掷骰子模拟 有一个骰子模拟器会每次投掷的时候生成一个 1 到 6 的随机数. 不过我们在使用它时有个约束,就是使得投掷骰子时,连续 掷出数字 i 的次数不能超过 rollMax[ ...

  6. Python绘制直方图 Pygal模拟掷骰子

    #coding=utf-8 from random import randint class Die(): """骰子类""" def __ ...

  7. Python 使用matplotlib模块模拟掷骰子

    掷骰子 骰子类 # die.py 骰子类模块 from random import randint class Die(): """骰子类""&quo ...

  8. 使用python实现模拟掷骰子数据分析

    Data:2020/4/8 主题:模拟实现掷骰子数据分析 编译环境:pycharm 库:pygal 说明: code 1:创建一个掷骰子类对象,类方法获得掷骰子随机数1-6,默认6个面,模拟20次将结 ...

  9. Python Tkinter小实例——模拟掷骰子

    什么是Tkinter? Tkinter 是 Python 的标准 GUI 库.Python 使用 Tkinter 可以快速的创建 GUI 应用程序. 由于 Tkinter 是内置到 python 的安 ...

随机推荐

  1. Hibernate初了解

    Hibernate是一个开放源代码的对象关系映射框架, 它对JDBC进行了非常轻量级的对象封装,使得Java程序员可以随心所欲的使用对象编程思维来操纵数据库. Hibernate可以应用在任何使用JD ...

  2. mnist数据集下载——mnist数据集提供百度网盘下载地址

    mnist数据集是由深度学习大神 LeCun等人制作完成的数据集,mnist数据集也常认为是深度学习的“ Hello World!”. 官网:http://yann.lecun.com/exdb/mn ...

  3. OUTLOOK、foxmail等无法直接打开邮件中的超级链接问题

         部分电脑,在OUTLOOK或Foxmail收到隔离邮件通知时,点击发送或删除时,提示“一般性错误,*******************,找不到应用程序”.或打开其它HTML格式的邮件正文中 ...

  4. 【VS开发】【智能语音处理】特定人语音识别算法—DTW算法

    DTW(动态时间弯折)算法原理:基于动态规划(DP)的思想,解决发音长短不一的模板匹配问题.相比HMM模型算法,DTW算法的训练几乎不需要额外的计算.所以在孤立词语音识别中,DTW算法仍得到广泛的应用 ...

  5. hive排错

    找出错的那个hive实例,看错误日志: 点下面stdout,找Error

  6. CMDB 理论

    TIL即IT基础架构库(Information Technology Infrastructure Library, ITIL,信息技术基础架构库)由英国政府部门CCTA(Central Comput ...

  7. 初相识|performance_schema全方位介绍

    初相识|performance_schema全方位介绍 |导 语 很久之前,当我还在尝试着系统地学习performance_schema的时候,通过在网上各种搜索资料进行学习,但很遗憾,学习的效果并不 ...

  8. linux:文件/目录权限设置

    一.权限修改命令:     chmod [OPTION]... MODE[,MODE]... FILE...二.文件/目录MODE指定方式: 1).八进制(r=4,w=2,x=1) [root@pin ...

  9. [Vue] vue的一些面试题

    1.v-model 的原理 v-model 是一个语法糖,它即可以支持原生表单元素,也可以支持自定义组件.v-model 在内部为不同的输入元素使用不同的属性并抛出不同的事件. text 和 text ...

  10. jackson json序列化 首字母大写 第二个字母需小写

    有这样一个类: @Setter @Getter @JsonNaming(value = PropertyNamingStrategy.UpperCamelCaseStrategy.class) pub ...