OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
MATHEMATICA
sfp[n_]:=Times@@Power@@@({#[[1]], Mod[#[[2]], 2]}&/@FactorInteger[n]); Select[ Range[7000], PrimeOmega[#]==sfp[#]&] (* Harvey P. Dale, Aug 20 2019 *)
PROG
(PARI) for(n=1, 10000, if(core(n)==bigomega(n), print1(n, ", ")))
(Python)
from sympy.ntheory.factor_ import core, primeomega
def ok(n): return n > 0 and core(n, 2) == primeomega(n)
print([k for k in range(6301) if ok(k)]) # Michael S. Branicky, Dec 27 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Apr 17 2002
STATUS
approved