We have a new paper in which we prove that if N is an OPN,
then \Omega(N) >= max(101, 2\omega(N)+51, (18\omega(N)-31)/7).
Michael Rao wrote a program producing a tree similar to the one in
[Brent, Cohen, te Riele]
that proves that an odd perfect number (OPN for short) is > 10^300.
We use it to prove that:
- an OPN is greater than 10^1500 (since then, we have pushed the computation to 10^2000),
- the total number of prime factors (counting multiplicities) of an OPN is at least 101,
- the largest component (i.e. factor p^a with p prime) of an OPN is > 10^62.
A draft of our paper is here.
The proof tree looks like this (26Mb for a bound of 10^540).
The slides of a talk in Montpellier (2014).
Differences with [BCR]:
- We forbid the factors 127,19,7,11,331,31,97,61,13,398581,1093,3,5,307,17 instead of 127,19,7,11,31,13,3,5.
- We branch on the overall largest available prime factor instead of branching on the largest available prime factor "from the previous line".
- We do not use the arguments of the B25 and B30 bounds described in [BCR].
- We use instead a method to circumvent the roadblocks similar to the one used by Kevin Hare
to show that an OPN has at least 75 not necessarily distinct prime factors. It is described below.
- In the end, we have to conclude that an odd perfect number not divisible by any of the forbidden factors must be greater than the target bound.
To do this, we use an improved version of the argument in [BCR] that is described below.
t550 (07/2017) t600 (07/2017) t800 (07/2017)
t1000 (07/2017) t1200 (07/2017) t1600 (07/2017) t2100 (07/2017)
tXXXX contains composite numbers encountered when targetting the lower bound 10^XXXX.
The format is p q c where c is a composite factor of sigma(p^q).
There are all and only the composites that appear firstly on their respective branch:
If you factor one of them, then the proof tree will be reduced, and this factorization will not become useless
because of the factorization of some other number.
A list of the most difficult composites for the proof of N > 10^2100 (09/2017).
The format for p^q-1 is "p q-1 composite weight", where "weight" is the number of roadblocks involving the composite.
The format for prime factors in checkfacts.txt (01/2017, 2.5 Gb) is the same as in Brent's tables:
p q+1 f where f is a factor of sigma(p^q). p, q+1, and f are prime.
The compressed file checkfacts.txt.gz (01/2017) is 660 Mb.
If you want to help us finding factors
Look at the threads about odd perfect numbers or OPN
at the Mersenne forum.
Given a roadblock, we compute an upper bound on the smallest not-yet-considered factor of the OPN.
Then we branch on all prime factors below this bound (except the forbidden or already considered ones).
Example: Consider the roadblock sigma(7^4) = 2801, sigma(2801^82) = C283.
We check that C283 has no small factors, say less than 10^9.
So the abundancy of the known factors is less than (1+1/6)*(1+1/2800)*(1+1/10^9)^(283/9)<1.1671.
Suppose our target bound is 10^1500 and p is the smallest prime factor of the OPN distinct from 7 and 2801.
Suppose p=641. We must have at least 340 other prime factors >= 641 because 1.1671*(1+1/640)^340 < 2.
But then 7^4*2801^164*641^340 > 10^1523.
So we have to rule out every prime factor < 641 except 19, 127 (forbidden) and 7 (already considered).
When branching on a prime in order to circumvent a roadblock, we might encounter a new roadblock, so we have
to apply this method recursively.
Once enough small primes are forbidden, [BCR] considers the contribution (multiplicative increasing)
of each prime p to both the abundancy that has to reach the value 2 and the OPN that should remain smaller than 10^300.
Smaller primes p are prefered since their contribution to the abundancy (less than p/(p-1)) is higher and
their contribution to the OPN (at least p^2, except for the special prime) is lower.
With these rules, the best thing to do is to take every allowed prime in order.
Then, the OPN gets bigger than 10^300 before that the abundancy reaches 2.
As the target bound increases, we forbid more small primes so that this argument still works.
We also refine the argument. We consider the component (i.e., the prime together with its exponent)
rather than the prime alone, and we use a bit more the fact that some factors are forbidden.
The efficiency of the component p^a is defined as eff(p,a) = ln(sigma(p^a)/p^a)/ln(p^a).
We deal with multiplicative increasings, which expains the logarithms, and then make the ratio between the contribution
of the component p^a to the abundancy and to the OPN.
a < b implies eff(p,a) > eff(p,b)
p < q implies eff(p,a) > eff(q,a)
- for each allowed prime p, we find the smallest a such that sigma(p^a) is not divisible by 4 or some forbidden factor.
Example: Consider p=23. sigma(23^1), sigma(23^2), sigma(23^3), are respectively divisible by 4, 7, 4. So the exponent of 23 is at least 4.
- we sort these components p^a by decreasing efficiency.
- we compute the product of such components as long as the abundancy of the product remains smaller than 2.
This product is greater than 10^1735, as shown by this maple code.
Older stuff (not used in this work):
abundant.txt is an idea to reduce a little the proof tree. E_q_k.c is a C/GMP code to compute the set E(q,k) in [BCR].