11 votos

Una pregunta muy difícil de probabilidad

En un cierto juego con 2 jugadores, el ganador es determinado por un laminado de solo 6 colindado mueren en turno, hasta que un 6 se muestra, al punto que el juego termina inmediatamente. Ahora, supongamos que k dice ahora se rodó simultáneamente por cada jugador en su turno, y el primer jugador para obtener un total de k (o más) 6, acumulado a lo largo de toda su lanza, gana el juego. (Por ejemplo, si k = 3, entonces el jugador 1 se lanzan 3 dados, y realizar el seguimiento de las 6 de la que se muestran. Si el jugador 1 no obtener todos los 6 entonces el jugador 2 va a hacer lo mismo. Suponiendo que el jugador 1 se presenta otra vez, volverá a tirar 3 dados, y cualquiera de los 6 que será añadido a su anterior total). Calcule la cantidad de vueltas que será necesario para completar el juego.

He analizado este problema de la siguiente manera: El problema puede ser modelada por una distribución binomial negativa con la probabilidad de $p=\frac{1}{6}$. Ahora, X es una variable aleatoria que representa el número de dados lanzados. Necesito encontrar el cdf $Pr[X\geq k]$, y, a continuación, encontrar la expectativa de la siguiente manera $E[X] = \int_0^\infty kPr[X\geq k]$. El problema aquí es que la cdf de una negativa binomal de distribución es una regularización de la función beta, y esto es bastante complicado de tratar. Me pregunto si hay otra forma de abordar este problema que no implican que?

8voto

Steve Kass Puntos 5967

Nota: La respuesta inicial aquí era incorrecta. Gracias a Markus, que resolvió el problema de una manera diferente y obtuvo diferentes resultados de mis primeras, he encontrado el error y reescribió la respuesta. (Si hay alguna razón para mí para volver a publicar el original, la respuesta es incorrecta en algún lugar, que me haga saber.)


Usted puede configurar una periodicidad para resolver si se consideran las etapas intermedias en el juego y considerar el número esperado de movimientos para el juego final a partir de ahí.

Deje $E(a,b,k)$ ser el esperado número de vueltas en el juego necesario para terminar si el jugador cuyo turno se necesita $a$ más de par de seis a ganar, el otro jugador necesidades de $b$ más, y $k$ dados son lanzados cada turno.

Si $a\le0$ o $b\le0$, el juego ya está terminado, y $E(a,b,k)=0$.

Ahora considere los posibles resultados de los próximos dos rollos. Deje que el número de $6$s en los rollos de ser$i$$j$, respectivamente.

La particular combinación de $(i,j)$ que sucede con la probabilidad de ${\tbinom{k}{i}}{\tbinom{k}{j}}\left(\frac{5}{6}\right)^{2k-i-j}\left(\frac{1}{6}\right)^{i+j}$. El esperado juego de la longitud de un caso determinado depende de si o no $i\ge a$: Si $i\ge a$, el juego termina después de la primera de las dos vueltas, y la (esperada) juego de la longitud de la es $1$. De lo contrario, es $2+E(a-i,b-j,k)$.

Esto conduce a la siguiente fórmula:

$$\begin{align} E(a,b,k)&= \sum_{i=0}^{a-1}\sum_{j=0}^k{\tbinom{k}{i}}{\tbinom{k}{j}}\left(\tfrac{5}{6}\right)^{2k-i-j}\left(\tfrac{1}{6}\right)^{i+j}\left(2+E(a-i,b-j,k)\right)\\ &+\sum_{i=a}^{k}{\tbinom{k}{i}}\left(\tfrac{5}{6}\right)^{k-i}\left(\tfrac{1}{6}\right)^{i},\\ \end{align}\\ \mbox{donde aquí y abajo } E(a,b,k)=0 \mbox{ si $a\le0$ o $b\le0$.} $$

Para obtener una útil recursiva fórmula, primero tire de la $i=0$, $j=0$ término de la suma.

