login
A126974
a(1)=0. a(2)=1. a(n) = a(d(n)) + a(phi(n)), where d(n) = A000005(n), phi(n) = A000010(n).
1
0, 1, 2, 3, 4, 4, 5, 6, 6, 6, 7, 7, 8, 7, 9, 10, 11, 8, 9, 10, 10, 9, 10, 12, 12, 10, 11, 11, 12, 12, 13, 14, 13, 13, 15, 13, 14, 11, 15, 16, 17, 13, 14, 14, 16, 12, 13, 16, 15, 14, 17, 16, 17, 14, 19, 18, 16, 14, 15, 17, 18, 15, 17, 19, 19, 16, 17, 18, 17, 18, 19, 19, 20, 16, 20
OFFSET
1,3
COMMENTS
Question: With what numbers the gaps (horizontal stripes) in the scatter plot are associated? - Antti Karttunen, Jul 12 2017
LINKS
MAPLE
with(numtheory): a[1]:=0: a[2]:=1: for n from 3 to 90 do a[n]:=a[tau(n)]+a[phi(n)] od: seq(a[n], n=1..90); # Emeric Deutsch, Mar 24 2007
MATHEMATICA
f[l_List] := Block[{n = Length[l] + 1}, Append[l, l[[DivisorSigma[0, n]]] + l[[EulerPhi[n]]]]]; Nest[f, {0, 1}, 75] (* Ray Chandler, Mar 24 2007 *)
a[n_] := a[n] = If[n <= 2, n - 1, a[DivisorSigma[0, n]] + a[EulerPhi@ n]]; Array[a, 75] (* Michael De Vlieger, Jul 12 2017 *)
PROG
(PARI) A126974(n) = if(n<3, n-1, A126974(numdiv(n))+ A126974(eulerphi(n))); \\ Antti Karttunen, Jul 12 2017
CROSSREFS
Sequence in context: A375626 A335572 A268680 * A354459 A089058 A282717
KEYWORD
nonn,look
AUTHOR
Leroy Quet, Mar 20 2007
EXTENSIONS
More terms from Emeric Deutsch and Ray Chandler, Mar 24 2007
STATUS
approved