|
|
A155043
|
|
a(0)=0; for n >= 1, a(n) = 1 + a(n-d(n)), where d(n) is the number of divisors of n (A000005).
|
|
55
|
|
|
0, 1, 1, 2, 2, 3, 2, 4, 3, 3, 3, 4, 3, 5, 4, 5, 5, 6, 4, 7, 5, 7, 5, 8, 6, 6, 6, 9, 6, 10, 6, 11, 7, 11, 7, 12, 10, 13, 8, 13, 8, 14, 8, 15, 9, 14, 9, 15, 9, 10, 10, 16, 10, 17, 10, 17, 10, 18, 11, 19, 10, 20, 12, 19, 19, 21, 12, 22, 13, 22, 13, 23, 11, 24, 14, 23, 14, 25, 14, 26, 14, 15, 15
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,4
|
|
COMMENTS
|
Number of steps needed to reach zero when starting from k = n and repeatedly applying the map that replaces k by k - d(k), where d(k) is the number of divisors of k (A000005).
The original name was: a(n) = 1 + a(n-sigma_0(n)), a(0)=0, sigma_0(n) number of divisors of n.
(End)
|
|
LINKS
|
|
|
FORMULA
|
a(0) = 0; for n >= 1, a(n) = 1 + a(A049820(n)).
Other identities. For all n >= 0:
(End)
|
|
MAPLE
|
with(numtheory): a := proc (n) if n = 0 then 0 else 1+a(n-tau(n)) end if end proc: seq(a(n), n = 0 .. 90); # Emeric Deutsch, Jan 26 2009
|
|
MATHEMATICA
|
a[0] = 0; a[n_] := a[n] = 1 + a[n - DivisorSigma[0, n]]; Table[a@n, {n, 0, 82}] (* Michael De Vlieger, Sep 24 2015 *)
|
|
PROG
|
(PARI)
v155043 = vector(uplim);
v155043[1] = 1; v155043[2] = 1;
for(i=3, uplim, v155043[i] = 1 + v155043[i-numdiv(i)]);
A155043 = n -> if(!n, n, v155043[n]);
for(n=0, uplim, write("b155043.txt", n, " ", A155043(n)));
(Haskell)
import Data.List (genericIndex)
a155043 n = genericIndex a155043_list n
a155043_list = 0 : map ((+ 1) . a155043) a049820_list
(Python)
from sympy import divisor_count as d
def a(n): return 0 if n==0 else 1 + a(n - d(n))
|
|
CROSSREFS
|
Cf. A262518, A262519 (bisections, compare their scatter plots), A262521 (where the latter is less than the former).
Cf. A262507 (number of times n occurs in total), A262508 (values occurring only once), A262509 (their indices).
Cf. A263265 (nonnegative integers arranged by the magnitude of a(n)).
Cf. A264893 (first differences), A264898 (where repeating values occur).
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|