$$ \begin{align}E(a,b,k) &= \left(\tfrac{5}{6}\right)^{2k}(2+E(a,b,k))\\ &+ \sum_{i={\color{red}{1}}}^{a-1}\sum_{j={\color{red}{0}}}^{{\color{red}{0}}}{\tbinom{k}{i}}{\tbinom{k}{j}}\left(\tfrac{5}{6}\right)^{2k-i-j}\left(\tfrac{1}{6}\right)^{i+j}\left(2+E(a-i,b-j,k)\right)\\ &+ \sum_{i=0}^{a-1}\sum_{j={\color{red}{1}}}^{k}{\tbinom{k}{i}}{\tbinom{k}{j}}\left(\tfrac{5}{6}\right)^{2k-i-j}\left(\tfrac{1}{6}\right)^{i+j}\left(2+E(a-i,b-j,k)\right)\\ &+\sum_{i=a}^{k}{\tbinom{k}{i}}\left(\tfrac{5}{6}\right)^{k-i}\left(\tfrac{1}{6}\right)^{i}\\ \end{align}$$

Luego resuelve $E(a,b,k)$. $$\begin{align} \left({1 - \left(\tfrac{5}{6}\right)^{2k}}\right)E(k,k,k) &= 2\left(\tfrac{5}{6}\right)^{2k}\\ &+ \sum_{i=1}^{a-1}{\tbinom{k}{i}}\left(\tfrac{5}{6}\right)^{2k-i}\left(\tfrac{1}{6}\right)^{i}\left(2+E(a-i,b,k)\right)\\ &+ \sum_{i=0}^{a-1}\sum_{j=1}^{k}{\tbinom{k}{i}}{\tbinom{k}{j}}\left(\tfrac{5}{6}\right)^{2k-i-j}\left(\tfrac{1}{6}\right)^{i+j}\left(2+E(a-i,b-j,k)\right)\\ &+\sum_{i=a}^{k}{\tbinom{k}{i}}\left(\tfrac{5}{6}\right)^{k-i}\left(\tfrac{1}{6}\right)^{i}\\ \end{align} $$

El uso de este para evaluar $E(k,k,k)$ en Mathematica (suponiendo que no hay errores tipográficos en mi transcripción de lo que he usado) da los mismos resultados que Markus poner en su respuesta. He aquí el comienzo de la secuencia de $E(k,k,k)$, a partir de $k=1$.

\begin{array}{ccc} \it{k} & \it{E(k,k,k)} \\ \hline \\ 1 & 6.00000000000 \\ 2 & 7.84416992031 \\ 3 & 8.69584550140 \\ 4 & 9.20010963516 \\ 5 & 9.54353874272 \\ 6 & 9.79634774936 \\ 7 & 9.99197435248\\ 8 & 10.1488645407\\ 9 & 10.2781431548\\ \cdots\\ 30 & 11.2061391171\\ \cdots\\ \end{array}

Estos valores deben ser correctos para todos los decimales, ya que no hay aproximación va a excepción de la final de la conversión de una fracción a decimal.

6voto

Markus Scheuer Puntos 16133

Nota: Esta pregunta parecía ser un buen reto para mí. Aquí doy una fórmula explícita, que realmente no es muy atractivo, porque es bastante complicado. De hecho, me gustó mucho la metodología elegante de Steve y yo aún lo hacen. Pero, mis cálculos con la fórmula explícita de mostrar diferentes valores de $k>1$, y la cosa se puso interesante para mí. He comprobado con un pequeño programa en el que las cifras presentadas por $E(k,k,k)$ de Steve fueron consistentes con su fórmula de recursión. Pero sin embargo mi cifras eran diferentes a las de aquellos. Escribí pequeños programas de cálculo de primer y el último expresiones de mis cálculos para verificar, que las transformaciones que estaban en lo correcto. Luego hice una alternativa de verificación, escribió algunas código R para producir muestras aleatorias y simulada de los rollos con $1$ $6$dados. La noticia agradable para mí fue, que la simulación se basa en muestras aleatorias confirmado de manera convincente las cifras calculadas con mi fórmula explícita. Así que, para mí una pregunta abierta todavía está a la izquierda, es decir, ¿cuál es la razón para la diferencia de resultados entre el explícito fórmula y la fórmula de recursión?

@Steve, yo realmente apreciaría si usted puede revisar tanto las respuestas y ayudar a aclarar.

Añadido 2014-03-31: Nota: la Recursividad corregido! En el mientras tanto Steve encontró una inteligente de corrección de su respuesta, y ahora ambos, recursiva y solución directa muestran resultados consistentes.

