login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A329879
Numbers k such that k and nusigma(k) have the same set of prime divisors, where nusigma(k) is the sum of nonunitary divisors of k (A048146).
2
4, 9, 24, 25, 49, 54, 112, 121, 150, 169, 289, 294, 361, 480, 529, 726, 750, 841, 961, 1014, 1369, 1681, 1734, 1849, 1984, 2058, 2166, 2209, 2430, 2520, 2688, 2809, 3174, 3481, 3721, 3780, 4489, 5041, 5046, 5329, 5760, 5766, 6241, 6889, 7921, 7986, 8214, 8700
OFFSET
1,1
COMMENTS
Numbers k such that rad(nusigma(k)) = rad(k), where rad(k) is the squarefree kernel of k (A007947).
LINKS
MATHEMATICA
rad[n_] := Times @@ (First@# & /@ FactorInteger@ n); usigma[1] = 1; usigma[n_] := Times @@ (1 + Power @@@ FactorInteger[n]); nusigma[n_] := DivisorSigma[1, n] - usigma[n]; Select[Range[10^4], rad[nusigma[#]] == rad[#] &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 23 2019
STATUS
approved