OFFSET
1,1
REFERENCES
József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, Chapter 1, p. 53.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
J. Hanumanthachari, V. V. Subrahmanya Sastri, and V. Srinivasan, On e-perfect numbers, Math. Student, Vol. 46, No. 1 (1978), pp. 71-80; entire issue.
FORMULA
300 is in the sequence since esigma(300) = 540, and esigma(540) = 900 = 3*300.
MATHEMATICA
f[p_, e_] := DivisorSum[e, p^# &]; esigma[1] = 1; esigma[n_] := Times @@ f @@@ FactorInteger[n]; espQ[n_] := esigma[esigma[n]] == 3n; Select[Range[50000], espQ]
PROG
(PARI) esigma(n) = {my(f = factor(n)); prod(k = 1, #f~, sumdiv(f[k, 2], d, f[k, 1]^d)); }
isok(k) = esigma(esigma(k)) == 3*k; \\ Amiram Eldar, Jan 09 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 04 2019
STATUS
approved