OFFSET
0,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
MAPLE
with(numtheory):
d:= n-> n*add(i[2]*pi(i[1])/i[1], i=ifactors(n)[2]):
a:= proc() local t, a; t, a:= -1, proc() -1 end;
proc(n) local h;
while a(n) = -1 do
t:= t+1; h:= d(t);
if a(h) = -1 then a(h):= t fi
od; a(n)
end
end():
seq(a(n), n=0..100);
MATHEMATICA
A258851[n_] := If[n == 0, 0, n*Total[Last[#]*PrimePi[First[#]]/First[#]& /@ FactorInteger[n]]];
a[n_] := For[m = 0, True, m++, If[A258851[m] == n, Return[m]]];
Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Sep 10 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jun 12 2015
STATUS
approved