Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A048520
Primes expressible as the sum of a prime plus its digit sum.
14
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,1
LINKS
M. F. Hasler, Table of n, a(n) for n = 1..10000, Nov 08 2018
EXAMPLE
a(15) = 181 which is 167 + (1+6+7).
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 *)
PROG
(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
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Patrick De Geest, May 15 1999
EXTENSIONS
Offset corrected to 1 by M. F. Hasler, Nov 08 2018
STATUS
approved