OFFSET
1,3
COMMENTS
Similar to A125087, but instead of exponents, we use number of proper divisors.
LINKS
Katarzyna Matylla, Table of n, a(n) for n = 1..1000
EXAMPLE
a(12)=0 because 12 has 5 proper divisors (1, 2, 3, 4 and 6) and there is no 5 in a(1), a(2), ..., a(11).
PROG
(Maxima) max:1000; f:makelist(0, i, 1, max); apr:makelist(0, i, 0, max); f[1]:1; apr[2]:1; for n:2 through max do block(f[n]:apr[divsum(n, 0)], apr[f[n]+1]:apr[f[n]+1]+1);
CROSSREFS
KEYWORD
nonn
AUTHOR
Katarzyna Matylla, Feb 25 2008
STATUS
approved