OFFSET
1,2
COMMENTS
The data was derived from the data in A007691.
The Mersenne primes A000668 are a subsequence; in fact a prime is in this sequence if and only if it is a Mersenne prime.
Note that John Voight's 'spoof odd perfect number' 22017975903 is included. - Peter Luschny, Oct 30 2020
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1600
BYU Computational Number Theory Group, Odd, spoof perfect factorizations, arXiv:2006.10697 [math.NT], 2020.
J. Voight, On the nonexistence of odd perfect numbers.
MAPLE
a_list := proc(len) local L, n; L := NULL; for n from 1 to len do
if numtheory:-sigma(n) mod n = 0 then L := L, n/2^padic[ordp](n, 2) fi od; L end:
MATHEMATICA
A007691 = Select[Range[1000000], Divisible[DivisorSigma[1, #], #] &]; Table[A007691[[n]]/2^IntegerExponent[A007691[[n]], 2], {n, 1, Length[A007691]}] (* Vaclav Kotesovec, Oct 14 2018 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Oct 03 2018
STATUS
approved