OFFSET
1,1
COMMENTS
The next term has more than 6655555 digits, because s(a(4)) >= a(3) = 59899999 and 59899999/9 > 6655555.
s(a(2)) = a(1) and s(a(3)) = a(2). Can one prove that s(a(n+1)) = a(n), for all n > 0? (Cf. formula below.) If yes, then a(n+1) is the smallest prime p with s(p) = a(n).
A046704 is primes p with s(p) also prime. A207294 is primes p with s(p) and s(s(p)) also prime. A070027 is primes p with all s(p), s(s(p)), s(s(s(p))), ... also prime. A104213 is primes p with s(p) not prime. A207293 is primes p with s(p) also prime, but not s(s(p)). A213354 is primes p with s(p) and s(s(p)) also prime, but not s(s(s(p))).
FORMULA
a(n) <= s(a(n+1)). (Proof: a(n) and s(a(n+1)) share the same property, but a(n) is minimal.)
EXAMPLE
s(13) = 1+3 = 4 is not prime, and s(p) is prime if p < 13, so a(1) = 13.
s(67) = 6+7 = 13 is prime, but s(s(67)) = s(13) = 1+3 = 4 is not prime, and no p < 67 has this property, so a(2) = 67.
s(59899999) = 5+9+8+9+9+9+9+9 = 67 and s(s(59899999)) = s(67) = 6+7 = 13 are prime, but s(s(s(59899999))) = s(13) = 1+3 = 4 is not prime, and no p < 59899999 has this property, so a(3) = 59899999.
CROSSREFS
KEYWORD
nonn,base,more,bref
AUTHOR
Jonathan Sondow, Jun 10 2012
STATUS
approved