|
|
A081213
|
|
Let r(n,k) = if k=0 then n, else r(A081210(n),k-1), then a(n)=r(n, A081212(n)).
|
|
5
|
|
|
1, 2, 3, 3, 5, 6, 7, 7, 7, 10, 11, 7, 13, 14, 15, 15, 17, 14, 19, 15, 21, 22, 23, 21, 23, 26, 26, 21, 29, 30, 31, 31, 33, 34, 35, 21, 37, 38, 39, 35, 41, 42, 43, 33, 35, 46, 47, 35, 47, 46, 51, 39, 53, 39, 55, 47, 57, 58, 59, 35, 61, 62, 47, 62, 65, 66, 67, 51, 69, 70, 71, 47, 73
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
MAPLE
|
A081212r := proc(n, k)
option remember ;
if k =0 then
n;
else
end if;
end proc:
local i ;
for i from 0 do
if A081212r(n, i) = A081212r(n, i+1) then
return i ;
end if;
end do:
end proc:
end proc:
|
|
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]];
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|