Lo siguiente es producto de leer en diagonal 
Creía que Samir buscaba hacer un modelo más sofisticado, que le diese la distribución de probabilidad sobre el valor total de la suma de éxitos y fracasos.
Hice un modelo de un juego similar hace años, cuando apenas sabía nada de matemáticas (de hecho empecé a estudiar matemáticas para poder hacer modelos de juegos probabilísticos como ése), así que puede estar mal:
https://www.geogebra.org/m/AT9hsghROtro juego parecido se discute aquí:
https://math.stackexchange.com/questions/2195516/rolling-dice-or-cumulative-binomial-distribution-with-a-twist/2195730Para muchos juegos de dados, dependiendo del número de dados, se puede sacar la probabilidad por "fuerza bruta" o, si la cosa se complica mucho, habría que buscar una aproximación analítica (si existiese), o aproximarle una después de una simulación de Monte-Carlo.
Para tu juego tenemos una serie de probabilidades multinomiales asociadas a una determinada composición de un número (creo que un número pequeño), me explico: hay una probabilidad de sacar un éxito, otra de sacar un fracaso y otra probabilidad digamos "neutra" que son probabilidades que en conjunto siguen una distribución trinomial.
Además de eso hay que considerar todas las sumas/restas posibles de números cero, uno o menos uno para dar lugar a un valor total determinado sobre el que calcular la probabilidad, esto es un tipo de composición de un número, en este caso "relativamente sencillo".
Por ejemplo, si tiramos X dados de 10 caras, ¿cuántas formas hay de obtener un total Y, y qué probabilidad hay asociada a cada forma?
Una forma de aproximar el problema es usando una función generadora, es decir, supongamos que tenemos 6 dados de 10 caras (con probabilidades iguales para sacar cada cara, y también se asume independencia de cada dado) donde sacar 6 ó más es un éxito y sacar un 1 un fracaso, entonces un dado se modela como el polinomio
\( \displaystyle f(x):=1\cdot x^{-1}+4\cdot x^0+5\cdot x^1 \)
Entonces la probabilidad de sacar un total de puntos \( T \) sería \( [x^T]\frac{f(x)^6}{f(1)^6} \), suponiendo que tiramos 6 dados, donde \( [x^T]f(x) \) significa "coeficiente de la potencia \( T \) en el polinomio \( f \)".
Como el polinomio es "feo" no creo que se puedan simplificar mucho las cuentas, es decir, usando una función generadora como la de arriba estamos usando un tipo de "fuerza bruta" (de hecho incluso aunque el polinomio fuese "simplificable" tampoco se iba a reducir demasiado el cálculo necesario, ya que se involucran sumas de coeficientes multinomiales), aunque las probabilidades dejadas son exactas. Por eso decía que cuando el número de dados (o combinaciones posibles) es alto es mejor usar directamente una simulación de Monte-Carlo y ya está.
Pero en este caso quizá no es práctico recurrir a funciones generadoras ni nada de eso ya que la composición del valor total es muy simple: que el total sea \( m>0 \) significa que el número de éxitos menos el de fracasos es justamente \( m \), con un número variable de valores neutros, así que posiblemente se pueda usar directamente una distribución trinomial para hallar estos valores de manera más sencilla, aunque computacionalmente no parece haber una gran ventaja respecto al uso de funciones generadoras, al menos para números de dados en juego pequeños.
