Estatísticas de permutação aleatória - Random permutation statistics

As estatísticas de permutações aleatórias , como a estrutura do ciclo de uma permutação aleatória, são de fundamental importância na análise de algoritmos , especialmente de algoritmos de ordenação, que operam em permutações aleatórias. Suponha, por exemplo, que estejamos usando quickselect (um primo de quicksort ) para selecionar um elemento aleatório de uma permutação aleatória. O Quickselect realizará uma classificação parcial no array, pois ele particiona o array de acordo com o pivô. Conseqüentemente, uma permutação ficará menos desordenada depois que a seleção rápida for realizada. A quantidade de desordem que permanece pode ser analisada com funções geradoras. Essas funções geradoras dependem fundamentalmente das funções geradoras das estatísticas de permutação aleatória. Portanto, é de vital importância computar essas funções geradoras.

O artigo sobre permutações aleatórias contém uma introdução às permutações aleatórias.

A relação fundamental

Permutações são conjuntos de ciclos rotulados. Usando o caso rotulado do teorema fundamental Flajolet-Sedgewick e escrevendo para o conjunto de permutações e para o conjunto singleton, temos

Traduzindo em funções geradoras exponenciais (EGFs), temos

onde usamos o fato de que o EGF das espécies combinatórias de permutações (há n ! permutações de n elementos) é

Esta única equação permite derivar um grande número de estatísticas de permutação. Em primeiro lugar, retirando os termos de , ou seja, exp, podemos restringir o número de ciclos que uma permutação contém, por exemplo, restringindo o EGF para obter permutações contendo dois ciclos. Em segundo lugar, observe que o EGF dos ciclos marcados, ou seja , é

porque existem k! / k ciclos rotulados. Isso significa que, eliminando termos dessa função geradora, podemos restringir o tamanho dos ciclos que ocorrem em uma permutação e obter um EGF das permutações contendo apenas ciclos de um determinado tamanho.

Em vez de remover e selecionar ciclos, também se pode colocar pesos diferentes em ciclos de tamanhos diferentes. Se é uma função de peso que depende apenas do tamanho k do ciclo e por brevidade escrevemos

definindo o valor de b para uma permutação como a soma de seus valores nos ciclos, então podemos marcar ciclos de comprimento k com u b ( k ) e obter uma função geradora de duas variáveis

Esta é uma função geradora "mista": é uma função geradora exponencial em ze uma função geradora comum no parâmetro secundário u. Diferenciando e avaliando em u  = 1, temos

Esta é a função geradora de probabilidade da expectativa de b . Em outras palavras, o coeficiente de nesta série de potências é o valor esperado de b nas permutações em , dado que cada permutação é escolhida com a mesma probabilidade .

Este artigo usa o operador de extração de coeficiente [ z n ], documentado na página para séries de potências formais .

Número de permutações que são involuções

Uma involução é uma permutação σ de forma que σ 2 = 1 sob a composição de permutação. Segue-se que σ pode conter apenas ciclos de comprimento um ou dois, ou seja, a função geradora exponencial g ( z ) dessas permutações é

Isso dá a fórmula explícita para o número total de involuções entre as permutações σ ∈  S n :

Dividindo por n ! produz a probabilidade de que uma permutação aleatória seja uma involução. Esses números são conhecidos como números de telefone .

Número de permutações que são m th raízes da unidade

Isso generaliza o conceito de uma involução. Uma m- ésima raiz da unidade é uma permutação σ de modo que σ m = 1 sob a composição de permutação. Agora, toda vez que aplicamos σ, movemos um passo em paralelo ao longo de todos os seus ciclos. Um ciclo de comprimento d aplicado d vezes produz a permutação de identidade em d elementos ( d pontos fixos) ed é o menor valor para isso. Portanto, m deve ser um múltiplo de todos os tamanhos de ciclo d , ou seja, os únicos ciclos possíveis são aqueles cujo comprimento d é um divisor de m . Conclui-se que o EGF g ( x ) dessas permutações é

Quando m = p , onde p é primo, isso simplifica para

Número de permutações de ordem exatamente k

Isso pode ser feito por inversão de Möbius . Trabalhando com o mesmo conceito da entrada anterior, notamos que as espécies combinatórias de permutações cuja ordem divide k é dada por

