OFFSET
1,4
COMMENTS
From Antti Karttunen, Mar 30 2021: (Start)
The original definition was: a(n) = Max_{k=1..n} gcd(n, A098371(k)).
but this generates terms that differ from ones in the data section at least at n=14, 26, 38, 74, 106.
The original author also commented that: A098371(n) = #{k: a(n) = gcd(n, A098371(k))}, but please see my notes in that sequence.
(End)
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..8192
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
PROG
(PARI) A098372list(up_to) = { my(v=vector(up_to), y=vector(up_to)); v[1] = 1; y[1] = 1; for(n=2, up_to, my(m=1); for(k=1, n-1, m = max(m, gcd(n, v[k]))); y[n] = m; v[n] = sum(k=1, n-1, m == gcd(n, v[k]))); (y); }; \\ Antti Karttunen, Mar 30 2021
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Reinhard Zumkeller, Sep 05 2004
EXTENSIONS
Definition corrected to match with the data. - Antti Karttunen, Mar 30 2021
STATUS
approved