A156850 a(1)=2, a(n+1) is the smallest prime > (a(n) + sum of digits of a(n)).

2, 5, 11, 17, 29, 41, 47, 59, 79, 97, 127, 139, 157, 173, 191, 211, 223, 233, 251, 263, 277, 307, 331, 347, 367, 389, 419, 439, 457, 479, 503, 521, 541, 557, 577, 599, 631, 643, 659, 683, 701, 719, 739, 761, 787, 811, 823, 839, 863, 881, 907, 929, 953, 971

Offset

1,1

Links

Table of n, a(n) for n=1..54.

Example

2, 5(>4=2+2), 11(>10=5+5), 17(>13=11+1+1), 29(>25=17+1+7), etc.

Maple

A007953 := proc(n) add(d, d=convert(n, base, 10)) ; end: A156850 := proc(n) option remember; if n = 1 then 2; else procname(n-1)+A007953(procname(n-1)) ; nextprime(%) ; fi; end: seq(A156850(n), n=1..120) ; # R. J. Mathar, Feb 20 2009

Mathematica

NestList[NextPrime[#+Total[IntegerDigits[#]]]&, 2, 60] (* Harvey P. Dale, Aug 21 2011 *)

Crossrefs

Cf. A000040, A156672.

Sequence in context: A278049 A007491 A124850 * A156611 A143509 A111166

Adjacent sequences: A156847 A156848 A156849 * A156851 A156852 A156853

Keyword

nonn,base

Author

Juri-Stepan Gerasimov, Feb 17 2009

Extensions

Definition clarified by Harvey P. Dale, Aug 21 2011

Status

approved