OFFSET
1,4
FORMULA
a(n) = Sum_{ i <= n | Omega(i) = 2} i * (1 - mu(i)).
EXAMPLE
a(10) = 13; 4 and 9 are the only squared primes less than or equal to 10 and their sum is 13.
MAPLE
with(numtheory); a:=n->add('if'(bigomega(i)=2, i*(1-mobius(i)), 0), i=1..n); seq(a(n), n=1..100);
MATHEMATICA
With[{sp=Prime[Range[10]]^2}, Table[Total[Select[sp, #<=n&]], {n, 100}]] (* Harvey P. Dale, Nov 11 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Nov 09 2013
STATUS
approved