Estimating $\pi$ with a Coin

Este artigo descreve um novo método de Monte Carlo para estimar o valor de $\pi$ através do lançamento de moedas, fundamentado em identidades de séries de números de Catalan que, embora presentes implicitamente na literatura de teoria das probabilidades, ganham aqui uma interpretação inédita para $\frac{\pi}{4}$.

O essencial

Imagine que você tem um objetivo muito difícil: descobrir o valor exato de π (o famoso número 3,14159...), mas você não pode usar régua, compasso ou calculadora. Você só tem uma moeda e um pouco de paciência.

O artigo de Jim Propp descreve um método divertido e um pouco "maluco" para fazer exatamente isso, usando apenas o acaso de jogar moedas. Vamos explicar como funciona, passo a passo, com analogias do dia a dia.

1. A Regra do Jogo: A Corrida até a Vitória

Imagine que você e um amigo estão numa corrida. Você é "Cara" (Heads) e seu amigo é "Coroa" (Tails).

• Vocês começam empatados (0 a 0).
• Cada vez que a moeda cai, um ponto vai para um lado.
• A regra é: Você deve continuar jogando até que você (Cara) fique à frente do seu amigo (Coroa) pela primeira vez.

Assim que você ganha a liderança (ex: 1 a 0, ou 2 a 1, ou 3 a 2), você para imediatamente.

2. O Truque Matemático

Agora vem a parte mágica. Quando você para, você olha para o placar final e calcula uma fração:

Quantos "Caras" você teve dividido pelo Total de jogadas.

• Exemplo 1: A moeda caiu "Coroa", depois "Cara". Você parou.

  • Placar: 1 Cara, 1 Coroa? Não, espere! Se caiu Coroa primeiro, você ainda está perdendo (0 a 1). Você tem que continuar.
  • Vamos tentar de novo: Coroa, Cara, Cara.
  • Placar: 2 Caras, 1 Coroa. Você venceu!
  • A fração é: 2 dividido por 3 (2/3). Você anota esse número.

• Exemplo 2: Coroa, Coroa, Cara, Coroa, Cara, Cara.

  • Placar: 3 Caras, 3 Coroa? Não, no último passo você venceu.
  • Placar final: 3 Caras, 3 Coroa? Não, espere. Vamos contar: C, C, Ca, C, Ca, Ca.
  • No passo 1: C (0-1). Passo 2: C (0-2). Passo 3: Ca (1-2). Passo 4: C (1-3). Passo 5: Ca (2-3). Passo 6: Ca (3-3). Ainda empatado? Não, a regra é "exceder".
  • Vamos simplificar: O jogo para no momento exato em que o número de Caras é maior que o de Coroas.
  • Se o jogo parou com 3 Caras e 2 Coroas, a fração é 3/5 (0,6).

3. O Segredo: Repetir e Média

Se você fizer isso apenas uma vez, o resultado será aleatório. Pode ser 0,6, 0,55, 0,7... nada que pareça com 3,14.

Mas, se você repetir esse jogo milhares de vezes (ou pedir para mil pessoas jogarem ao mesmo tempo), e depois calcular a média de todas essas frações que você anotou, algo incrível acontece:

O resultado da média vai se aproximar magicamente de π/4 (ou seja, 3,14159... dividido por 4, que é cerca de 0,785).

Para descobrir o valor de π, basta multiplicar esse resultado por 4.

4. Por que isso funciona? (A Analogia da Floresta)

O autor explica isso usando a teoria de "caminhadas aleatórias". Imagine que você está numa floresta e dá passos para a direita (Cara) ou para a esquerda (Coroa).

• A matemática por trás disso envolve números chamados Números de Catalan (que aparecem em muitos problemas de contagem, como formas de desenhar parênteses ou caminhos em grades).
• Existe uma conexão profunda entre esses caminhos que nunca cruzam uma linha e a função matemática chamada Arcosseno.
• Quando o autor soma todas as probabilidades de todos os caminhos possíveis onde você vence pela primeira vez, a matemática "se encaixa" perfeitamente na fórmula do Arcosseno de 1, que é igual a π/2. Dividindo por 2 (devido à forma como a fração é calculada), chegamos a π/4.

É como se o universo tivesse escondido o número π dentro da estrutura básica do acaso, e essa brincadeira de moeda é a chave para abri-lo.

5. A Realidade: É Prático?

Aqui vem o "choque de realidade" do artigo.

• A precisão é lenta: Para obter um resultado bom, você precisa de muitas jogadas. O erro diminui muito devagar.
• O exemplo do Matt Parker: Um youtuber famoso de matemática tentou isso com 10.000 jogadas e obteve 3,22. O valor real é 3,14. O erro foi grande.
• O tempo necessário: Para chegar perto de 3,14159 com precisão, o autor diz que você precisaria de um trilhão de jogadas. Se você jogasse uma moeda por segundo, sem parar, levaria mais de 30.000 anos para terminar o experimento!

Conclusão

Este artigo não é um método para calcular π na prática (é muito lento e ineficiente). É, na verdade, uma prova matemática elegante e divertida.

Ele mostra que, mesmo em algo tão simples quanto jogar uma moeda, a estrutura do universo esconde o número π. É uma maneira de "ver" a matemática pura acontecendo no mundo real, provando que π está entre 3 e 4 apenas olhando para a média de vitórias em uma corrida de moedas.

Resumo em uma frase: Se você jogar moedas até ganhar pela primeira vez, anotar sua porcentagem de vitórias e repetir isso infinitas vezes, a média desses resultados revelará o valor de π, escondido no coração do acaso.

Resumo técnico

Aqui está um resumo técnico detalhado do artigo "Estimating π with a Coin" (Estimando π com uma Moeda), de Jim Propp, apresentado em português:

Resumo Técnico: Estimativa de π através de Caminhadas Aleatórias e Paradas Otimizadas

1. O Problema

O artigo aborda a questão clássica de estimar o valor de $\pi$ utilizando métodos probabilísticos (Monte Carlo). Embora o método de agulha de Buffon seja bem conhecido, o autor propõe uma abordagem alternativa baseada em lançamentos de moedas. O objetivo é demonstrar que, ao lançar uma moeda justa repetidamente até que o número de "caras" (Heads) exceda pela primeira vez o número de "coroas" (Tails), o valor esperado da proporção de caras no momento da parada converge para $\pi/4$.

2. Metodologia

A prova baseia-se na teoria de caminhadas aleatórias simples (random walks) e tempos de parada (stopping times).

• Definição do Processo:

  • Seja $H_n$ e $T_n$ o número de caras e coroas nos primeiros $n$ lançamentos.
  • Define-se o passeio aleatório $S_n = H_n - T_n$, onde $S_0 = 0$.
  • O tempo de parada $\tau$ é definido como o primeiro instante em que o passeio atinge $+1$:
    $$\tau = \min\{n \ge 1 : S_n > 0\}$$
  • Como o passeio muda em $\pm 1$ a cada passo, $\tau$ assume apenas valores ímpares ($2k+1$).

• Análise Combinatória:

  • Para que $\tau = 2k+1$, o passeio deve permanecer $\le 0$ nos primeiros $2k$passos, retornar a 0 no tempo$2k$e, em seguida, dar um passo para$+1$.
  • O número de tais trajetórias é dado pelo número de Catalan $C_k = \frac{1}{k+1}\binom{2k}{k}$.
  • A probabilidade de $\tau = 2k+1$ é:
    $$P(\tau = 2k+1) = \frac{1}{2} \cdot \frac{1}{4^k} \cdot \frac{1}{k+1} \binom{2k}{k}$$

• Cálculo do Valor Esperado:

  • No momento da parada $\tau = 2k+1$, temos $H_\tau = k+1$ e $T_\tau = k$.
  • A proporção de caras é $H_\tau/\tau = \frac{k+1}{2k+1}$.
  • O valor esperado $E[H_\tau/\tau]$ é calculado somando sobre todos os $k \ge 0$:
    $$E\left[\frac{H_\tau}{\tau}\right] = \sum_{k \ge 0} P(\tau = 2k+1) \cdot \frac{k+1}{2k+1}$$
  • Substituindo a probabilidade e simplificando, a série resultante é:
    $$\frac{1}{2} \sum_{k \ge 0} \frac{1}{4^k} \frac{1}{2k+1} \binom{2k}{k}$$

• Conexão com Séries de Potência:

  • O autor identifica que a série acima corresponde à expansão em série de Taylor da função arco-seno ($\arcsin(x)$) avaliada em $x=1$:
    $$\arcsin(x) = \sum_{n=0}^{\infty} \frac{1}{2^{2n}} \binom{2n}{n} \frac{x^{2n+1}}{2n+1}$$
  • Como $\arcsin(1) = \pi/2$, o valor esperado torna-se:
    $$E\left[\frac{H_\tau}{\tau}\right] = \frac{1}{2} \cdot \frac{\pi}{2} = \frac{\pi}{4}$$

3. Contribuições Chave

• Nova Interpretação Probabilística: Embora a identidade envolvendo a série de Catalan e $\pi$ fosse conhecida implicitamente na literatura de teoria da probabilidade (relacionada a $E[1/\tau] = \pi/2 - 1$), a interpretação de $\pi/4$ como o valor esperado da proporção de caras no tempo de parada é apresentada como uma contribuição nova.
• Prova Probabilística de Limites de $\pi$: O artigo observa que, como a proporção $H_\tau/\tau$ é sempre maior que $1/2$(e igual a 1 em metade dos casos), o valor esperado deve estar estritamente entre$3/4$e$1$. Isso fornece uma prova probabilística elegante de que$3 < \pi < 4$.
• Generalização para Outros Excedentes: O autor explora brevemente o caso onde a moeda é lançada até que as caras excedam as coroas por um valor $a$.
  • Para $a=2$, o valor esperado é $\ln(2)$.
  • Para $a$ ímpar, surgem expressões envolvendo $\pi$; para $a$ par, envolvem $\ln(2)$. (Nota: O autor menciona que essas generalizações foram verificadas via IA e permanecem como conjecturas não totalmente validadas no texto).

4. Resultados e Eficiência

• Simulação Prática: O artigo relata um teste realizado por Matt Parker com 10.000 lançamentos, resultando em uma estimativa de $\pi \approx 3.2266$.
• Convergência Lenta: A precisão do método é limitada. O erro diminui na ordem de $1/N^{1/4}$, onde$N$ é o número total de lançamentos.
  • Para obter uma precisão próxima de 3.14, seriam necessários cerca de um trilhão de lançamentos, o que levaria mais de 30.000 anos se feito a uma taxa de um lançamento por segundo.
• Comparação: O autor reconhece que algoritmos baseados em moedas explorados por Flajolet, Pelletier e Soria são muito mais eficientes que este método específico.

5. Significado

O trabalho destaca a beleza da interseção entre combinatória (números de Catalan), análise (séries de arco-seno) e teoria da probabilidade (caminhadas aleatórias). Embora não seja um método prático para calcular $\pi$ devido à sua lentidão de convergência, serve como um exemplo pedagógico poderoso e uma prova matemática rigorosa de como constantes fundamentais podem emergir de processos estocásticos simples. Além disso, oferece uma nova perspectiva sobre a relação entre o tempo de parada de um passeio aleatório e a constante $\pi$.