Ahora, he aquí una prueba de lo que resulta en una fórmula explícita para la espera: Vamos a $A^{(k)}_{N,j}$ denotar el número de maneras reproductor $A$ puede conseguir $j(\ge0)$ veces $6$ $k$ dados en $N$ rollos. Deje $A^{(k)}(x,y)=(5x+y)^{Nk}$ la correspondiente generación de función, con el exponente de $y$ dando el número de $6$ después $N$ rollos con $k$ dados y el coeficiente de dar el número de diferentes posibilidades de este evento. El coeficiente de $x$ marca el número de posibilidades de todos los demás laminados en pips con $k$ dados. Deje que el coeficiente de operador $[y^j]$ el valor del coeficiente de $y^j$$A^{(k)}(x,y)$. Así,

\begin{align} A^{(k)}_{N,j}=[y^j](5x+y)^{Nk}=[y^j]\sum_{l=0}^{Nk}\binom{Nk}{l}y^l(5x)^{Nk-l}=\binom{Nk}{j}5^{Nk-j} \end{align}

Igualmente, os $A^{(k)}_{N,\leq j}$ denotar el número de posibilidades de reproductor $A$ para obtener no más de $j$ $6$ con $k$ dados en $N$ rollos.

$$A^{(k)}_{N,\leq j}=\sum_{l=0}^{j}\binom{Nk}{l}5^{Nk-l}$$ También tenemos $A^{(k)}_{N,\geq j}=6^{Nk} - A^{(k)}_{N,< j}$, ya que el $6^{Nk}$ es el número de todos los resultados posibles después de $N$ rollos con $k$ dados. Para el jugador $B$ utilizamos la correspondiente notación $B^{(k)}_{N,j}$, etc.

La resultante de las probabilidades de estos eventos son: \begin{align} p(A^{(k)}_{N,j})&=p(B^{(k)}_{N,j})=\left(\frac{5}{6}\right)^{Nk}\binom{Nk}{j}\frac{1}{5^j}\\ p(A^{(k)}_{N,\leq j})&=p(B^{(k)}_{N,\leq j})=\left(\frac{5}{6}\right)^{Nk}\sum_{l=0}^{j}\binom{Nk}{l}\frac{1}{5^l}\\ p(A^{(k)}_{N,\geq j})&=p(B^{(k)}_{N,\geq j})=1-p(A^{(k)}_{N,< j})\\ \end{align}

Para calcular el $E^{(k)}(X)$, la expectativa de cuando se juega con $k$ dados, vamos a $P(X=N)$ denotar la probabilidad de que Un jugador o el jugador B gana después de un total de $N$ rollos. En orden a ello, consideramos dos escenarios. Cualquiera de las $A$ ha rodado $k$ (o más) $6$ con su $(N+1)$-st rollo de $k$ dados, de modo que $A$ gana después de un total de $2N+1$ rollos, o $B$ gana después de la $A$ ha lanzado la $k$ dados $N$ veces, por lo que han llegado a un total de $2N$ veces. Simbólicamente

\begin{align} A\ wins:\ \ \ &(A B) (A B)\ldots (A B) (A B)A = (AB)^{N}A&(2N+1)\ rolls\\ B\ wins:\ \ \ &(A B) (A B)\ldots (A B) (A B) = (AB)^{N}&(2N)\ rolls \end{align}

Así, la expectativa $E^{(k)}(X)$ se compone de dos partes. El primer sumando corresponde a la situación con el jugador de la $A$ como ganador de la partida. Esto implica que foreach $N\ge0$ player $B$ ha alcanzado menos de $k$ $6$ en $N$ rollos, mientras que el jugador $A$ ha alcanzado $j<k$ veces $6$ $N$ rollos y con su $(N+1)$-st rollo ha $k$ o más $6$. Del mismo modo la segunda suma con $B$ como el ganador.

\begin{align} E^{(k)}(X)&=\sum_{N\ge0}NP(X=N)\\ &=\sum_{N\ge0}(2N+1)P(X=2N+1) + \sum_{N\ge0}(2N)P(X=2N)\\ &=\sum_{N\ge0}(2N+1)p(B_{N,<k})\sum_{j=0}^{k-1}p(A_{N,j})p(A_{1,\ge k-j})\\ &+\sum_{N\ge0}(2N)p(A_{N,<k})\sum_{j=0}^{k-1}p(B_{N-1,j})p(B_{1,\ge k-j}) \end{align}

La sustitución de las probabilidades de las fórmulas indicadas arriba da

