OFFSET
1,1
COMMENTS
The asymptotic density of this sequence is A065463 - A059956 = 0.09651509914... . - Amiram Eldar, Sep 27 2023
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
FORMULA
EXAMPLE
24 = 2^3 * 3^1 is in this sequence because it has 2 distinct prime factors whose multiplicities are odd and one such multiplicity exceeds 1.
MATHEMATICA
Select[Select[Range[1000], Nor[SquareFreeQ[#], PrimePowerQ[#]] &], Times @@ FactorInteger[#][[All, 1]] == (Sqrt[#] /. (c_ : 1)*a_^(b_ : 0) :> (c*a^b)^2) &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael De Vlieger, Sep 08 2023
STATUS
approved