TOPICS
Search

Göbel's Sequence


Göbel's sequence is defined by the recurrence relation

 x_n=(x_0^2+x_1^2+...+x_(n-1)^2)/(n-1)
(1)

for n>=2, with x_0=x_1=1. This indexing gives x_(n+1)=g_(2,2)(n) for n>=1 in the notation of Gima et al. (2024). The first few terms of x_n for n=0, 1, ... are 1, 1, 2, 3, 5, 10, 28, 154, 3520, ... (OEIS A003504). The terms grow extremely rapidly, but are given by the asymptotic formula

 x_(n+1)∼(n+2-n^(-1)+4n^(-2)-21n^(-3)+138n^(-4)
 -1091n^(-5)+...)C^(2^n)
(2)

(OEIS A116603; correcting Finch 2003, p. 446), where

 C=1.04783144757641122955990946274313755459...
(3)

(OEIS A115632; Finch 2003, p. 446; Zagier).

The terms x_n are integers up to x_(43), but x_(44) is a nonintegral rational number (Gima et al. 2024). The numbers of decimal digits in the numerators of x_n for n=0, 1, ... are 1, 1, 1, 1, 1, 2, 2, 3, 4, 7, 12, 22, 43, 85, ... (OEIS A389769). While an elegant approach using congruences can show that x_(44) is not an integer (see below), the exact value x_(44)=N/43 was found explicitly by E. Weisstein (Oct. 13, 2025), following smaller initial computations by S. Wagon, in a calculation requiring 370.9 gigabytes of memory. The exact value of N is given by N=23260...77568, where N has 178485291570 decimal digits which require 74.1 gigabytes to represent internally.

It is more convenient to work with the transformed sequence

 s_n=2+sum_(j=2)^nx_j^2=nx_(n+1),
(4)

which gives the new recurrence relation

 s_(n+1)=s_n+(s_n^2)/(n^2)=s_n+((s_n)/n)^2
(5)

for n>=1, with initial term s_1=2. The first few terms for n=1, 2, ... are 2, 6, 15, 40, 140, 924, 24640, ... (OEIS A061322). If s_n is an integer, then s_(n+1) is nonintegral iff ns_n. For a prime number p, a nonzero residue of s_p modulo p shows that x_(p+1)=s_p/p is nonintegral. Checking only prime numbers does not by itself rule out an earlier failure at a composite index; prime power divisibility must also be checked (Zagier).

For example, the first few sequences {s_n (mod p)}_(n=1)^p are summarized in the following table. Residues of rational numbers whose denominators are not divisible by p are computed using the multiplicative inverse of the denominator modulo p.

p{s_n (mod p)}_(n=1)^p
20, 0
32, 0, 0
52, 1, 0, 0, 0
72, 6, 1, 5, 0, 0, 0
112, 6, 4, 7, 8, 0, 0, 0, 0, 0, 0
132, 6, 2, 1, 10, 1, 5, 1, 0, 0, 0, 0, 0
172, 6, 15, 6, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 0

While s_n (and therefore their values modulo p) can be computed explicitly for small values of n and p, the fractions quickly grow too large to represent exactly. However, computing terms modulo p directly avoids term growth, and the residue for p=43 proves that x_(44) is nonintegral. As expected, direct verification of this fact is difficult since

 x_(44) approx 5.4093×10^(178485291567)
(6)

(calculated using the asymptotic formula) is very large to be computed and stored explicitly. The first few values of p for which s_p≢0  (modp) are 43, 61, 67, 83, 103, 107, 109, 157, ... (OEIS A378851).

A sequence even more striking for assuming integer values only for finitely many terms is the 3-Göbel sequence

 g_n=(1+g_0^3+g_1^3+...+g_(n-1)^3)/n
(7)

for n>=1, with g_0=1. The first few terms are 1, 2, 5, 45, 22815, ... (OEIS A005166).

The Göbel sequences can be generalized to integer powers k>=2 by

 g_n=(1+g_0^k+g_1^k+...+g_(n-1)^k)/n
(8)

for n>=1, with g_0=1.

Göbel's sequences can be generalized to (k,l)-Göbel sequences, defined by the recursion

 (n+1)g_(k,l)(n+1)=g_(k,l)(n)(n+g_(k,l)(n)^(k-1))
(9)

for integers k,l>=2 and with initial value g_(k,l)(1)=l (Ibstedt 1990, Gima et al. 2024).


See also

Somos Sequence

Explore with Wolfram|Alpha

References

Finch, S. R. Mathematical Constants. Cambridge, England: Cambridge University Press, p. 446, 2003.Gima, H.; Matsusaka, T.; Miyazaki, T.; and Yara, S. "On Integrality and Asymptotic Behavior of the (k,l)-Göbel Sequences." 14 Feb 2024. https://arxiv.org/abs/2402.09064.Guy, R. K. "The Strong Law of Small Numbers." Amer. Math. Monthly 95, 697-712, 1988.Guy, R. K. "A Recursion of Göbel." §E15 in Unsolved Problems in Number Theory, 2nd ed. New York: Springer-Verlag, pp. 214-215, 1994.Ibstedt, H. "Some Sequences of Large Integers." Fib. Quart. 28, No. 3, 200-203, 1990.Sloane, N. J. A. Sequences A003504/M0728, A005166/M1551, A061322, A115632, A116603, A378851, and A389769 in "The On-Line Encyclopedia of Integer Sequences."Stone, A. "The Astonishing Behavior of Recursive Sequences." Quanta. Nov. 16, 2023. https://www.quantamagazine.org/the-astonishing-behavior-of-recursive-sequences-20231116/.Zagier, D. "Solution: Day 5, Problem 3." https://mathshistory.st-andrews.ac.uk/EMS/Zagier/Solution53/.

Referenced on Wolfram|Alpha

Göbel's Sequence

Cite this as:

Weisstein, Eric W. "Göbel's Sequence." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/GoebelsSequence.html

Subject classifications