login
A069550
Squarefree part of n equals bigomega(n).
1
12, 27, 75, 80, 96, 147, 180, 216, 363, 405, 448, 486, 500, 507, 600, 867, 980, 1008, 1083, 1125, 1176, 1350, 1587, 2205, 2268, 2420, 2523, 2560, 2646, 2800, 2883, 2904, 3125, 3380, 3750, 4056, 4107, 5043, 5103, 5445, 5488, 5547, 5760, 5780, 6125, 6300
OFFSET
1,1
LINKS
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
Sequence in context: A003346 A047723 A151542 * A079705 A068721 A228181
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Apr 17 2002
STATUS
approved