reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
Daniel Suteu and Robert Price, <a href="/A138963/b138963_1.txt">Table of n, a(n) for n = 1..63</a> (terms a(1)..a(45) from Robert Price)
reviewed
approved
proposed
reviewed
editing
proposed
(PARI)
(PARI)f(n) = my(D = Vec(concat(apply(s->Str(s), [1..n])))); eval(concat(vector(#D, k, D[#D-k+1]))); \\ A138793
a(1) = 1; , a(n+1) = the largest prime divisor of A138793(n+1).
1, 7, 107, 149, 953, 218107, 402859, 4877, 379721, 4349353, 169373, 182473, 1940144339383, 2184641, 437064932281, 5136696159619, 67580875919190833, 1156764458711, 464994193118899, 4617931439293, 1277512103328491957510030561, 8177269604099
Daniel Suteu and Robert Price, <a href="/A138963/b138963_1.txt">Table of n, a(n) for n = 1..63</a> (terms a(1)..a(45</a>) from Robert Price)
a(n) = A006530(A138793(n)). - Daniel Suteu, May 26 2022
(PARI)f(n) = my(D = Vec(concat(apply(s->Str(s), [1..n])))); eval(concat(vector(#D, k, D[#D-k+1]))); \\ A138793
a(n) = if(n == 1, 1, vecmax(factor(f(n))[, 1])); \\ Daniel Suteu, May 26 2022
approved
editing