El calendario gregoriano favorece a cinco de los siete días de la semana durante los años bisiestos. Por lo tanto, la posibilidad no es precisamente .2/7
Este fue esencialmente el problema B3 en la competencia de matemáticas de Putnam de 1950 :
n se elige al azar de los números naturales. Demuestre que la probabilidad de que el 25 de diciembre en el año sea un miércoles no es 1/7.n
En el calendario gregoriano , los años que son múltiplos de son años bisiestos (con días), pero los años que son múltiplos de no son años bisiestos (y por lo tanto tienen días), con la excepción de que los años que son múltiplos de son años bisiestos. (Muchos de nosotros recordamos la excepción más reciente en ). Esto crea un ciclo de años que contiene años bisiestos.47×52+2=3661007×52+1=3654002000400400 / 4 - 400 / 100 + 400 / 400 = 97
Lo que es especialmente interesante es que el número total de días en este ciclo es un múltiplo entero de siete:
400 × ( 7 × 52 + 1 ) + 97 × 1 ≡ 400 + 97 ≡ 71 × 7 ≡ 0 mod7)
Esto muestra que el ciclo de años comprende un número entero de semanas. En consecuencia, el patrón de días de la semana es exactamente el mismo de un ciclo a otro.400
Por lo tanto, podemos interpretar la pregunta como la posibilidad de domingos al tomar muestras de manera aleatoria y uniforme de cualquier ciclo de años bisiestos. Un cálculo de la fuerza bruta (usando, por ejemplo, el hecho de que el 1 de enero de 2001 era un lunes) muestra que de los años bisiestos en cada ciclo tienen domingos. Por lo tanto, la posibilidad es53400289753
Pr(53 Sundays)=2897.
Tenga en cuenta que esto no es igual a : es ligeramente mayor. Por cierto, existe la misma probabilidad de miércoles, viernes, sábados o lunes y solo una probabilidad de de martes o jueves.28/98=2/75327/9753
Para aquellos que desean hacer cálculos más detallados (y pueden desconfiar de cualquier simplificación matemática), aquí hay un código de fuerza bruta que computa y examina todos los días de la semana durante un conjunto de años. Al final muestra el número de años con apariciones de cada día de la semana. Está escrito en . 53R
Aquí está su salida para el ciclo :2001−2400
Friday Monday Saturday Sunday Thursday Tuesday Wednesday
28 28 28 28 27 27 28
Aquí está el código en sí.
leapyear <- function(y) {
(y %% 4 == 0 & !(y%% 100 == 0)) | (y %% 400 == 0)
}
leapyears <- seq(2001, length.out=400)
leapyears <- leapyears[leapyear(leapyears)]
results <- sapply(leapyears, function(y) {
table(weekdays(seq.Date(as.Date(paste0(y, "-01-01")), by="1 day", length.out=366)))
})
rowSums(results==53)