A069550 Squarefree part of n equals bigomega(n).

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

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

Cf. A001222, A007913.

Sequence in context: A003346 A047723 A151542 * A079705 A068721 A228181

Adjacent sequences: A069547 A069548 A069549 * A069551 A069552 A069553

Keyword

nonn

Author

Benoit Cloitre, Apr 17 2002

Status

approved