Traduzindo para funções geradoras exponenciais, obtemos o EGF de permutações cuja ordem divide k , que é

Agora podemos usar esta função geradora para contar permutações de ordem exatamente k . Seja o número de permutações em n cuja ordem é exatamente d e o número de permutações em n a contagem de permutação cuja ordem divide k . Então nós temos

Segue-se pela inversão de Möbius que

Portanto, temos o EGF

A contagem desejada é então fornecida por

Esta fórmula produz, por exemplo, para k  = 6 o EGF

com a sequência de valores começando em n  = 5

(sequência A061121 no OEIS )

Para k  = 8, obtemos o EGF

com a sequência de valores começando em n  = 8

(sequência A061122 no OEIS )

Finalmente, para k  = 12, obtemos o EGF

com a sequência de valores começando em n  = 7

(sequência A061125 no OEIS )

Número de permutações que são desarranjos

Suponha que haja n pessoas em uma festa, cada uma trazendo um guarda-chuva. No final da festa, todos pegam um guarda-chuva da pilha de guarda-chuvas e folhas. Qual é a probabilidade de que ninguém tenha saído com seu guarda-chuva? Este problema é equivalente a contar permutações sem pontos fixos (chamadas de desarranjos ) e, portanto, o EGF, onde subtraímos pontos fixos (ciclos de comprimento 1) removendo o termo z da relação fundamental é

A multiplicação por soma dos coeficientes de , portanto , o número total de desarranjos, é dado por:

Portanto, há cerca de perturbações e a probabilidade de que uma permutação aleatória seja uma perturbação é

Este resultado também pode ser comprovado por inclusão-exclusão . Usando os conjuntos onde para designar o conjunto de permutações que fix p , temos

Esta fórmula conta o número de permutações que possuem pelo menos um ponto fixo. As cardinalidades são as seguintes:

Portanto, o número de permutações sem ponto fixo é

ou

e nós temos a reclamação.

Há uma generalização desses números, que é conhecida como números rencontres , ou seja, o número de permutações contendo m pontos fixos. O EGF correspondente é obtido marcando ciclos de tamanho um com a variável u , ou seja, escolhendo b ( k ) igual a um para e zero caso contrário, o que produz a função geradora do conjunto de permutações pelo número de pontos fixos:

Segue que

e, portanto

Isso imediatamente implica que

para n grande, m fixo.

Ordem de uma permutação aleatória

Se P é uma permutação, a ordem de P é o menor inteiro positivo n para o qual é a permutação de identidade. Este é o menor múltiplo comum dos comprimentos dos ciclos de P .

Um teorema de Goh e Schmutz afirma que se é a ordem esperada de uma permutação aleatória de tamanho n , então

onde a constante c é

Perturbações contendo um número par e ímpar de ciclos

Podemos usar a mesma construção da seção anterior para calcular o número de desarranjos contendo um número par de ciclos e o número contendo um número ímpar de ciclos. Para fazer isso, precisamos marcar todos os ciclos e subtrair pontos fixos, dando

Agora, um raciocínio muito básico mostra que o EGF de é dado por

Portanto, temos

qual é

Subtraindo de , encontramos

A diferença entre esses dois ( e ) é

Cem prisioneiros

Um diretor da prisão quer abrir espaço em sua prisão e está considerando a libertação de cem prisioneiros, libertando assim cem celas. Ele, portanto, reúne cem prisioneiros e pede-lhes que joguem o seguinte jogo: ele alinha cem urnas em uma fileira, cada uma contendo o nome de um prisioneiro, onde o nome de cada prisioneiro ocorre exatamente uma vez. O jogo é jogado da seguinte forma: cada prisioneiro tem permissão para olhar dentro de cinquenta urnas. Se ele ou ela não encontrar seu nome em uma das cinquenta urnas, todos os presos serão executados imediatamente, caso contrário o jogo continua. Os presos têm alguns momentos para decidir a estratégia, sabendo que uma vez iniciado o jogo, eles não poderão se comunicar, marcar as urnas de forma alguma ou movimentar as urnas ou os nomes dentro delas. Escolhendo urnas ao acaso, suas chances de sobrevivência são quase zero, mas existe uma estratégia que lhes dá 30% de chance de sobrevivência, supondo que os nomes sejam atribuídos às urnas aleatoriamente - o que é?

