Hofstadter G-sequence

(defun G (n)
"Emacs Lisp function
for OEIS sequence number A005206.
Hofstadter G-sequence: a(n) = n - a(a(n-1)),
where
a(n)=floor((n+1)*tau)-n-1 where tau=(1+sqrt(5))/2;
a(0)=0, a(1)=1, a(n)=n-a(floor(n/tau))
or
phi=(1+sqrt(5))/2
G(n)=(n=0 -> 0, n=1 -> 1, t -> n-G([n/phi]))"
  (let ((phi (/ (+ 1 (sqrt 5)) 2)))
    (cond ((eq n 0) 0)
          ((eq n 1) 1)
          (t (- n (G (floor (/ n phi))))))))


(defun table-Hofstadter-G-sequence (n)
"Table of k, G(k) for k=0...n"
  (let ((k 0))
    (while (<= k n)
      (princ k )
      (princ " ")
      (princ (G k))
      (princ "\n")
      (setq k (+  1 k)))))

A tabela dos 10 primeiros termos da sequência de Hofstadter obtida com a instrução é (table-Hofstadter-G-sequence 10) é

Um cálculo mais exaustivo com 10000 termos está aqui.

Refs:

http://www.research.att.com/~njas/sequences/A005206

Palavras chave/keywords: Hofstadter G-sequence, OEIS A005206, Lisp, Emacs

Criado/Created: NaN

Última actualização/Last updated: 10-10-2022 [14:47]

Voltar à página inicial.