OFFSET
1,6
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
f:= proc(t) if not isprime(t) then (numtheory:-tau(t) +`if`(issqr(t), 1, 0))/2 - 1 fi end proc:
map(f, [$4..1000]); # Robert Israel, Feb 15 2018
PROG
(PARI) f(n) = (numdiv(n) + issquare(n))/2 - 1;
lista(nn) = forcomposite(n=1, nn, print1(f(n), ", ")) \\ Michel Marcus, Feb 16 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Corrected by Robert Israel, Feb 15 2018
STATUS
approved