Em primeiro lugar, a probabilidade de sobrevivência usando escolhas aleatórias é

então esta definitivamente não é uma estratégia prática.

A estratégia de sobrevivência de 30% é considerar o conteúdo das urnas como uma permutação dos prisioneiros e ciclos transversais. Para manter a notação simples, atribua um número a cada prisioneiro, por exemplo, classificando seus nomes em ordem alfabética. Posteriormente, as urnas podem ser consideradas como contendo números em vez de nomes. Agora, claramente, o conteúdo das urnas define uma permutação. O primeiro prisioneiro abre a primeira urna. Se ele encontrar seu nome, ele terminou e sobreviveu. Caso contrário, ele abre a urna com o número que encontrou na primeira urna. O processo se repete: o preso abre uma urna e sobrevive se encontrar seu nome, caso contrário, abre a urna com o número que acabou de recuperar, até o limite de cinquenta urnas. O segundo prisioneiro começa com a urna número dois, o terceiro com a urna número três e assim por diante. Essa estratégia é precisamente equivalente a uma travessia dos ciclos da permutação representada pelas urnas. Cada prisioneiro começa com a urna com seu número e segue percorrendo seu ciclo até o limite de cinquenta urnas. O número da urna que contém seu número é a pré-imagem desse número sob a permutação. Conseqüentemente, os prisioneiros sobrevivem se todos os ciclos da permutação contiverem no máximo cinquenta elementos. Temos que mostrar que essa probabilidade é de pelo menos 30%.

Observe que isso pressupõe que o diretor escolhe a permutação aleatoriamente; se o diretor antecipar essa estratégia, ele pode simplesmente escolher uma permutação com um ciclo de 51. Para superar isso, os prisioneiros podem concordar previamente com uma permutação aleatória de seus nomes.

Consideramos o caso geral de presos e urnas sendo abertas. Primeiro calculamos a probabilidade complementar, ou seja, que existe um ciclo de mais de elementos. Com isso em mente, apresentamos

ou

de modo que a probabilidade desejada seja

porque o ciclo de mais do que elementos será necessariamente único. Usando o fato de que , descobrimos que

que produz

Finalmente, usando uma estimativa integral, como a soma de Euler-Maclaurin , ou a expansão assintótica do n- ésimo número harmônico , obtemos

de modo a

ou pelo menos 30%, conforme reivindicado.

Um resultado relacionado é que assintoticamente, o comprimento esperado do ciclo mais longo é λn, onde λ é a constante de Golomb-Dickman , aproximadamente 0,62.

Este exemplo é devido a Anna Gál e Peter Bro Miltersen; consulte o artigo de Peter Winkler para obter mais informações e veja a discussão em Les-Mathematiques.net . Consulte as referências sobre 100 prisioneiros para obter links para essas referências.

O cálculo acima pode ser executado de uma forma mais simples e direta, como segue: primeiro observe que uma permutação de elementos contém no máximo um ciclo de comprimento estritamente maior que . Assim, se denotarmos

então

Pois , o número de permutações que contêm um ciclo de comprimento exatamente é

Explicação: é a quantidade de maneiras de escolher os elementos que compõem o ciclo; é o número de maneiras de organizar os itens em um ciclo; e é o número de maneiras de permutar os elementos restantes. Não existe dupla contagem aqui porque não há no máximo um ciclo de comprimento quando . Desse modo,

Concluimos que

Uma variação do problema dos 100 prisioneiros (chaves e caixas)

Existe um problema intimamente relacionado que se ajusta perfeitamente ao método apresentado aqui. Digamos que você tenha n caixas encomendadas. Cada caixa contém uma chave para alguma outra caixa ou possivelmente ela mesma dando uma permutação das chaves. Você tem permissão para selecionar k dessas n caixas todas de uma vez e abri-las simultaneamente, obtendo acesso às teclas k . Qual é a probabilidade de que usando essas chaves você possa abrir todas as n caixas, onde você usa uma chave encontrada para abrir a caixa a que pertence e repetir.

O enunciado matemático desse problema é o seguinte: escolha uma permutação aleatória em n elementos e k valores na faixa de 1 a n , também aleatoriamente, chame essas marcas. Qual é a probabilidade de haver pelo menos uma marca em cada ciclo da permutação? A alegação é que essa probabilidade é k / n .

