Probabilidade de que o arranjo do Papai Noel Secreto resulte em pares perfeitos

11

Então, nós tivemos o Papai Noel Secreto no trabalho.

Somos 8 pessoas. Cada um de nós se revezou e puxou um pequeno pedaço de papel de uma tigela com um nome. A única regra: se você digitar seu nome, precisará colocar o pedaço de papel de volta na tigela e tentar novamente.

Vamos chamar as pessoas de A, B, C, D, E, F, G, H, que também é a ordem na qual eles pegaram seu pedaço de papel.

Fizemos a troca de presentes ontem à noite.

A era o Papai Noel secreto de F.
B era o Papai Noel secreto de E.
C era o Papai Noel secreto de D.
D era o Papai Noel secreto de C.
E era o Papai Noel secreto de B.
F era o Papai Noel secreto de A.
G era o Papai Noel secreto de H.
Ele era o Papai Noel secreto de G.

Veja o que aconteceu? Nós fizemos casais.

A e F eram o Papai Noel secreto um do outro.
B e E eram o Papai Noel secreto um do outro.
C e D eram o Papai Noel secreto um do outro.
G e H eram o Papai Noel secreto um do outro.

Qual é a probabilidade disso acontecer e como você o calcula?

Hermann
fonte
1
"Se você digitar seu nome, precisará colocar o pedaço de papel de volta na tigela e tentar novamente." O que acontece se você é o último a escolher o seu próprio nome?
Juho Kokkala
Se a pessoa A desenha o rótulo C (digamos) e a pessoa B desenha o rótulo B, a pessoa A também coloca o rótulo C de volta no chapéu e desenha novamente? É isso que as respostas parecem sugerir, mas entendo que as palavras significam que A mantém os rótulos C e B redesenhados do chapéu que contém rótulos (A, B, D, E, F, G, H).
Juho Kokkala

Respostas:

14

