OFFSET
1,4
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16384
FORMULA
a(1) = 1; for n > 1, if n = p^k for some prime p and exponent k >= 1, then a(n) = -k, otherwise [when n is not a prime power], a(n) = 1+A085970(n) = running count from 2 onward.
PROG
(PARI)
up_to = 65537;
partialsums(f, up_to) = { my(v = vector(up_to), s=0); for(i=1, up_to, s += f(i); v[i] = s); (v); }
v065515 = partialsums(n -> (omega(n)<=1), up_to);
A065515(n) = v065515[n];
CROSSREFS
KEYWORD
sign
AUTHOR
Antti Karttunen, Jul 02 2018
STATUS
approved