Matemático e IA provam existência de hexágonos mágicos de qualquer ordem

O desenvolvedor e pesquisador Konstantin Gukov publicou a solução para um problema aberto na teoria dos números e na geometria combinatória: a prova de que existem hexágonos mágicos anômalos para qualquer ordem n > 3. A descoberta, desenvolvida com o apoio do modelo de inteligência artificial GPT-5.6 Sol e do agente de prova formal Aristotle, supera o limite anterior estabelecido em 2024 pelo pesquisador Klaus Meffert, que havia encontrado computacionalmente uma solução isolada para a ordem n = 9.

A investigação começou durante as celebrações de 19 anos da YSDA (Yandex School of Data Analysis), quando ex-alunos discutiram as propriedades do número 19. Além de ser um primo gêmeo, 19 é a quantidade exata de células que compõem o único hexágono mágico normal não trivial possível na matemática, cuja estrutura de ordem n = 3 foi demonstrada no passado como única devido a restrições de divisibilidade impostas pelas linhas da grelha.

Para contornar a impossibilidade teórica de criar hexágonos normais maiores, Konstantin Gukov concentrou seu trabalho nos hexágonos mágicos anômalos. Nesses arranjos, os números contidos nas células permanecem obrigatoriamente consecutivos, mas não precisam iniciar do valor 1. Utilizando uma abordagem que combinou simetria, campos de potencial e algoritmos de simulated annealing otimizados com Numba em 24 núcleos de CPU, o projeto obteve primeiro soluções numéricas até a ordem n = 21 e, em seguida, uma demonstração construtiva completa para ordens arbitrárias, gerando exemplos de validação até a ordem n = 500.

O problema dos hexágonos mágicos

Diferentemente dos quadrados mágicos — grelhas numéricas conhecidas há milênios para as quais existem algoritmos determinísticos de construção normal para qualquer ordem n > 2 —, os arranjos hexagonais impõem severas limitações geométricas. Em um quadrado mágico normal de ordem n, os inteiros consecutivos variam de 1 a n², e a soma de cada linha, coluna e duas diagonais principais resulta em um mesmo valor constante.

Em uma grelha hexagonal de ordem n, as células alinham-se em três direções retas distintas no plano. A quantidade total de células no hexágono é dada pela fórmula 3n² - 3n + 1, dispostas em 2n - 1 linhas retas ao longo de cada um dos três eixos. Para que a estrutura seja considerada mágica, a soma de todos os números contidos em qualquer uma dessas linhas deve ser exatamente idêntica.

A prova de que não existem hexágonos mágicos normais para ordens superiores a n = 3 baseia-se em um teste de divisibilidade simples. Como as células são particionadas em 2n - 1 linhas em cada uma das três direções, a soma total da sequência de números de 1 a 3n² - 3n + 1 precisa ser divisível por 2n - 1. Para a ordem n = 3 (que possui 19 células e 5 linhas por direção), a divisão resulta em um número inteiro; contudo, para qualquer ordem n > 3, o somatório total falha sistematicamente no teste de divisibilidade.

A única alternativa para buscar estruturas hexagonais de maior escala foi relaxar o ponto de partida da sequência numérica. Ao permitir que os números fossem consecutivos mas começassem em valores arbitrários (positivos ou negativos), a comunidade matemática passou a investigar os hexágonos anômalos. No entanto, até julho de 2026, como constava nos registros da Wikipedia, o maior arranjo desse tipo conhecido era o hexágono de ordem n = 9 descoberto por Klaus Meffert em 2024 por meio de busca computacional intensiva.

Redução do espaço por antissimetria

Ao analisar a complexidade do problema, Konstantin Gukov observou que a busca direta por força bruta em arranjos genéricos tornava-se inviável devido ao crescimento exponencial do espaço de estados. Em vez de recorrer a solucionadores genéricos de restrições como Z3 ou OR-Tools, o autor decidiu reduzir o domínio de busca aplicando restrições de antissimetria sobre a grelha.

