|
|
A073751
|
|
Prime numbers that when multiplied in order yield the sequence of colossally abundant numbers A004490.
|
|
17
|
|
|
2, 3, 2, 5, 2, 3, 7, 2, 11, 13, 2, 3, 5, 17, 19, 23, 2, 29, 31, 7, 3, 37, 41, 43, 2, 47, 53, 59, 5, 61, 67, 71, 73, 11, 79, 2, 83, 3, 89, 97, 13, 101, 103, 107, 109, 113, 127, 131, 137, 139, 2, 149, 151, 7, 157, 163, 167, 17, 173, 179, 181, 191, 193, 197, 199, 19, 211, 3
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
The Mathematica program presents a very fast method of computing the factors of colossally abundant numbers. The 100th number has a sigma(n)/n ratio of 10.5681.
This calculation assumes that the ratio of consecutive colossally abundant numbers is always prime, which is implied by a conjecture mentioned in Lagarias's paper.
The ratio of consecutive colossally abundant numbers is prime for at least the first 10^7 terms. The (10^7)-th term is a 77908696-digit number which has a sigma(n)/n value of 33.849.
Alaoglu and Erdős's paper proves that the quotient of two consecutive colossally abundant numbers is either a prime or the product of two distinct primes.
First occurrence of the n-th prime: 1, 2, 4, 7, 9, 10, 14, 15, 16, 18, 19, 22, 23, 24, 26, 27, 28, 30, 31, 32, ..., .
Positions of 2: 1, 3, 5, 8, 11, 17, 25, 36, 51, 77, 114, 178, 282, 461, 759, 1286, 2200, 3812, 6664, ..., .
Positions of 3: 2, 6, 12, 21, 38, 68, 132, 271, 595, 1356, 3191, 7775, ..., . (End)
|
|
LINKS
|
|
|
MATHEMATICA
|
pFactor[f_List] := Module[{p=f[[1]], k=f[[2]]}, N[Log[(p^(k+2)-1)/(p^(k+1)-1)]/Log[p]]-1]; maxN=100; f={{2, 1}, {3, 0}}; primes=1; lst={2}; x=Table[pFactor[f[[i]]], {i, primes+1}]; For[n=2, n<=maxN, n++, i=Position[x, Max[x]][[1, 1]]; AppendTo[lst, f[[i, 1]]]; f[[i, 2]]++; If[i>primes, primes++; AppendTo[f, {Prime[i+1], 0}]; AppendTo[x, pFactor[f[[ -1]]]]]; x[[i]]=pFactor[f[[i]]]]; lst
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|