|
|
A323757
|
|
Modified exponential perfect numbers: numbers k such that A241405(k) = 2*k.
|
|
4
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Also in the sequence are 1028004440830371164160, 20546724596095746048000, and 146361946186458562560000 (corresponding to the 3 additional terms of A054980 given by Andrew Lelechenko). - Amiram Eldar, Jul 18 2019
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
f[p_, e_] := DivisorSum[e+1, p^(#-1)&]; mesigma[1]=1; mesigma[n_] := Times @@ f @@@FactorInteger@n; mePerfectQ[n_] := mesigma[n]==2n; Select[Range[10000], mePerfectQ]
|
|
PROG
|
(PARI) f(n) = {my(f=factor(n)); prod(i=1, #f[, 1], sumdiv(f[i, 2]+1, d, f[i, 1]^(d-1))); } \\ A241405
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|