As espécies de permutações por ciclos com algum subconjunto não vazio de cada ciclo sendo marcado tem a especificação

O índice na soma interna começa em um porque devemos ter pelo menos uma marca em cada ciclo.

Traduzindo a especificação para funções geradoras, obtemos a função geradora bivariada

Isso simplifica para

ou

A fim de extrair coeficientes desta reescrita assim

Segue-se agora que

e, portanto

Divida por para obter

Não precisamos dividir por n! porque é exponencial em z .

Número de permutações contendo m ciclos

Aplicando o teorema fundamental de Flajolet-Sedgewick , ou seja, o teorema de enumeração rotulado com , ao conjunto

nós obtemos a função geradora

O termo

produz os números de Stirling assinados do primeiro tipo e é o EGF dos números de Stirling não assinados do primeiro tipo, ou seja,

Podemos calcular o OGF dos números de Stirling assinados para n fixos, ou seja,

Começar com

que produz

Somando isso, obtemos

Usando a fórmula envolvendo o logaritmo para à esquerda, a definição de à direita e o teorema binomial , obtemos

Comparando os coeficientes de , e usando a definição do coeficiente binomial , finalmente temos

um fatorial decrescente . O cálculo do OGF dos números de Stirling não assinados do primeiro tipo funciona de maneira semelhante.

Número esperado de ciclos de um determinado tamanho m

Neste problema, usamos uma função geradora bivariada g ( zu ) conforme descrito na introdução. O valor de b para um ciclo de tamanho m é zero e um para um ciclo de tamanho m . Nós temos

ou

Isso significa que o número esperado de ciclos de tamanho m em uma permutação de comprimento n menor que m é zero (obviamente). Uma permutação aleatória de comprimento de pelo menos m contém em média 1 / m ciclos de comprimento m . Em particular, uma permutação aleatória contém cerca de um ponto fixo.

O OGF do número esperado de ciclos de comprimento menor ou igual a m é, portanto,

onde H m é o m- ésimo número harmônico . Conseqüentemente, o número esperado de ciclos de comprimento no máximo m em uma permutação aleatória é cerca de ln  m .

Momentos de pontos fixos

O GF misto do conjunto de permutações pelo número de pontos fixos é

Seja a variável aleatória X o número de pontos fixos de uma permutação aleatória. Usando números de Stirling de segundo tipo , temos a seguinte fórmula para o m- ésimo momento de X :

onde está um fatorial decrescente . Usando , temos

que é zero quando , e um caso contrário. Portanto, apenas os termos com contribuem para a soma. Isso produz

Número esperado de pontos fixos em permutação aleatória elevado a alguma potência k

Suponha que você escolha uma permutação aleatória e a eleve a alguma potência , com um inteiro positivo, e pergunte sobre o número esperado de pontos fixos no resultado. Denote esse valor por .

Para cada divisor de um ciclo de comprimento se divide em pontos fixos quando elevado à potência. Portanto, precisamos marcar esses ciclos com. Para ilustrar isso, considere

Nós temos

qual é

Mais uma vez continuando conforme descrito na introdução, encontramos

qual é

A conclusão é que para e há quatro pontos fixos em média.

O procedimento geral é

Mais uma vez continuando como antes, encontramos

Mostramos que o valor de é igual a (o número de divisores de ) assim que começa em for e aumenta em um cada vez que atinge um divisor de até e incluindo a si mesmo.

Número esperado de ciclos de qualquer comprimento de uma permutação aleatória

Construímos a função geradora bivariada usando , onde é um para todos os ciclos (cada ciclo contribui com um para o número total de ciclos).

Observe que tem o formulário fechado

e gera os números de Stirling não assinados do primeiro tipo .

Nós temos

Conseqüentemente, o número esperado de ciclos é o número harmônico , ou cerca de .

Número de permutações com um ciclo de comprimento maior que n / 2

(Observe que a Seção Cem prisioneiros contém exatamente o mesmo problema com um cálculo muito semelhante, além de uma prova elementar mais simples.)

Mais uma vez, comece com a função geradora exponencial , desta vez da classe de permutações de acordo com o tamanho onde os ciclos de comprimento mais do que são marcados com a variável :

