OFFSET
1,5
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000
MATHEMATICA
Table[Count[Table[n i, {i, n}], _?(DivisorSigma[0, #]==4&)], {n, 90}] (* Harvey P. Dale, Jun 14 2022 *)
PROG
(PARI) a(n) = sum(i=1, n, numdiv(i*n)==4); \\ Michel Marcus, Nov 11 2020
(Python)
from collections import Counter
from sympy import factorint
def A338713(n):
f = Counter(factorint(n))
return sum(1 for i in range(1, n+1) if (l:=tuple((f+Counter(factorint(i))).values()))==(1, 1) or l==(3, )) # Chai Wah Wu, Oct 14 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Nov 11 2020
STATUS
approved