OFFSET
1,1
COMMENTS
LINKS
Ely Golden, Table of n, a(n) for n = 1..7144
PROG
(PARI) gpf(n) = vecmax(factor(n)[, 1]);
f(n) = if (n==1, 2, n*gpf(n+1)/gpf(n)); \\ A126288
incycle(n, list) = {my(v=Vec(list)); #select(x->(x==n), v); }
cycle(n) = {my(list = List(), repeat=1); while(repeat, n = f(n); if (incycle(n, list), repeat=0); listput(list, n); ); list; }
isok(n) = {my(list = cycle(n)); incycle(n, list); } \\ Michel Marcus, Sep 08 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Ely Golden, Sep 05 2020
STATUS
approved