OFFSET
1,1
COMMENTS
If s is a squarefree number, then s^2 is a term if Product_{p | s} (1 + 1/p) > 2.
The primitive terms of A391427: if k is a term, then k*m is a term in A391427 for any squarefree number m that is coprime to k.
Any term in A391427 is of the form k*m where k is a term in this sequence and m is a squarefree number coprime to k. Therefore, A391427 can be generated from this sequence by multiplying with coprime squarefree numbers, and the asymptotic density of A391427 can be evaluated from the terms in this sequence (see the Comments section of A391427).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = A087248(n)^2.
MATHEMATICA
Select[Range[700], SquareFreeQ[#] && DivisorSigma[-1, #] > 2 &]^2
PROG
(PARI) isok(k) = {my(r); if(!issquare(k, &r), 0, my(f = factor(r)); issquarefree(f) && sigma(f, -1) > 2); }
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Dec 09 2025
STATUS
approved