O número total de tarefas entre pessoas, onde ninguém é designado para si, é d ( 2 n ) = ( 2 n ) ! ( 1 / 2 - 1 / 6 + + ( - 1 ) k / k ! + + 1 / ( 2 N ) ! ) . (Estes são chamados de desarranjos .) O valor está muito próximo de (2n

d(2n)=(2n)!(1/21/6++(1)k/k!++1/(2n)!).
.(2n)!/e

Se eles correspondem a pares perfeitos, então são um produto de transposições disjuntas . Isso implica que sua estrutura de ciclo é da forma

(a11a12)(a21a22)(an1an2).

O número de padrões distintos é a ordem do grupo de todas as permutações dos nomes divididos pela ordem do estabilizador do padrão. Um elemento estabilizador pode trocar qualquer número de pares e também pode permutar n ! pares, de onde existem 2 n n ! elementos estabilizadores. Portanto, existem2nn!2nn!

p(2n)=(2n)!2nn!

esses pares.

Como todos esses pares perfeitos são desarranjos, e todos os desarranjos são igualmente prováveis, a chance é igual a

p(2n)d(2n)=12nn!(11/2+1/6+(1)k/k!++1/(2n)!)e2nn!.

Para pessoas, por conseguinte, a resposta exacta é 15 / 2119 0,00707881 enquanto que a aproximação é de e / ( 2 42n=815/21190.00707881 : eles concordam com cinco números significativos.e/(244!)0.00707886


Para verificar, esta Rsimulação desenha um milhão de permutações aleatórias de oito objetos, retém apenas aqueles que são desarranjos e conta aqueles que são pares perfeitos. Ele gera sua estimativa, o erro padrão da estimativa e um escore Z para compará-lo ao valor teórico. Sua saída é

       p.hat           se            Z 
 0.006981031  0.000137385 -0.711721705

0.00660.0073

paired <- function(x) crossprod(x[x] - 1:length(x))==0
good <- function(x) sum(x==1:length(x)) == 0

n <- 8
set.seed(17)
x <- replicate(1e6, sample(1:n, n))
i.good <- apply(x, 2, good)
i.paired <- apply(x, 2, paired)

n.deranged <- sum(i.good)
k.paired <- sum(i.good & i.paired)
p.hat <- k.paired / n.deranged
se <- sqrt(p.hat * (1-p.hat) / n.deranged)
(c(p.hat=p.hat, se=se, Z=(p.hat - 15/2119)/se))
whuber
fonte
+1 para a cara de guaxinim boba e os óculos ... Peguei um atalho no conceito de "elemento estabilizador" porque não sei por onde começar a procurá-lo, mas faz muito sentido mesmo levar esse pouco intuitivamente.
Antoni Parellada
@Antoni Veja en.wikipedia.org/wiki/Burnside's_lemma, por exemplo.
whuber
1
@Amoeba, pensei em fazer isso, mas decidi manter o foco no problema atual, pois os distúrbios são muito conhecidos. O artigo da Wikipedia ao qual vinculei fornece vários métodos para derivar essa fórmula. O método mais evidente usa o Princípio da inclusão-exclusão, como é evidente na expressão de soma alternada.
whuber
1
Você supõe que o desenho dos rótulos é iniciado do zero se alguém desenha seu próprio rótulo (veja meus comentários à pergunta). Caso contrário, não acho que todas as perturbações sejam igualmente prováveis.
Juho Kokkala
1
@ Juho Essa é uma boa pergunta, digna de mais reflexão. Respondi com base na intenção implícita do procedimento de desenho, que seria criar todos os desarranjos com igual probabilidade, mas não está claro exatamente qual procedimento foi seguido ou se ele geraria desarranjos com uma distribuição uniforme (ou se é garantido mesmo para produzir um desarranjo!).
whuber
7

Fiquei bastante impressionado com a elegância na resposta @whuber. Para ser sincero, tive que me familiarizar com novos conceitos para seguir as etapas de sua solução. Depois de gastar muito tempo nisso, decidi postar o que consegui. Então, o que se segue é uma nota exegética à sua resposta já aceita. Dessa forma, não há tentativa de originalidade, e meu único objetivo é fornecer alguns pontos de ancoragem adicionais para seguir algumas das etapas envolvidas.

Então aqui vai ...

2n

2. Podemos derivar a fórmula para desarranjos?

n

d(n)=(n1)[d(n2)+d(n1)]=

=nd(n2)d(n2)+nd(n1)d(n1)

d(n)nd(n1)=[d(n1)(n1)d(n2)]

Agora, observando o paralelismo entre o LHS dessa equação e a parte do RHS entre parênteses, podemos continuar recursivamente:

d(n)nd(n1)=[d(n1)(n1)d(n2)]=

=(1)2[d(n2)(n2)d(n3)]==(1)n2d(2)2d(1)

d(n)=nd(n1)+(1)n

Trabalhando para trás:

d(2)=1

d(3)=3d(2)1=311

d(4)=4d(3)+1=4314+1

d(5)=5d(4)1=543154+51

d(6)=6d(5)+1=65431654+656+1=

=6!(12132+143215432+16!)=

=6!(16!15!+14!13!+12!11!+1)

Então, em geral,

d(n)=n!(11+12!13!+14!++1n!)

exx=1

d(n)n!e

a,b,c,d,e,fb,d,a,c,f,ea -> b -> d -> c after which it returns to ae -> f(a b d c)(e f)

4

(2n)!2n2nn!p(2n)=(2n)!2nn!


Para a Rsimulação:

1 paired <- function(x) crossprod(x[x] - 1:length(x))==0

x[x]8 vetor de elementos representando as atribuições atuais e determinar se ele é composto de loops de 2 elementos, como no problema do Papai Noel. Enquanto as permutações corresponderem a elementos de transposição, de modo que, se Paulo deveria dar um presente a Maria ( Paul -> Mariae vice-versa Maria -> Paul) e Max a John ( Max -> John/ John -> Max) inicialmente, a transposição resultante resultaria em um novo emparelhamento perfeito ( Max -> Maria/ Maria -> Maxe Paul -> John/ John -> Paul) estamos cumprindo a condição inicial de emparelhamento perfeito: insira a descrição da imagem aqui

Em outras palavras, se voltarmos ao exemplo dos chapéus em i 1

2) good <- function(x) sum(x==1:length(x)) == 0

x(1,2,3,4,5,6,7,8)

3.k.paired <- sum(i.good & i.paired) existe para excluir permutações emparelhadas como a acima no diagrama, que não são perturbações:

v <- c(1,2,3,4,5,6,7,8)
w <- c(1,2,3,5,4,6,7,8)

(c("is v paired?" = paired(v), "is w paired?" = paired(w),
   "is v a derang?" = good(w), "is w a derang?" = good(w)))

# not all paired permutations are derangements.
Antoni Parellada
fonte
1
e=
1
11
@whuber Obrigado. Eu realmente brinquei lá. Não sou bom em tarefas repetitivas e indexadas ... sabia que havia algo errado. Agora deve ser consertado.
Antoni Parellada