OFFSET
1,1
COMMENTS
The asymptotic density of this sequence is 1 - Product_{p prime} (1 - 1/(p^2*(p+1))) = 0.1184861602... (= 1 - A065465). - Amiram Eldar, May 18 2022
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
MATHEMATICA
Select[Range[540], Count[FactorInteger[#][[All, -1]], _?(And[OddQ@ #, # > 1] &)] > 0 &] (* Michael De Vlieger, Nov 28 2017 *)
PROG
(Python)
from sympy import factorint
def ok(n):
return max((e for e in factorint(n).values() if e%2), default=-1) > 1
print(list(filter(ok, range(541)))) # Michael S. Branicky, Aug 24 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Nov 28 2017
STATUS
approved