OFFSET
1,4
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..10000
EXAMPLE
a(10) = 3, 10 -> 6 -> 3 -> 1.
a(100) = 9; f(100) = 100 - 25 = 75, f(75) = 75 - 21 = 54, f(54) = 54 - 16 = 38, f(38) = 38 - 12 = 26, f(26) = 26 - 9 = 17, f(17) = 17 - 7 = 10, f(10) = 10 - 4 = 6, f(6) = 6 - 3 = 3, f(3) = 3 - 2 = 1.
MATHEMATICA
pi[n_] := pi[n] = PrimePi[n];
A090532[n_] := Length[FixedPointList[# - pi[#] &, n]] - 2;
Array[A090532, 100] (* Paolo Xausa, Dec 13 2025 *)
PROG
(PARI) a(n) = my(nb=0); while(n != 1, nb++; n -= primepi(n)); nb; \\ Michel Marcus, Dec 14 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Dec 07 2003
EXTENSIONS
Corrected and extended by Sam Handler (sam_5_5_5_0(AT)yahoo.com), Dec 11 2004
a(1) = 0 prepended by Paolo Xausa, Dec 13 2025
STATUS
approved