\begin{align} E^{(k)}(X)&=\sum_{N\ge0}(2N+1) \left(\sum_{i=0}^{k-1}\binom{Nk}{i}\left(\frac{5}{6}\right)^{Nk}\frac{1}{5^i}\right)\\ &\qquad\qquad\cdot\sum_{j=0}^{k-1}\left(\binom{Nk}{j}\left(\frac{5}{6}\right)^{Nk}\frac{1}{5^j} \sum_{l=k-j}^{k}\binom{k}{l}\left(\frac{5}{6}\right)^{k}\frac{1}{5^l}\right)\\ &+\sum_{N\ge0}(2N) \left(\sum_{i=0}^{k-1}\binom{Nk}{i}\left(\frac{5}{6}\right)^{Nk}\frac{1}{5^i}\right)\\ &\qquad\qquad\cdot\sum_{j=0}^{k-1}\left(\binom{(N-1)k}{j}\left(\frac{5}{6}\right)^{(N-1)k}\frac{1}{5^j} \sum_{l=k-j}^{k}\binom{k}{l}\left(\frac{5}{6}\right)^{k}\frac{1}{5^l}\right)\\ &=\sum_{N\ge0}(2N+1) \left(\frac{5}{6}\right)^{(2N+1)k}\sum_{i=0}^{k-1}\binom{Nk}{i}\frac{1}{5^i} \sum_{j=0}^{k-1}\binom{Nk}{j}\frac{1}{5^j} \sum_{l=k-j}^{k}\binom{k}{l}\frac{1}{5^l}\\ &+\sum_{N\ge0}(2N) \left(\frac{5}{6}\right)^{2Nk}\sum_{i=0}^{k-1}\binom{Nk}{i}\frac{1}{5^i} \sum_{j=0}^{k-1}\binom{(N-1)k}{j}\frac{1}{5^j} \sum_{l=k-j}^{k}\binom{k}{l}\frac{1}{5^l}\\ \end{align}

Después de algunas simplificaciones, el valor resultante para $E^{(k)}(x)$ es \begin{align} E^{(k)}(X)&=\frac{1}{6^k}\sum_{N\ge0}(2N+1) \left(\frac{5}{6}\right)^{2Nk}\sum_{i=0}^{k-1}\binom{Nk}{i}\frac{1}{5^i} \sum_{j=0}^{k-1}\binom{Nk}{j} \sum_{l=0}^{j}\binom{k}{j-l}\frac{1}{5^l}\\ &+\frac{1}{6^k}\sum_{N\ge0}(2N) \left(\frac{5}{6}\right)^{(2N-1)k}\sum_{i=0}^{k-1}\binom{Nk}{i}\frac{1}{5^i} \sum_{j=0}^{k-1}\binom{(N-1)k}{j} \sum_{l=0}^{j}\binom{k}{j-l}\frac{1}{5^l} \end{align}

He calculado las expectativas de $k=1\ldots 6$ con un pequeño programa. Los valores calculados aproximar $E^{(k)}(x)$ convergen rápidamente con el aumento de $N$. La tabla siguiente muestra las cifras resultantes con $7$ dígitos significativos.

Expectativa $E^{(k)}(X)$: El valor en la columna $N$ indica que el valor más pequeño, donde los dígitos significativos ya no cambio por el aumento de $N$. $$ \begin{array}{ccc} \it{k} & \it{E^{(k)}(X)} & \it{N}\\ \hline \\ 1 & 6.000 000 & 57 \\ 2 & 7.844 169 & 32 \\ 3 & 8.695 845 & 26 \\ 4 & 9.200 109 & 22 \\ 5 & 9.543 538 & 19 \\ 6 & 9.796 347 & 18 \\ \end{array} $$

En la siguiente me dan valores explícitos para $k=1$$2$. Con la ayuda de funciones de generación y $D$ (formal) operador diferencial que tenemos para $k=1$:

\begin{align} E^{(1)}(x)&=\frac{1}{6}\sum_{N\ge0}(2N+1)\left(\frac{5}{6}\right)^{2N} + \frac{1}{6}\sum_{N\ge0}(2N)\left(\frac{5}{6}\right)^{2N-1}\\ &=\frac{1}{6}\left.\left(D\frac{1}{1-x}\right)\right\rvert_{x=\frac{5}{6}}=\frac{1}{6}\left.\frac{1}{(1-x)^2}\right\rvert_{x=\frac{5}{6}}\\ &=6 \end{align}

Cálculo de $k=2$ le da:

