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.

Showing 1-10 of 27 results. Next

A266201 Goodstein numbers: a(n) = G_n(n), where G is the Goodstein function.

Original entry on oeis.org

0, 0, 1, 2, 83, 1197, 187243, 37665879, 20000000211, 855935016215, 44580503598539, 2120126221988686, 155568095557812625, 6568408355712901455, 295147905179358418247, 14063084452070776884879
Offset: 0

Views

Author

Natan Arie Consigli, Jan 22 2016

Keywords

Comments

To write an integer n in base-k hereditary representation, write n in ordinary base-k representation, and then do the same recursively for all exponents which are greater than k.
For example, the hereditary representation of 132132 in base-2 is:
132132 = 2^17 + 2^10 + 2^5 + 2^2
= 2^(2^4 + 1) + 2^(2^3 + 2) + 2^(2^2 + 1) + 2^2
= 2^(2^(2^2) + 1) + 2^(2^(2+1) + 2) + 2^(2^2 + 1) + 2^2.
Define B_k(n) to be the function that substitutes k+1 for all the bases of the base-k hereditary representation of n.
E.g., B_2(101) = B_2(2^(2^2 + 2) + 2^(2^2 + 1) + 2^2 + 1) = 3^(3^3 + 3) + 3^(3^3 + 1) + 3^3 + 1 = 228767924549638.
(Sometimes B_k(n) is referred to as n "bumped" from base k.)
The Goodstein function is defined as: G_k(n) = B_{k+1}(G_{k-1}(n)) - 1 with G_0(n) = n, i.e., iteration of bumping the number to the next larger base and subtracting one; see example section for instances.
Goodstein's theorem says that for any nonnegative n, the sequence G_k(n) eventually stabilizes and then decreases by 1 in each step until it reaches 0. (The subsequent values of G_k(n) < 0 are not part of the sequence.)
Named after the English mathematician Reuben Louis Goodstein (1912-1985). - Amiram Eldar, Jun 19 2021

Examples

			Compute a(5) = G_5(5):
G_0(5) = 5;
G_1(5) = B_2(G_0(5))-1 = B_2(2^2+1)-1 = (3^3+1)-1 = 27 = 3^3;
G_2(5) = B_3(G_1(5))-1 = B_3(3^3)-1 = 4^4-1 = 255 = 3*4^3+3*4^2+3*4+3;
G_3(5) = B_4(G_2(5))-1 = B_4(3*4^3+3*4^2+3*4+3)-1 = 467;
G_4(5) = B_5(G_3(5))-1 = B_5(3*5^3+3*5^2+3*5+2)-1 = 775;
G_5(5) = B_6(G_4(5))-1 = B_6(3*6^3+3*6^2+3*6+1)-1 = 1197.
		

Crossrefs

Cf. Goodstein sequences: A056004: G_1(n); A057650: G_2(n); A059934: G_3(n); A059935: G_4(n); A059936: G_5(n); A215409: G_n(3); A056193: G_n(4); A266204: G_n(5); A266205: G_n(6); A222117: G_n(15); A059933: G_n(16); A211378: G_n(19).
Weak Goodstein sequences: A137411: g_n(11); A265034: g_n(266); A266202: g_n(n); A266203: a(n) = k such that g_k(n)=0;
Bumping Sequences: A222112: B_2(n);
Other sequences: A222113.

