OFFSET
1,19
LINKS
Wikipedia, Squarefree integer.
EXAMPLE
a(19) = 3 since A005117(19) = 30 and 30 has 3 semiprime divisors, namely {6, 10, 15}.
MATHEMATICA
Table[If[PrimeNu[n] == PrimeOmega[n], PrimeNu[n] (PrimeNu[n] - 1)/2, {}], {n, 200}] // Flatten
PROG
(PARI) apply(x->(sumdiv(x, d, bigomega(d)==2)), select(issquarefree, [1..300])) \\ Michel Marcus, Nov 22 2023
(Python)
from math import isqrt
from sympy import mobius, primenu
def A367452(n):
def f(x): return n+x-sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))
m, k = n, f(n)
while m != k:
m, k = k, f(k)
return (c:=primenu(m))*(c-1)>>1 # Chai Wah Wu, Aug 12 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Nov 18 2023
STATUS
approved