OFFSET
3,2
COMMENTS
a(17) = 2.066...*10^4574. - Pontus von Brömssen, Sep 25 2020
LINKS
Pontus von Brömssen, Table of n, a(n) for n = 3..16
Wikipedia, Goodstein's theorem
EXAMPLE
Compute G_9(10):
G_1(10)= B_2(10)-1 = B_2(2^(2+1)+2)-1 = 3^(3+1)+3-1 = 83;
G_2(10) = B_3(3^(3+1)+2)-1 = 4^(4+1)+2-1 = 1025;
G_3(10) = B_4(4^(4+1)+1)-1 = 5^(5+1)+1-1 = 15625;
G_4(10) = B_5(5*5^(5+1))-1 = 6^(6+1)-1= 279935;
G_5(10) = B_6(5*6^6+5*6^5+5*6^4+5*6^3+5*6^2+5*6+5)-1 = 5*7^7+5*7^5+5*7^4+5*7^3+5*7^2+5*7+5-1 = 4215754;
G_6(10) = B_7(5*7^7+5*7^5+5*7^4+5*7^3+5*7^2+5*7+4)-1 = 5*8^8+5*8^5+5*8^4+5*8^3+5*8^2+5*8+4-1 = 84073323;
G_7(10) = B_8(5*8^8+5*8^5+5*8^4+5*8^3+5*8^2+5*8+3)-1 = 5*9^9+5*9^5+5*9^4+5*9^3+5*9^2+5*9+3-1 = 1937434592;
G_8(10) = B_9(5*9^9+5*9^5+5*9^4+5*9^3+5*9^2+5*9+2)-1 = 5*10^10+5*10^5+5*10^4+5*10^3+5*10^2+5*10+2-1 = 50000555551;
G_9(10) = B_10(5*10^10+5*10^5+5*10^4+5*10^3+5*10^2+5*10+1)-1 = 5*11^11+5*11^5+5*11^4+5*11^3+5*11^2+5*11+1-1 = 1426559238830.
PROG
(Python)
from sympy.ntheory.factor_ import digits
def bump(n, b):
s=digits(n, b)[1:]
l=len(s)
return sum(s[i]*(b+1)**bump(l-i-1, b) for i in range(l) if s[i])
def A271985(n):
if n==3: return 0
for i in range(2, 11):
n=bump(n, i)-1
return n # Pontus von Brömssen, Sep 25 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Natan Arie Consigli, Apr 30 2016
EXTENSIONS
Incorrect program and terms removed by Pontus von Brömssen, Sep 25 2020
STATUS
approved