login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A319955
a(1)=1; thereafter a(n+1) = a(n) + Product_{b=2..n+1} b^(1+floor(log_b(a(i)))).
0
1, 3, 39, 331815, 114126085737676800331815
OFFSET
1,2
LINKS
Jeffrey Shallit, Automaticity IV: Sequences, sets, and diversity, J. Théor. Nombres Bordeaux 8.2 (1996): 347-367. See Theorem 10 p. 353, but b<=i+2 is a typo.
PROG
(PARI) a(n) = if (n==1, 1, a(n-1) + prod(b=2, n, b^(1+logint(a(n-1), b)))); \\ Michel Marcus, Jan 04 2021
CROSSREFS
Sequence in context: A097421 A180418 A166999 * A317484 A196159 A084186
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Oct 04 2018
EXTENSIONS
Name corrected by Michel Marcus, Jan 04 2021
STATUS
approved