Dadas N monedas, colocadas en una fila, indexadas 1 a N de izquierda a derecha. Inicialmente todas las monedas muestran cabeza. En cada turno, dos enteros, no necesariamente distintos, A y B entre 1 y N (inclusive) se eligen de manera uniforme al azar. Todas las monedas con un índice de A a B (inclusive) se voltean. ¿Cuál es el número esperado de monedas que muestran la cabeza después de que M gira?

Las monedas que mostrarán cabeza después de M turnos, serían las que se voltearon varias veces. Para cada moneda hay una cierta probabilidad asociada a que sea seleccionada / volteada en cualquier turno dado. Para la moneda [i] con índice i, esto es:
p = 2 * (i) / N * (N-i + 1) / N – 1 / N * 1 / N donde i = [1, N]
(el segundo término negativo es eliminar el exceso de conteo)

para cualquier moneda, la probabilidad de que se seleccione / voltee el número par de veces sería = mC0 (p) ^ 0 (1-p) ^ m + mC2 (p) ^ 2 (1-p) ^ m-2 + … que podría calcularse utilizando la fórmula de expansión binomial:
(p + 1-p) ^ m = mC0 p ^ 0 (1-p) ^ m + mC1 p ^ 1 (1-p) ^ m-1…. + mCm p ^ m (1-p) ^ 0
solo necesitamos la suma de términos con mC0, mC2, mC4 … es decir, términos pares, primero reemplace p + 1-p con -p + 1-p para obtener:
( -p + 1-p) ^ m = mC0 ( -p ) ^ 0 (1-p) ^ m + mC1 ( -p ) ^ 1 (1-p) ^ m-1…. + mCm ( -p ) ^ m (1-P) ^ 0
agregue estas dos series, los términos impares se cancelarían:
mC0 (p) ^ 0 (1-p) ^ m + mC2 (p) ^ 2 (1-p) ^ m-2 +… = 1/2 + 1/2 * (1-2p) ^ m
así que para cualquier índice “i”, P (moneda de lanzamiento uniforme [i]) = P (EFi) = 1/2 + 1/2 * (1-2p (i)) ^ m

Número esperado de moneda mostrando cara = Número esperado de moneda con lanzamiento uniforme
(utilizando la linealidad de la expectativa E (A + B + C + ..) = E (A) + E (B) + E (C) + …)
= P (EF1) + P (EF2) +…. P (EFN)
= 1/2 * N + 1/2 * [(1-2p (1)) ^ m + (1-2p (2)) ^ m +…. 1-2 (p (N)) ^ m)]

donde p [i] = 2 * (i) / N * (N-i + 1) / N – 1 / N * 1 / N

No es demasiado difícil ver que a medida que m se incrementa a inf, la probabilidad de lanzar una moneda de forma pareja es de 1/2, ya que en esa etapa es igualmente probable seleccionar una moneda impar o par de veces … en ese caso, se espera El número de monedas con cabeza sería N / 2.

La respuesta a su problema es:
[matemática] E (X) = \ frac {1} {2} \ sum \ limites_ {i = 1} ^ {n} (1+ (1-2p_ {i}) ^ {m}) [/ matemática]
Donde [matemáticas] p_i = \ frac {2i (n-i + 1) -1} {n ^ 2} [/ matemáticas]

Para ver esto, fijemos nuestra atención en una moneda específica en el índice [math] i [/ math].

¿Cuál es la probabilidad de que esta moneda en particular en el índice [matemáticas] i [/ matemáticas] se voltee?

Bueno, para que la moneda se voltee, debe estar entre nuestros índices elegidos al azar [matemática] A [/ matemática] y [matemática] B [/ matemática]. Ahora fije su atención en estos índices [matemáticas] A [/ matemáticas] y [matemáticas] B [/ matemáticas]. Tenga en cuenta que estamos interesados ​​en [matemáticas] 2 [/ matemáticas] opciones de [matemáticas] A [/ matemáticas] y [matemáticas] B [/ matemáticas] es decir

  1. [matemáticas] A \ le i [/ matemáticas] y [matemáticas] B \ ge i [/ matemáticas]
  2. [matemáticas] B \ le i [/ matemáticas] y [matemáticas] A \ ge i [/ matemáticas]

