|
%I #16 Dec 27 2021 08:27:31
%S 12,27,75,80,96,147,180,216,363,405,448,486,500,507,600,867,980,1008,
%T 1083,1125,1176,1350,1587,2205,2268,2420,2523,2560,2646,2800,2883,
%U 2904,3125,3380,3750,4056,4107,5043,5103,5445,5488,5547,5760,5780,6125,6300
%N Squarefree part of n equals bigomega(n).
%H Harvey P. Dale, <a href="/A069550/b069550.txt">Table of n, a(n) for n = 1..1000</a>
%t sfp[n_]:=Times@@Power@@@({#[[1]],Mod[#[[2]],2]}&/@FactorInteger[n]); Select[ Range[7000],PrimeOmega[#]==sfp[#]&] (* _Harvey P. Dale_, Aug 20 2019 *)
%o (PARI) for(n=1,10000,if(core(n)==bigomega(n),print1(n,",")))
%o (Python)
%o from sympy.ntheory.factor_ import core, primeomega
%o def ok(n): return n > 0 and core(n, 2) == primeomega(n)
%o print([k for k in range(6301) if ok(k)]) # _Michael S. Branicky_, Dec 27 2021
%Y Cf. A001222, A007913.
%K nonn
%O 1,1
%A _Benoit Cloitre_, Apr 17 2002
|
