A048520 Primes expressible as the sum of a prime plus its digit sum.

13, 17, 29, 47, 61, 73, 79, 83, 103, 107, 113, 137, 163, 173, 181, 191, 241, 271, 281, 293, 307, 317, 383, 397, 409, 433, 439, 443, 499, 521, 523, 563, 577, 607, 631, 641, 709, 743, 757, 809, 821, 859, 877, 881, 967, 1019, 1063, 1069, 1103, 1163, 1181

Offset: 1

Author

Patrick De Geest, May 15 1999

Keywords

• nonn
• base

Examples

			a(15) = 181 which is 167 + (1+6+7).
		

Links

• M. F. Hasler, Table of n, a(n) for n = 1..10000, Nov 08 2018

Crossrefs

Cf. A007953, A047791, A048519.

Programs

Mathematica

Sort[Select[Table[p=Prime[n];p+Total[IntegerDigits[p]],{n,195}],PrimeQ]] (* Jayanta Basu, May 03 2013 *)
Select[#+Total[IntegerDigits[#]]&/@Prime[Range[200]],PrimeQ]//Sort (* Harvey P. Dale, Sep 02 2023 *)

PARI

is_A048520(n)=#select(p->p+sumdigits(p)==n,primes([n-9*#digits(n),n-2]))&&isprime(n) \\ M. F. Hasler, Nov 08 2018

Extensions

Offset corrected to 1 by M. F. Hasler, Nov 08 2018