A primeira modificação consistiu em definir que os números do hexágono pertenceriam ao intervalo simétrico -K, ..., K, fixando o número 0 no centro exato da estrutura. A partir desse ponto, estabeleceu-se a regra de antissimetria: duas células opostas entre si por uma rotação de 180 graus (células antipodais) devem conter valores opostos. Assim, se uma célula armazena o valor x, sua célula antipodal armazena automaticamente -x.

A imposição dessa antissimetria eliminou grande parte das restrições do sistema de equações. Qualquer linha reta que passa pelo centro do hexágono tem sua soma reduzida a 0 de forma automática, pois seus pares de células antipodais cancelam-se mutuamente. Para as linhas retas que não cruzam o centro, a soma dos elementos de sua linha antipodal correspondente é o inverso exato, o que significa que, se uma linha somar zero, sua correspondente também somará zero.

Embora a inclusão dessa simetria rotacional corresse o risco de inviabilizar a existência de soluções, ela permitiu restringir as buscas exclusivamente ao espaço de soluções onde a constante mágica de todas as linhas retas é exatamente igual a 0, simplificando significativamente a complexidade computacional para as etapas seguintes.

Representação por campo de potencial

Durante a análise dos hexágonos de soma zero, Konstantin Gukov identificou uma estrutura algébrica subjacente fundamentada em anéis alternados de 6 pontos. Em qualquer grelha hexagonal, ao tomar as seis células que cercam um determinado ponto interior e somar a elas o padrão alternado [-1, +1, -1, +1, -1, +1] — mantendo a célula central inalterada —, a soma de todas as linhas retas que cruzam esse anel permanece inalterada.

Esse fenômeno ocorre porque qualquer linha reta que intercepta um anel local de 6 pontos recebe exatamente duas contribuições opostas (um valor +1 e um valor -1) ou nenhuma contribuição. Como resultado, é possível adicionar qualquer múltiplo desse padrão alternado sobre a grelha sem alterar o somatório de nenhuma linha reta do hexágono.

Esses anéis locais formam uma base linear para o sistema. Dessa forma, qualquer hexágono de soma zero de ordem n pode ser representado de forma equivalente por um campo de potencial de ordem n - 1, que registra os coeficientes de cada anel local aplicado. Por construção, qualquer configuração definida dentro desse campo de potencial satisfaz automaticamente todas as equações de soma de linha igual a zero.

Com a representação por campo de potencial, o problema computacional mudou de figura: a exigência de que as linhas somem zero passou a ser garantida pela própria estrutura do campo. O novo desafio consistiu inteiramente em encontrar configurações no campo de potencial que resultassem em valores visíveis de células que fossem estritamente inteiros, distintos e consecutivos.

Desenvolvimento do resolvedor com IA

Para construir um algoritmo capaz de navegar no espaço dos campos de potencial, Konstantin Gukov utilizou técnicas experimentadas durante sua participação no processo de preparação do Midnight Code Cup 2026, uma competição de programação que incentiva explicitamente o uso de modelos de linguagem (LLMs) para o desenvolvimento de otimizadores específicos de domínio.

O modelo GPT-5.6 Sol foi utilizado para formular a arquitetura do solucionador. Ao analisar o problema dos hexágonos e a representação por campos de potencial, o modelo conectou a estrutura matemática ao conceito de arranjos de Heffter (Heffter arrays), que são configurações combinatórias de inteiros com sinal e condições de soma nula. Essa conexão forneceu heurísticas para a permuta de valores entre células mantendo o controle das somas globais.

Em vez de utilizar solucionadores formais genéricos como Z3 ou OR-Tools, o programa gerado pelo GPT-5.6 Sol implementou um algoritmo customizado de simulated annealing (têmpera simulada). Para garantir alta eficiência de processamento, o código passou por otimização técnica utilizando a biblioteca Numba para compilar os laços críticos de execução diretamente em código de máquina.

