I am reading this paper written by Jeffrey C. Lagarias "An Elementary Problem Equivalent to the Riemann Hypothesis".
On page 5, it says "Fluctuations in the distribution of primes will be reflected in fluctuations in the growth rate of $\sigma(n)/n$ taken over the set of colossally abundant numbers." It is kind of soft, but I don't follow this sentence.
I understand all before this sentence. After this sentence, the author uses distribution of primes to study $\sigma(n)/n$. So I am wondering how this sentence follows from the context, but any explanation is welcome.
Anyway, the program(s) I am describing use:
$$ \delta = \frac{\log(p^{k+1} - 1) - \log(p^{k+1} - p)}{\log p} $$ given some $\delta > 0,$ the correct exponent for some prime $p$ is $$ \left\lfloor \frac{\log (p^{1 + \delta} - 1) - \log(p^\delta - 1)}{\log p} \right\rfloor \; - \; 1. $$
Take some bound $B$ on any prime power that will be used. For each prime $p$ and exponent $k \geq 1$ with $p^k < B,$ calculate $ \delta = \frac{\log(p^{k+1} - 1) - \log(p^{k+1} - p)}{\log p}.$ Make a list/file, each line has three numbers $ \delta \; \; p \; \; k \; \; .$ The printing out of $\delta$ must be done carefully... Part two, sort that file by decreasing $\delta.$ Step three, read that file, for each line print out the "bumped" prime $p,$ then the prime factorization of that CA number, and the CA number itself if it is not too large.
This method says each CA number is a prime times its predecessor. That holds true as high as anyone has computed, and follows if some conjectures of Siegel are true.
Maybe I can do this again in a couple of hours.
Alright, explicitly bounding the primes was clumsy. Here, instead, I enforce a lower bound $\delta > 0.01$
$$ \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc $$
Next morning, self-contained program. Compare OEIS
$$ \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc \bigcirc $$