OFFSET
1,2
COMMENTS
LINKS
R. Zumkeller, Table of n, a(n) for n = 1..10000
MAPLE
A081212r := proc(n, k)
option remember ;
if k =0 then
n;
else
procname(A081210(n), k-1) ;
end if;
end proc:
A081212 := proc(n)
local i ;
for i from 0 do
if A081212r(n, i) = A081212r(n, i+1) then
return i ;
end if;
end do:
end proc:
A081213 := proc(n)
A081212r(n, A081212(n)) ;
end proc:
seq(A081213(n), n=1..84) ; # R. J. Mathar, May 25 2023
MATHEMATICA
gsf[n_] := For[k = n, True, k--, If[SquareFreeQ[k], Return[k]]];
A081210[n_] := (cnt++; Times @@ gsf /@ Power @@@ FactorInteger[n]);
r[n_, k_] := r[n, k] = If[k == 0, n, r[A081210[n], k - 1]];
a[n_] := r[n, A081212[n]];
Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Sep 12 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Mar 10 2003
STATUS
approved