\begin{align} E^{(2)}(x)&=\frac{1}{6^2}\sum_{N\ge0}(2N+1)\left(\frac{5}{6}\right)^{4N} \sum_{i=0}^{1}\binom{2N}{i}\frac{1}{5^i} \sum_{j=0}^{1}\binom{2N}{j}\sum_{l=0}{j}\binom{2}{j-l}\frac{1}{5^l}\\ &+\frac{1}{6^2}\sum_{N\ge0}(2N)\left(\frac{5}{6}\right)^{4N-2} \sum_{i=0}^{1}\binom{2N}{i}\frac{1}{5^i} \sum_{j=0}^{1}\binom{2N-2}{j}\sum_{l=0}{j}\binom{2}{j-l}\frac{1}{5^l}\\ &=\frac{1}{5^26^2}\sum_{N\ge0}(2N+1)(2N+5)(22N+5)\left(\frac{5}{6}\right)^{4N}\\ &+\frac{1}{5^26^2}\sum_{N\ge0}(2N)(2N+5)(22N-17)\left(\frac{5}{6}\right)^{4N-2}\\ &=\frac{1}{5^46^2}\sum_{N\ge0}(5368N^3+12572N^2-1870N+625)\left(\frac{5}{6}\right)^{4N}\\ \end{align}

Utilizamos el siguiente bien conocida de las relaciones:

\begin{align} \sum_{N\ge0}Nx^N&=(xD)\frac{1}{1-x}=\frac{x}{(1-x)^2}\\ \sum_{N\ge0}N^2x^N&=(xD)^2\frac{1}{1-x}=\frac{x(x+1)}{(1-x)^3}\\ \sum_{N\ge0}N^3x^N&=(xD)^3\frac{1}{1-x}=\frac{x(x^2+4x+1)}{(1-x)^4}\\ \end{align}

para obtener: \begin{align} E^{(2)}(x)&=\frac{1}{5^46^2}\left.\left(5368\frac{x(x^2+4x+1)}{(1-x)^4} +12572\frac{x(x+1)}{(1-x)^3}\right.\right.\\ &\qquad\qquad\qquad\left.\left.-1870\frac{x}{(1-x)^2} +625\frac{1}{1-x}\right)\right|_{x=\left(\frac{5}{6}\right)^4}\\ &=\frac{1}{5^46^2}\left.\frac{-9699x^3+27087x^2+14195x+625}{(1-x)^4}\right|_{x=\left(\frac{5}{6}\right)^4}\\ &=\frac{636876}{81191}\approx7.844169 \end{align}

Resumen de la expectativa $k=1,2$: \begin{align} E^{(1)}(x)&=6\\ E^{(2)}(x)&=\frac{636876}{81191}\approx7.844169 \end{align}

Las fórmulas para $k>2$ puede ser calculado de una manera similar, pero tenga en cuenta que el grado de $N$ aumenta por $2$ siempre $k$ se incrementa por $1$. Así, el cálculo manual podría convertirse realmente engorroso. :-)

0voto

JSS Puntos 1

Podría ser más fácil de abordar este problema con la generación de funciones.

Cuando se lanza un dado, usted tiene 1/6 oportunidad de rodar un 6, y una 5/6 oportunidad de rodar cualquier otro de la cara. Vamos a representar a la rodadura de un 6, y deje b representan el balanceo de cualquier otra cara.

Tirar un solo dado, los resultados pueden ser representados como $({a\over 6} + {5b\over 6})$

Para k dice, tenemos $p(a)=({a\over 6} + {5b\over 6})^k$

Cada monomio de la expansión será de la forma $a^{r_1}b^{r_2}$, de tal manera que $r_1+r_2=k$. El coeficiente da la probabilidad de obtener exactamente $r_1$ 6 $r_2$ no-6. Ya que no ganan puntos cuando no lo hacemos rodar un 6, podemos ignorar estos resultados mediante el establecimiento $b=1$.

$p(a)=({a\over 6} + {5\over 6})^k$

Podemos obtener el número esperado de 6 por turno por primera diferenciación con respecto a una, y luego de evaluar a=1. Si queremos puntuación w puntos para ganar, dividimos por la espera puntos por turno para obtener el número esperado de personal se convierte en un jugador debe esperar antes de ganar.

Si usted quiere cuántos acumulada se convierte debemos esperar a completar el juego: si cada jugador espera $t$ vuelve a ganar, y un jugador va primero, es de suponer que el juego debería ir $2t-1$ vueltas en total para completar. Alternativamente, manipular la generación de la función de la cuenta que se espera de los totales acumulados.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X