Does there exist a power-summer of order $+ \infty$?
We start with some number $n$ and sum its digits (we can denote sum-of-digits function as $S_d$) to obtain number $S_d(n)$.
If $S_d(n)$ is prime then we calculate number $n^2$ and sum its digits to obtain $S_d(n^2)$. If $S_d(n^2)$ is prime then we calculate number $n^3$ and sum its digits to obtain $S_d(n^3)$, and so on...
We can call number $n$ a power-summer of order m if the numbers $S_d(n),...S_d(n^m)$ are all primes.
We can call number $n$ a power-summer of order $+ \infty$ if $n$ is power-summer of order m for every $m \in \mathbb N$
A question is:
Does there exist a power-summer of order $+ \infty$?
Are you of the opinion that there is some global maximum, that is, a natural number $W$ such that order of every $n$ is less than $W$?
An answer is not in my reach, I do not know much about sum-of-digits functions, but maybe someone has some good ideas.
Peter found a number of order $14$, a number $20619661$ and calculated that upto $n=10^9$ there is no number with an order greater than $14$.
Solution 1:
Basic probabilistic heuristics suggest that there is no power-summer of order $+\infty$ nor is there some global maximum $W$.
If we view $n^m$ as a random number between $0$ and $9\log_{10}(n^m) = 9m\log(n)$, then $Pr(S_d(n^m) \text{ is prime}) \approx \frac{1}{S_d(n^m)} \approx \frac{1}{4.5m\log(n)}$. So assuming any slight independence between the events $S_d(n^m)$ (for $n$ fixed as $m$ ranges), which is reasonable, will give that the probability that each $S_d(n^m)$ is prime is $0$.
However, for any given positive integer $W$, $Pr(S_d(n^1),\dots,S_d(n^W) \text{ are prime}) \approx \frac{1}{4.5^W W!}\frac{1}{(\log n)^W}$, so since $\sum_{n \ge 1} \frac{1}{(\log n)^W} = +\infty$, Borel-Cantelli suggests that infinitely many positive integers are power-summer of order $W$.