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