Programs

  • PARI
    (B(n,b)=sum(i=1,#n=digits(n,b),n[i]*(b+1)^if(#nA266201(n)=for(k=1,n,n=B(n,k+1)-1);n \\ M. F. Hasler, Feb 12 2017

Extensions

Edited by M. F. Hasler, Feb 12 2017
Incorrect a(16) deleted (the correct value is ~ 2.77*10^861) by M. F. Hasler, Feb 19 2017

A059933 Goodstein sequence starting with 16: to calculate a(n+1), write a(n) in the hereditary representation in base n+2, then bump the base to n+3, then subtract 1.

Original entry on oeis.org

16, 7625597484986, 50973998591214355139406377, 53793641718868912174424175024032593379100060, 19916489515870532960258562190639398471599239042185934648024761145811, 5103708485122940631839901111036829791435007685667303872450435153015345686896530517814322070729709
Offset: 0

Views

Author

Henry Bottomley, Feb 12 2001

Keywords

Comments

See A266201 for definitions of and key links for hereditary representation and Goodstein sequences.
Goodstein's theorem shows that the Goodstein sequence G_n(k) eventually stabilizes and then decreases by 1 at each step until it reaches 0. Thereafter the values of G_n(k) < 0 are not part of the sequence. By Goodstein's theorem we conclude that G_n(k) is a finite sequence.
In this case when a(0) = G_0(16) = 16, there seems little possibility of describing how incredibly large n must be for a(n) to reach 0.

Examples

			a(0) = 16 = 2^(2^2) so a(1) = 3^(3^3)-1 = 7625597484986.
So a(1) = 2*3^(2*3^2 + 2*3 + 2) + 2*3^(2*3^2 + 2*3 + 1) + 2*3^(2*3^2 + 2*3) + 2*3^(2*3^2 + 1*3 + 2) + 2*3^(2*3^2 + 1*3 + 1) + 2*3^(2*3^2 + 1*3) + 2*3^(2*3^2 + 2) + 2*3^(2*3^2 + 1) + 2*3^(2*3^2) + 2*3^(3^2 + 2*3 + 2) + 2*3^(3^2 + 2*3 + 1) + 2*3^(3^2 + 2*3) + 2*3^(3^2 + 1*3 + 2) + 2*3^(3^2 + 1*3 + 1) + 2*3^(3^2 + 1*3) + 2*3^(3^2 + 2) + 2*3^(3^2 + 1) + 2*3^(3^2) + 2*3^(2*3 + 2) + 2*3^(2*3 + 1) + 2*3^(2*3) + 2*3^(1*3 + 2) + 2*3^(1*3 + 1) + 2*3^(1*3) + 2*3^(2) + 2*3^(1) + 2,
leading to a(2) = 2*4^(2*4^2 + 2*4 + 2) + 2*4^(2*4^2 + 2*4 + 1) + 2*4^(2*4^2 + 2*4) + 2*4^(2*4^2 + 1*4 + 2) + 2*4^(2*4^2 + 1*4 + 1) + 2*4^(2*4^2 + 1*4) + 2*4^(2*4^2 + 2) + 2*4^(2*4^2 + 1) + 2*4^(2*4^2) + 2*4^(4^2 + 2*4 + 2) + 2*4^(4^2 + 2*4 + 1) + 2*4^(4^2 + 2*4) + 2*4^(4^2 + 1*4 + 2) + 2*4^(4^2 + 1*4 + 1) + 2*4^(4^2 + 1*4) + 2*4^(4^2 + 2) + 2*4^(4^2 + 1) + 2*4^(4^2) + 2*4^(2*4 + 2) + 2*4^(2*4 + 1) + 2*4^(2*4) + 2*4^(1*4 + 2) + 2*4^(1*4 + 1) + 2*4^(1*4) + 2*4^(2) + 2*4^(1) + 1 = 2*(4^32 + 4^16 + 1)*(4^8 + 4^4 + 1)*(4^2 + 4*1)-1 = 50973998591214355139406377.
		

Crossrefs

Cf. A266201: G_n(n).
Cf. A056193: G_n(4), A056004: G_1(n), A057650 G_2(n), A056041.
Cf. A215409: G_n(3), A222117: G_n(15), A211378: G_n(19), A266204: G_n(5), A266205: G_n(6).

Programs

  • Haskell
    -- See Link
    
  • PARI
    bump(a, n) = {if (a < n, return (a)); my(pd = Pol(digits(a, n)));  my(de = vector(poldegree(pd)+1, k, k--; polcoeff(pd, k))); my(bde = vector(#de, k, k--; bump(k, n))); my(q = sum(k=0, poldegree(pd), if (c=polcoeff(pd, k), c*x^bde[k+1], 0))); return(subst(q, x, n+1)); }
    lista(nn) = {print1(a = 16, ", "); for (n=2, nn, a = bump(a, n)-1; print1(a, ", "); ); } \\ Michel Marcus, Feb 28 2016
    
  • PARI
    (B(n,b)=sum(i=1,#n=digits(n,b),n[i]*(b+1)^if(#n1,B(a,n)-1,16)) \\ M. F. Hasler, Feb 12 2017

Formula

a(n) = G_n(16), where G is the function defined in A266201.

Extensions

Definition corrected by N. J. A. Sloane, Mar 06 2006
Missing a(5) inserted and wrong a(7) replaced by Reinhard Zumkeller, Feb 13 2013
Revised by Natan Arie Consigli, Jan 23 2016
Offset changed to 0 by Nicholas Matteo, Aug 21 2019

A211378 Goodstein sequence starting with 19.

Original entry on oeis.org

19, 7625597484990, 13407807929942597099574024998205846127479365820592393377723561443721764030073546976801874298166903427690031858186486050853753882811946569946433649006084099
Offset: 0

Views

Author

Reinhard Zumkeller, Feb 13 2013

Keywords

Comments

To calculate a(n), write a(n-1) in the hereditary representation base n+1, then bump the base to n+2, then subtract 1.

Examples

			The first terms are (see Wikipedia):
a(0) = 2^2^2 + 2^1 + 2^0 = 19
a(1) = 3^3^3 + 3^1 + 3^0 - 1 = 7625597484990
a(2) = 4^4^4 + 4^1 - 1                       (155 digits)
a(3) = 5^5^5 + 3 - 1                         (2185 digits)
a(4) = 6^6^6 + 2 - 1                         (36306 digits)
a(5) = 7^7^7 + 1 - 1                         (695975 digits)
a(6) = 8^8^8 - 1                             (15151336 digits).
		

Crossrefs

Cf. A215409 (start=3), A056193 (start=4), A222117 (start=15), A059933 (start=16).

Programs

  • Haskell
    -- See Link

Extensions

Offset changed to 0 by Nicholas Matteo, Aug 21 2019

A222117 Goodstein sequence starting with 15.

Original entry on oeis.org

15, 111, 1283, 18752, 326593, 6588344, 150994943, 3524450280, 100077777775, 3138578427934, 106993479003783, 3937376861542204, 155568096352467863, 6568408356994335930, 295147905181357143919, 14063084452070776884879, 708235345355342213988445
Offset: 0

Views

Author

Reinhard Zumkeller, Feb 13 2013

Keywords

Comments

To calculate a(n+1), write a(n) in the hereditary representation base n+2, then bump the base to n+3, then subtract 1;
Compare to A222113: the underlying variants to define Goodstein sequences are equivalent.

Examples

			The first terms are:
a(0) = 2^(2+1) + 2^2 + 2^1 + 2^0 = 15;
a(1) = 3^(3+1) + 3^3 + 3^1 + 3^0 - 1 = 111;
a(2) = 4^(4+1) + 4^4 + 4^1 - 1 = 4^(4+1) + 4^4 + 3*4^0 = 1283;
a(3) = 5^(5+1) + 5^5 + 3*5^0 - 1 = 5^(5+1) + 5^5 + 2*5^0 = 18752;
a(4) = 6^(6+1) + 6^6 + 2*6^0 - 1 = 6^(6+1) + 6^6 + 1 = 326593;
a(5) = 7^(7+1) + 7^7 + 1 - 1 = 6588344;
a(6) = 8^(8+1) + 8^8 - 1 = 150994943.
		

Crossrefs

Cf. A215409 (start=3), A056193 (start=4), A059933 (start=16), A211378 (start=19).

Programs

  • Haskell
    -- See Link
    
  • PARI
    lista(nn) = {print1(a = 15, ", "); for (n=2, nn, pd = Pol(digits(a, n)); q = sum(k=0, poldegree(pd), if (c=polcoeff(pd, k), c*x^subst(Pol(digits(k, n)), x, n+1), 0)); a = subst(q, x, n+1) - 1; print1(a, ", "); ); } \\ Michel Marcus, Feb 24 2016

Extensions

Offset changed to 0 by Nicholas Matteo, Aug 21 2019

A215409 The Goodstein sequence G_n(3).

Original entry on oeis.org

3, 3, 3, 2, 1, 0
Offset: 0

Views

Author

Jonathan Sondow, Aug 10 2012

Keywords

Comments

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.
Decimal expansion of 33321/100000. - Natan Arie Consigli, Jan 23 2015

Examples

			a(0) = 3 = 2^1 + 1;
a(1) = 3^1 + 1 - 1 = 3^1 = 3;
a(2) = 4^1 - 1 = 3;
a(3) = 3 - 1 = 2;
a(4) = 2 - 1 = 1;
a(5) = 1 - 1 = 0.
		

Crossrefs

Programs

  • Haskell
    -- See Link
    
  • Mathematica
    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 *)
  • PARI
    B(n, b)=sum(i=1, #n=digits(n, b), n[i]*(b+1)^if(#nIain Fox, Dec 13 2017
    
  • 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
    
  • PARI
    first(n) = Vec(3 + 3*x + 3*x^2 + 2*x^3 + x^4 + O(x^n)) \\ Iain Fox, Dec 12 2017
    
  • PARI
    a(n) = floor(2 - (4/Pi)*atan(n-3)) \\ Iain Fox, Dec 12 2017

Formula

a(0) = a(1) = a(2) = 3; a(3) = 2; a(4) = 1; a(n) = 0, n > 4;
From Iain Fox, Dec 12 2017: (Start)
G.f.: 3 + 3*x + 3*x^2 + 2*x^3 + x^4.
E.g.f.: 3 + 3*x + (3/2)*x^2 + (1/3)*x^3 + (1/24)*x^4.
a(n) = floor(2 - (4/Pi)*arctan(n-3)), n >= 0.
(End)

Extensions

Corrected by Natan Arie Consigli, Jan 23 2015

A266204 a(n) = G_n(5), where G_n(k) is the Goodstein function defined in A266201.

Original entry on oeis.org

5, 27, 255, 467, 775, 1197, 1751, 2454, 3325, 4382, 5643, 7126, 8849, 10830, 13087, 15637, 18499, 21691, 25231, 29137, 33427, 38119, 43231, 48781, 54787, 61267, 68239, 75721, 83731, 92287, 101407, 111108, 121409, 132328, 143883, 156092, 168973, 182544, 196823
Offset: 0

Views

Author

Natan Arie Consigli, Jan 22 2016

Keywords

Examples

			G_0(5) = 5;
G_1(5) = B_2(5) - 1 = B_2(2^2 + 1) - 1 = 27;
G_2(5) = B_3(3^3) - 1 = 4^4 - 1 = 255;
G_3(5) = B_4(3*4^3 + 3*4^2 + 3*4 + 3) - 1 = 3*5^3 + 3*5^2 + 3*5 + 3 - 1 = 467.
		

Crossrefs

Cf. A056193: G_n(4), A059933: G_n(16), A211378: G_n(19), A215409: G_n(3), A222117: G_n(15), A266204: G_n(5), A266205: G_n(6), A059936: G_5(n), A266201: G_n(n).

Programs

  • PARI
    bump(a, n) = {if (a < n, return (a)); my(pd = Pol(digits(a, n)));  my(de = vector(poldegree(pd)+1, k, k--; polcoeff(pd, k))); my(bde = vector(#de, k, k--; bump(k, n))); my(q = sum(k=0, poldegree(pd), if (c=polcoeff(pd, k), c*x^bde[k+1], 0))); return(subst(q, x, n+1)); }
    lista(nn) = {print1(a = 5, ", "); for (n=2, nn, a = bump(a, n)-1; print1(a, ", "); ); } \\ Michel Marcus, Feb 28 2016

A056004 Initial step in Goodstein sequences: write n in hereditary representation base 2, bump to base 3, then subtract 1.

Original entry on oeis.org

0, 2, 3, 26, 27, 29, 30, 80, 81, 83, 84, 107, 108, 110, 111, 7625597484986, 7625597484987, 7625597484989, 7625597484990, 7625597485013, 7625597485014, 7625597485016, 7625597485017, 7625597485067, 7625597485068, 7625597485070, 7625597485071, 7625597485094
Offset: 1

Views

Author

Henry Bottomley, Aug 04 2000

Keywords

Comments

To write an integer n in base-k hereditary representation, write n in ordinary base-k representation, and then do the same recursively for all exponents which are greater than k: e.g., 2^18 = 2^(2^4 + 2) = 2^(2^(2^2) + 2). "Bump to base 3" means to replace all the 2's in that representation by 3. - M. F. Hasler, Feb 19 2017

Examples

			a(18)=7625597484989 since 18=2^(2^2)+2^1 which when bumped from 2 to 3 becomes 3^(3^3)+3^1=76255974849890 and when 1 is subtracted gives 7625597484989.
		

Crossrefs

Using G_k to denote the k-th step, this is the first in the following list: A056004: G_1(n), A057650: G_2(n), A059934: G_3(n), A059935: G_4(n), A059936: G_5(n); A266201: G_n(n); A056041.
Cf. A215409: G_n(3), A056193: G_n(4), A266204: G_n(5), A266205: G_n(6), A222117: G_n(15), A059933: G_n(16), A211378: G_n(19).
See A222112 for an alternate version.

Programs

  • Haskell
    -- See Link
    
  • PARI
    A056004(n)=sum(i=1,#n=binary(n),if(n[i],3^if(#n-i<2,#n-i,A056004(#n-i)+1)))-1 \\ See A266201 for more general code. - M. F. Hasler, Feb 19 2017

Extensions

Edited by M. F. Hasler, Feb 19 2017

A266205 a(n) = G_n(6), where G is the Goodstein function defined in A266201.

Original entry on oeis.org

6, 29, 257, 3125, 46655, 98039, 187243, 332147, 555551, 885775, 1357259, 2011162, 2895965, 4068068, 5592391, 7542974, 10003577, 13068280, 16842083, 21441506, 26995189, 33644492, 41544095, 50862597, 61783119, 74503901, 89238903, 106218405, 125689607, 147917229
Offset: 0

Views

Author

Natan Arie Consigli, Jan 23 2016

Keywords

Examples

			G_1(6) = B_2(6) - 1 = B_2(2^2 + 2) - 1 = 3^3 + 3 - 1 = 29;
G_2(6) = B_3(G_1(6)) - 1 = B_3(3^3 + 2) - 1 =  4^4 + 2 - 1 = 257;
G_3(6) = B_4(G_2(6)) - 1 = 5^5 + 1 - 1 = 3125;
G_4(6) = B_5(G_3(6)) - 1 = 6^6 - 1 = 46655;
G_5(6) = B_6(G_4(6)) - 1 = 5*7^5 + 5*7^4 + 5*7^3 + 5*7^2 + 5*7 + 5 - 1 = 98039.
		

Crossrefs

Cf. A056193: G_n(4), A059933: G_n(16), A211378: G_n(19), A215409: G_n(3), A222117: G_n(15), A266204: G_n(5), A266205: G_n(6), A059936: G_5(n), A266201: G_n(n).

Programs

  • PARI
    lista(nn) = {print1(a = 6, ", "); for (n=2, nn, pd = Pol(digits(a, n)); q = sum(k=0, poldegree(pd), if (c=polcoeff(pd, k), c*x^subst(Pol(digits(k, n)), x, n+1), 0)); a = subst(q, x, n+1) - 1; print1(a, ", "););} \\ Michel Marcus, Feb 22 2016

A271554 a(n) = G_n(7), where G is the Goodstein function defined in A266201.

Original entry on oeis.org

7, 30, 259, 3127, 46657, 823543, 16777215, 37665879, 77777775, 150051213, 273624711, 475842915, 794655639, 1281445305, 2004318063, 3051893870, 4537630813, 6604718946, 9431578931, 13238000758, 18291957825, 24917131658, 33501182551, 44504801406, 58471578053, 76038721330
Offset: 0

Views

Author

Natan Arie Consigli, Apr 10 2016

Keywords

Examples

			G_1(7) = B_2(7) - 1 = B[2](2^2 + 2 + 1) - 1 = 3^3 + 3 + 1 - 1 = 30;
G_2(7) = B_3(G_1(7)) - 1 = B[3](3^3 + 3) - 1 =  4^4 + 4 - 1 = 259;
G_3(7) = B_4(G_2(7)) - 1 = 5^5 + 3 - 1 = 3127;
G_4(7) = B_5(G_3(7)) - 1 = 6^6 + 2 - 1 = 46657;
G_5(7) = B_6(G_4(7)) - 1 = 7^7 + 1 - 1 = 823543;
G_6(7) = B_7(G_5(7)) - 1 = 8^8 - 1 = 16777215;
G_7(7) = B_8(G_6(7)) - 1 = 7*9^7 + 7*9^6 + 7*9^5 + 7*9^4 + 7*9^3 + 7*9^2 + 7*9 + 7 - 1 = 37665879.
		

Crossrefs

Cf. A056193: G_n(4), A059933: G_n(16), A211378: G_n(19), A215409: G_n(3), A222117: G_n(15), A266204: G_n(5), A266205: G_n(6), A266201: G_n(n).

Programs

  • PARI
    lista(nn) = {print1(a = 7, ", "); for (n=2, nn, pd = Pol(digits(a, n)); q = sum(k=0, poldegree(pd), if (c=polcoeff(pd, k), c*x^subst(Pol(digits(k, n)), x, n+1), 0)); a = subst(q, x, n+1) - 1; print1(a, ", "); ); }

A265034 Weak Goodstein sequence beginning with 266.

Original entry on oeis.org

266, 6590, 65601, 390750, 1679831, 5765085, 16777579, 43047173, 100000551, 214359541, 429982475, 815731628, 1475790101, 2562891818, 4294968647, 6975758960, 11019962273, 16983564926, 25600002083, 37822861652, 54875876045, 78310988018, 110075317151, 152587893847
Offset: 0

Views

Author

N. J. A. Sloane, Dec 09 2015, following a suggestion from Alexander R. Povolotsky

Keywords

Crossrefs

Extensions

More terms from Chai Wah Wu, Dec 09 2015
Showing 1-10 of 27 results. Next