A069551 Numbers k such that the squarefree part of k equals A062799(k).

27, 28, 40, 63, 175, 208, 250, 847, 1053, 1183, 1216, 1408, 1656, 2023, 2527, 3125, 3703, 4600, 5887, 6727, 8125, 8496, 9016, 9583, 11767, 12943, 13851, 15463, 19116, 19652, 19663, 22264, 23600, 24367, 26047, 26568, 31096, 31213, 31423, 31744

Offset: 1

Author

Benoit Cloitre, Apr 17 2002

Keywords

• nonn

Links

• Amiram Eldar, Table of n, a(n) for n = 1..1000

Crossrefs

Cf. A007913, A062799.

Programs

Mathematica

f[p_, e_] := If[OddQ[e], p, 1]; sqf[n_] := Times @@ (f @@@ FactorInteger[n]); imomega[n_] := DivisorSum[n, PrimeNu[#] &]; Select[Range[32000], sqf[#] == imomega[#] &] (* Amiram Eldar, Apr 06 2020 *)

PARI

for(n=1,50000,if(core(n)==sumdiv(n,d,omega(d)),print1(n,",")))

Formula

k such that A007913(k) = A062799(k).