OFFSET
1,4
COMMENTS
If n is A007369 (sigma(x) = n has no solution) then a(n) = 0.
Obviously a(n) >= A054973(n), number of solutions to sigma(x) = n.
The equality is obtained for terms of A007369, but not only: see a(n) for 1, 3, 6, 13, 18, 20, 30, 31, 36, 38 ...
The maxima for a(n) are : 1, 2, 3, 4, 5, 10, 14, 15 ... and are obtained for n: 1, 4, 7, 8, 15, 24, 60, 120, ...
LINKS
Jens Kruse Andersen, Table of n, a(n) for n = 1..10000
G. L. Cohen and H. J. J. te Riele, Iterating the sum-of-divisors function, Experimental Mathematics, 5 (1996), pp. 93-100 (see p. 97)
EXAMPLE
There is a single integer such that sigma(x) = 1 so a(1) = 1.
For n=4, we have only sigma(3) = 4 and sigma(sigma(2)) = 4, so a(4) = 2.
For n=7, we have only sigma(4) = 7, sigma(sigma(3)) = 7, and sigma(sigma(sigma(2))) = 7, so a(7) = 3.
PROG
(PARI) isok(i, n) = {j = i; while((k = sigma(j)) < n, j = k); k == n; }
a(n) = {if (n == 1, return (1)); nb = 0; for (i=2, n-1, nb += isok(i, n); ); nb; }
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Aug 09 2014
STATUS
approved