OFFSET
1,6
COMMENTS
Number of divisors of n that are not positive powers of primes (cf. A000961). - Benoit Cloitre, May 03 2003; corrected Dec 16 2008 at the suggestion of Ray Chandler.
a(n) = 1 iff n is in A000961. - Robert Israel, Nov 23 2015
a(n) = 2 iff n is in A006881. - Altug Alkan, Nov 23 2015
a(n) = 3 iff n is in A054753. - Michel Marcus, Nov 24 2015
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
After first statement in comment section, a(60) = 8 because we have: 1,6,10,12,15,20,30,60. The divisors 2,3,4,5 are excluded from the count. - Geoffrey Critzer, Nov 22 2015
MAPLE
seq(numtheory:-tau(n)-numtheory:-bigomega(n), n=1..300); # Robert Israel, Nov 23 2015
MATHEMATICA
Table[DivisorSigma[0, n] - PrimeOmega[n], {n, 1, 105}] (* Geoffrey Critzer, Nov 22 2015 *)
PROG
(PARI) first(m)=vector(m, n, numdiv(n) - bigomega(n)) \\ Anders Hellström, Nov 22 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladeta Jovovic, Apr 03 2003
STATUS
approved