|
|
A258862
|
|
Second pi-based antiderivative of n: the least m such that A258851^2(m) equals n.
|
|
6
|
|
|
0, 3, 5, 11, 4, 10, 41, 13, 15, 83, 109, 29, 19, 35, 191, 43, 30, 277, 74, 14, 8, 42, 77, 431, 461, 21, 22, 563, 66, 599, 26, 78, 12, 61, 141, 163, 877, 18, 214, 218, 226, 38, 114, 201, 105, 1201, 215, 1297, 302, 55, 1447, 1471, 89, 25, 103, 170, 58, 291, 51
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
FORMULA
|
a(n) = min { m >= 0 : A258851^2(m) = n }.
|
|
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(d(t));
if a(h) = -1 then a(h):= t fi
od; a(n)
end
end():
seq(a(n), n=0..100);
|
|
MATHEMATICA
|
d[n_] := d[n] = If[n == 0, 0, n*Total[Last[#]*PrimePi[First[#]]/First[#]& /@ FactorInteger[n]]];
A[n_, k_] := For[m = 0, True, m++, If[Nest[d, m, k] == n, Return[m]]];
a[n_] := A[n, 2];
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|