A330163 Even perfect numbers m from A000396 such that w = (m + 2^(k(m) - 1) - 1) * 2^(2*(k(m) - 1)) is also an even perfect number, where k(m) is the Mersenne exponent A000043(m).
6, 28, 8128, 2305843008139952128
Offset: 1
Programs
-
Magma
[(2^k - 1) * (2^(k - 1)): k in [1..100] | SumOfDivisors((2^k - 1) * (2^(k - 1))) / ( (2^k - 1) * (2^(k - 1))) eq 2 and SumOfDivisors(((2^k - 1) * (2^(k - 1)) + (2^(k - 1) - 1)) * (2^(2*(k - 1)))) / (((2^k - 1) * (2^(k - 1)) + (2^(k - 1) - 1)) * (2^(2*(k - 1)))) eq 2]
-
Mathematica
f[n_] := 2^(n - 1)*(2^n - 1); g[n_] := 2^n - 2^((n - 1)/2); mers = MersennePrimeExponent[Range[10]]; g /@ Select[mers, MemberQ[f /@ mers, g[#]] &] (* Amiram Eldar, Dec 06 2019 *)
Comments