Só pode haver um ciclo de comprimento a mais , portanto, a resposta à pergunta é dada por

ou

qual é

O expoente de no termo sendo elevado à potência é maior do que e, portanto, nenhum valor para pode possivelmente contribuir para

Conclui-se que a resposta é

A soma tem uma representação alternativa que se encontra, por exemplo, no OEIS OEISA024167 .

finalmente dando

Número esperado de transposições de uma permutação aleatória

Podemos usar a decomposição de ciclo disjunto de uma permutação para fatorá-la como um produto de transposições, substituindo um ciclo de comprimento k por k  - 1 transposições. Por exemplo, os fatores de ciclo como . A função para ciclos é igual a e obtemos

e

Portanto, o número esperado de transposições é

onde está o número Harmônico . Também poderíamos ter obtido essa fórmula observando que o número de transposições é obtido somando os comprimentos de todos os ciclos (o que dá n ) e subtraindo um para cada ciclo (o que dá pela seção anterior).

Observe que novamente gera os números de Stirling sem sinal do primeiro tipo , mas na ordem inversa. Mais precisamente, temos

Para ver isso, observe que o acima é equivalente a

e essa

que vimos ser o EGF dos números de Stirling sem sinal do primeiro tipo na seção sobre permutações consistindo em precisamente m ciclos.

Tamanho de ciclo esperado de um elemento aleatório

Selecionamos um elemento aleatório q de uma permutação aleatória e perguntamos sobre o tamanho esperado do ciclo que contém q . Aqui, a função é igual a , porque um ciclo de comprimento k contribui com k elementos que estão em ciclos de comprimento k . Observe que, ao contrário dos cálculos anteriores, precisamos calcular a média desse parâmetro depois de extraí-lo da função geradora (dividir por n ). Nós temos

Portanto, a duração esperada do ciclo que contém q é

Probabilidade de que um elemento aleatório se encontre em um ciclo de tamanho m

Este parâmetro médio representa a probabilidade de que se selecionarmos novamente um elemento aleatório de de uma permutação aleatória, o elemento encontra-se em um ciclo de tamanho m . A função é igual a para e zero caso contrário, porque apenas os ciclos de comprimento m contribuem, ou seja, m elementos que estão em um ciclo de comprimento m . Nós temos

Segue-se que a probabilidade de que um elemento aleatório se encontre em um ciclo de comprimento m é

Probabilidade de que um subconjunto aleatório de [ n ] esteja no mesmo ciclo

Selecione um subconjunto aleatório Q de [ n ] contendo m elementos e uma permutação aleatória e pergunte sobre a probabilidade de que todos os elementos de Q estejam no mesmo ciclo. Este é outro parâmetro médio. A função b ( k ) é igual a , porque um ciclo de comprimento k contribui com subconjuntos de tamanho m , onde para k < m . Isso produz

Calculando a média, obtemos que a probabilidade dos elementos de Q estarem no mesmo ciclo é

ou

Em particular, a probabilidade de que dois elementos p < q estejam no mesmo ciclo é 1/2.

Número de permutações contendo um número par de ciclos pares

Podemos usar o teorema fundamental de Flajolet-Sedgewick diretamente e calcular estatísticas de permutação mais avançadas. (Verifique nessa página uma explicação de como os operadores que usaremos são calculados.) Por exemplo, o conjunto de permutações contendo um número par de ciclos pares é dado por

Traduzindo para funções geradoras exponenciais (EGFs), obtemos

ou

Isso simplifica para

ou

Isso diz que há uma permutação de tamanho zero contendo um número par de ciclos pares (a permutação vazia, que contém zero ciclos de comprimento par), uma tal permutação de tamanho um (o ponto fixo, que também contém zero ciclos de comprimento par ), e para isso , existem tais permutações.

Permutações que são quadrados

Considere o que acontece quando elevamos ao quadrado uma permutação. Os pontos fixos são mapeados para pontos fixos. Ciclos ímpares são mapeados para ciclos ímpares em uma correspondência um a um, por exemplo, se transforma em . Mesmo os ciclos se dividem em dois e produzem um par de ciclos com metade do tamanho do ciclo original, por exemplo, se transforma em . Portanto, as permutações que são quadrados podem conter qualquer número de ciclos ímpares e um número par de ciclos de tamanho dois, um número par de ciclos de tamanho quatro etc., e são dadas por

