OFFSET
1,2
COMMENTS
A number n is in this sequence if the sequence defined by x(k+1) = A339541(k) with x(0)=n has more initial primes than the sequences for smaller n.
No more terms < 10^8.
EXAMPLE
MAPLE
sod:= (n.b) -> convert(convert(n, base, b), `+`):
f:= n -> n + sod(n, sod(n, 10)):
g:= proc(n) option remember;
if isprime(n) then 1 + procname(f(n))
else 0
fi
end proc:
R:= 1: vmax:= 0: p:= 1:
while p < 10^7 do
p:= nextprime(p);
v:= g(p);
if v > vmax then
R:= R, p; vmax:= v;
fi
od:
R;
CROSSREFS
KEYWORD
nonn,base,bref,more
AUTHOR
J. M. Bergot and Robert Israel, Dec 08 2020
STATUS
approved