Cada uno ocurre con probabilidad [matemática] \ frac {i (n-i + 1)} {n ^ 2} [/ matemática] (puede multiplicar las probabilidades). Pero también tenga en cuenta que para encontrar la unión de estos dos eventos, también debemos restar su intersección (es decir, principio de exclusión de inclusión)

Por lo tanto, [matemáticas] p_i = \ frac {2i (n-i + 1)} {n ^ 2} – \ frac {1} {n ^ 2} = \ frac {2i (n-i + 1) -1} {n ^ 2} [/ matemáticas]

Ahora que sabemos esto, definamos [matemática] H (i, j) [/ matemática] como la probabilidad de que la moneda [matemática] i [/ matemática] sea cara después de que [matemática] j [/ matemática] se voltee.

Tenga en cuenta que para que la moneda [math] i [/ math] sea cara después de que [math] j [/ math] se lance, debe caer en uno de estos dos casos

  1. La moneda [matemática] i [/ matemática] ya estaba en la cara de [matemática] (j-1) [/ matemática] y se mantuvo así (es decir, no se volcó)
  2. La moneda [math] i [/ math] era la cola en el [math] (j-1) [/ math] th flip y se cambió a head por el [math] j [/ math] th flip

Esto nos da la siguiente recurrencia
[matemáticas] H (i, j) = H (i, j-1) (1-p_i) + [/ matemáticas] [matemáticas] (1-H (i, j-1)) p_i [/ ​​matemáticas]
donde [matemáticas] H (i, 0) = 1 [/ matemáticas] (obviamente)

Tenga en cuenta que podemos multiplicar las probabilidades como se indica arriba porque estamos lidiando con un proceso sin memoria .

Podemos resolver esta recurrencia fijando [math] i [/ math] constante y obtenemos
[matemáticas] H (i, j) = \ frac {1} {2} (1+ (1-2p_i) ^ j) [/ matemáticas]

Ahora para obtener nuestro valor esperado, definamos una variable aleatoria indicadora [matemática] X_i [/ ​​matemática] que sea [matemática] 1 [/ matemática] si la moneda [matemática] i [/ matemática] es cara y [matemática] 0 [/ matemáticas] de lo contrario

Entonces, el número esperado de caras es simplemente la suma de los valores esperados de estas variables aleatorias, que es simplemente la suma de sus probabilidades.

Por lo tanto, [matemáticas] E (X) = \ sum \ limites_ {i = 1} ^ {n} H (i, m) [/ matemáticas]
[matemática] = \ frac {1} {2} \ sum \ limits_ {i = 1} ^ {n} (1+ (1-2p_i) ^ m) [/ matemática] como se esperaba.

Esta es exactamente la misma respuesta obtenida por Abhishek Shekhawat aunque usando un enfoque diferente

Este es un problema del Proyecto Euler. No responderé

More Interesting

¿Por qué la gente siempre critica un doctorado en CS? ¿Por qué existe el título cuando bloquea las oportunidades profesionales?

¿Es necesaria la lógica equitativa en informática?

¿Cuáles son las diversas formas en que puede resolver el siguiente laberinto con un robot seguidor de enlace negro basado en IR? ¿Cómo puede resolverlo con el mínimo número de sensores posible y el tiempo más rápido para llegar al final?

¿Qué algoritmo debo usar para crear un solucionador de Sudoku?

¿Puede un programa de computadora derivar las matemáticas?

¿Cuáles son los cursos matemáticos recomendados para el aprendizaje automático y el procesamiento de big data?

¿Qué significa shift / reduce en el análisis?

¿Qué libros leerías para poder resolver el problema P versus NP por tu cuenta, recién salido de la escuela secundaria?

¿Podría la funcionalidad de una computadora digital ser duplicada por una computadora mecánica (con engranajes, ruedas, palancas, etc.)?

¿Quiénes son los mejores profesores que trabajan en algoritmos de aproximación?

¿Puede una resta dar un resultado negativo usando un número sin signo?

¿Existe alguna arquitectura de computadora basada en el cálculo lambda (en lugar de la máquina de Turing)?

¿Qué entero decimal está representado por 0xE4 en una notación de complemento a dos de 8 bits?

Informática: ¿Son nerds los estudiantes de informática?

¿La investigación colaborativa dificulta el uso de mejores herramientas (por ejemplo, TeXmacs en lugar de TeX / LaTeX)?