que produz o EGF

Invariantes de ciclo ímpar

Os tipos de permutações apresentadas nas duas seções anteriores, ou seja, permutações contendo um número par de ciclos pares e permutações que são quadrados, são exemplos dos chamados invariantes de ciclo ímpar , estudados por Sung e Zhang (ver links externos ). O termo invariante de ciclo ímpar significa simplesmente que a pertença à respectiva classe combinatória é independente do tamanho e do número de ciclos ímpares que ocorrem na permutação. Na verdade, podemos provar que todos os invariantes de ciclo ímpar obedecem a uma recorrência simples, que derivaremos. Primeiro, aqui estão mais alguns exemplos de invariantes de ciclo ímpar.

Permutações onde a soma dos comprimentos dos ciclos pares é de seis

Esta classe tem a especificação

e a função geradora

Os primeiros valores são

Permutações onde todos os ciclos pares têm a mesma duração

Esta classe tem a especificação

e a função geradora

Há uma nuance semântica aqui. Poderíamos considerar permutações que não contêm ciclos pares como pertencentes a essa classe, uma vez que zero é par . Os primeiros valores são

Permutações onde o comprimento máximo de um ciclo par é quatro

Esta classe tem a especificação

e a função geradora

Os primeiros valores são

A recorrência

Observe cuidadosamente como as especificações do componente de ciclo uniforme são construídas. É melhor pensar neles em termos de árvores de análise. Essas árvores têm três níveis. Os nós no nível mais baixo representam somas de produtos de ciclos de comprimento par do singleton . Os nós no nível intermediário representam restrições do operador de conjunto. Finalmente, o nó no nível superior soma os produtos das contribuições do nível médio. Observe que as restrições do operador de conjunto, quando aplicadas a uma função geradora que seja par, preservarão esse recurso, ou seja, produzirão outra função geradora par. Mas todas as entradas para os operadores de conjunto são uniformes, uma vez que surgem de ciclos de comprimento uniforme. O resultado é que todas as funções geradoras envolvidas têm a forma

onde está uma função par. Isso significa que

é igual também e, portanto,

Deixando e extraindo coeficientes, descobrimos que

que produz a recorrência

Um problema da competição Putnam de 2005

Um link para o site da competição Putnam aparece na seção Links externos . O problema pede uma prova de que

onde a soma é sobre todas as permutações de , é o sinal de , isto é , se é par e se é ímpar, e é o número de pontos fixos de .

Agora o sinal de é dado por

onde o produto está em todos os ciclos c de , como explicado, por exemplo, na página em permutações pares e ímpares .

Portanto, consideramos a classe combinatória

onde marca um menos a duração de um ciclo de contribuição e marca pontos fixos. Traduzindo para funções geradoras, obtemos

ou

Agora temos

e, portanto, a quantidade desejada é dada por

Fazendo o cálculo, obtemos

ou

Extraindo os coeficientes, descobrimos que o coeficiente de é zero. A constante é um, o que não concorda com a fórmula (deveria ser zero). Para positivo, no entanto, obtemos

ou

qual é o resultado desejado.

Como um aparte interessante, observamos que pode ser usado para avaliar o seguinte determinante de uma matriz:

onde . Lembre-se da fórmula para o determinante:

Agora, o valor do produto à direita para uma permutação é , onde f é o número de pontos fixos de . Por isso

que produz

e finalmente

A diferença entre o número de ciclos em permutações pares e ímpares

Aqui procuramos mostrar que essa diferença é dada por

Lembre-se de que o sinal de uma permutação é dado por

onde o produto varia ao longo dos ciclos c da composição do ciclo disjunto de .

Segue-se que a espécie combinatória que reflete os sinais e a contagem do ciclo do conjunto de permutações é dada por

onde usamos para marcar sinais e para a contagem do ciclo.

Traduzindo para funções geradoras que temos

Isso simplifica para

qual é

Agora, as duas funções geradoras e de permutações pares e ímpares por contagem de ciclo são dadas por

e

Exigimos a quantidade

qual é

Finalmente, extraindo coeficientes desta função geradora, obtemos

qual é

que é por sua vez

Isso conclui a prova.

Veja também

Referências

links externos

100 prisioneiros