cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A215409 The Goodstein sequence G_n(3).

This page as a plain text file.
%I A215409 #59 Feb 16 2025 08:33:18
%S A215409 3,3,3,2,1,0
%N A215409 The Goodstein sequence G_n(3).
%C A215409 G_0(m) = m. To get the 2nd term, write m in hereditary base 2 notation (see links), change all the 2s to 3s, and then subtract 1 from the result. To get the 3rd term, write the 2nd term in hereditary base 3 notation, change all 3s to 4s, and subtract 1 again. Continue until the result is zero (by Goodstein's Theorem), when the sequence terminates.
%C A215409 Decimal expansion of 33321/100000. - _Natan Arie Consigli_, Jan 23 2015
%H A215409 R. L. Goodstein, <a href="http://www.jstor.org/stable/2268019">On the Restricted Ordinal Theorem</a>, J. Symb. Logic 9, 33-41, 1944.
%H A215409 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/HereditaryRepresentation.html">Hereditary Representation</a>
%H A215409 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/GoodsteinSequence.html">Goodstein Sequence</a>
%H A215409 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/GoodsteinsTheorem.html">Goodstein's Theorem</a>
%H A215409 Wikipedia, <a href="http://en.wikipedia.org/wiki/Goodstein&#39;s_theorem#Hereditary_base-n_notation">Hereditary base-n notation</a>
%H A215409 Wikipedia, <a href="http://en.wikipedia.org/wiki/Goodstein&#39;s_theorem#Goodstein_sequences">Goodstein sequence</a>
%H A215409 Wikipedia, <a href="http://en.wikipedia.org/wiki/Goodstein&#39;s_theorem">Goodstein's Theorem</a>
%H A215409 Reinhard Zumkeller, <a href="/A211378/a211378.hs.txt">Haskell programs for Goodstein sequences</a>
%F A215409 a(0) = a(1) = a(2) = 3; a(3) = 2; a(4) = 1; a(n) = 0, n > 4;
%F A215409 From _Iain Fox_, Dec 12 2017: (Start)
%F A215409 G.f.: 3 + 3*x + 3*x^2 + 2*x^3 + x^4.
%F A215409 E.g.f.: 3 + 3*x + (3/2)*x^2 + (1/3)*x^3 + (1/24)*x^4.
%F A215409 a(n) = floor(2 - (4/Pi)*arctan(n-3)), n >= 0.
%F A215409 (End)
%e A215409 a(0) = 3 = 2^1 + 1;
%e A215409 a(1) = 3^1 + 1 - 1 = 3^1 = 3;
%e A215409 a(2) = 4^1 - 1 = 3;
%e A215409 a(3) = 3 - 1 = 2;
%e A215409 a(4) = 2 - 1 = 1;
%e A215409 a(5) = 1 - 1 = 0.
%t A215409 PadRight[CoefficientList[Series[3 + 3 x + 3 x^2 + 2 x^3 + x^4, {x, 0, 4}], x], 6] (* _Michael De Vlieger_, Dec 12 2017 *)
%o A215409 (Haskell) -- See Link
%o A215409 (PARI) B(n, b)=sum(i=1, #n=digits(n, b), n[i]*(b+1)^if(#n<b+i, #n-i, B(#n-i, b)))
%o A215409 a(n) = my(x=3); for(i=1, n, x=B(x, i+1)-1; if(x==0, break())); x \\ (uses definition of sequence) _Iain Fox_, Dec 13 2017
%o A215409 (PARI) first(n) = my(res = vector(n)); res[1] = res[2] = res[3] = 3; res[4] = 2; res[5] = 1; res; \\ _Iain Fox_, Dec 12 2017
%o A215409 (PARI) first(n) = Vec(3 + 3*x + 3*x^2 + 2*x^3 + x^4 + O(x^n)) \\ _Iain Fox_, Dec 12 2017
%o A215409 (PARI) a(n) = floor(2 - (4/Pi)*atan(n-3)) \\ _Iain Fox_, Dec 12 2017
%Y A215409 Cf. A056004, A057650, A059934, A059935, A059936, A271977.
%Y A215409 Cf. A056041, A056193, A266204, A271554, A222117, A059933, A211378.
%K A215409 cons,easy,nonn,fini,full
%O A215409 0,1
%A A215409 _Jonathan Sondow_, Aug 10 2012
%E A215409 Corrected by _Natan Arie Consigli_, Jan 23 2015