login
A069551
Numbers k such that the squarefree part of k equals A062799(k).
1
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,1
LINKS
FORMULA
k such that A007913(k) = A062799(k).
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 *)
PROG
(PARI) for(n=1, 50000, if(core(n)==sumdiv(n, d, omega(d)), print1(n, ", ")))
CROSSREFS
Sequence in context: A332927 A031171 A106126 * A022397 A042470 A042468
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Apr 17 2002
STATUS
approved