Com o auxílio da ferramenta de perfilamento perf em ambiente Linux, o autor identificou gargalos de alocação de memória e geração de números aleatórios, obtendo um ganho adicional de 50% de desempenho. A aplicação foi colocada em execução em um servidor doméstico utilizando 24 núcleos de CPU ao longo de vários dias.

O sistema encontrou novas soluções para hexágonos mágicos anômalos de ordens inéditas até n = 21. Ao projetar visualmente os campos de potencial correspondentes a essas soluções de ordem n = 21, observou-se que, embora os números contidos nas células parecessem distribuídos de forma caótica, seus campos de potencial apresentavam superfícies contínuas e suaves, semelhantes a mapas topográficos com elevações, vales e cristas.

Prova matemática e agente Aristotle

A obtenção recorrente de soluções numéricas até a ordem n = 21 permitiu formular a conjectura de que existem hexágonos mágicos anômalos para todas as ordens n > 3. Para desenvolver a demonstração matemática formal dessa hipótese, Konstantin Gukov combinou a capacidade do modelo GPT-5.6 Sol (em suas modalidades high e max) com o agente Aristotle, especializado na linguagem de prova formal Lean.

Nas primeiras etapas, o GPT-5.6 Sol (high) decompôs a hipótese em lemas menores e propostas de construções combinatórias, enquanto o agente Aristotle tentava codificar os passos intermediários em Lean. O processo travou temporariamente após alguns dias, quando os modelos passaram a repetir argumentos conhecidos sem avançar na redução das cotas teóricas.

O autor acionou então o GPT-5.6 Sol (max) em uma sessão de raciocínio contínuo de várias horas. Embora o modelo não tenha concluído a prova inteira antes do término dos créditos de processamento, ele produziu abordagens conceituais inéditas que foram incorporadas ao contexto de trabalho das sessões seguintes.

Reutilizando esses insights no GPT-5.6 Sol (high), a equipe obteve a primeira prova construtiva para ordens n > 800 que fossem divisíveis por 16. A partir dessa conquista, o processo de generalização reduziu progressivamente os critérios de divisibilidade: primeiro para ordens divisíveis por 8, depois por 4, por 2 e, finalmente, eliminou qualquer exigência de divisibilidade.

O limite inferior da demonstração teórica caiu de n = 800 para n = 114, patamar no qual as desigualdades e seleções combinatórias da construção determinística puderam ser justificadas do ponto de vista matemático. O resultado final produziu um algoritmo determinístico que gera hexágonos mágicos para qualquer ordem elevada. Unindo as testemunhas encontradas por busca direta até n = 21 com a construção determinística válida acima desse patamar, a existência de hexágonos mágicos anômalos ficou demonstrada para todo n > 3, permitindo a geração de hexágonos de ordem n = 500 de forma instantânea.

Validação formal e impacto técnico

Apesar da conclusão da prova construtiva e do algoritmo determinístico, Konstantin Gukov ressaltou que a demonstração matemática generalizada ainda não havia sido inteiramente formalizada na linguagem Lean até a publicação do trabalho, nem submetida a processo formal de revisão por pares (peer review). A verificação completa dentro do ambiente Lean com apoio do agente Aristotle consta como o próximo passo do projeto.

A resolução do problema traz implicações estratégicas para o desenvolvimento de softwares de otimização combinatória. Os resultados observados na Midnight Code Cup 2026 e no projeto do hexágono sugerem que modelos de linguagem avançados como o GPT-5.6 Sol são capazes de superar solucionadores genéricos de restrições (como Z3 ou OR-Tools) ao projetar algoritmos de simulated annealing altamente customizados para a geometria do problema.

O trabalho exemplifica a evolução no papel do desenvolvedor de tecnologia, que passa da codificação manual de laços otimizados em Numba para a condução e orquestração de sistemas autônomos de inteligência artificial aplicados à resolução de problemas científicos complexos.