login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A048519 Prime plus its digit sum equals a prime. 29
11, 13, 19, 37, 53, 59, 71, 73, 97, 101, 103, 127, 149, 163, 167, 181, 233, 257, 271, 277, 293, 307, 367, 383, 389, 419, 431, 433, 479, 499, 509, 547, 563, 587, 617, 631, 701, 727, 743, 787, 811, 839, 857, 859, 947, 1009, 1049, 1061, 1087, 1153, 1171 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
For any prime p, p +- digitsum(p, base b) can't be prime unless the base b is even, since in an odd base, an odd number always has an odd digit sum (powers of b are congruent to b (mod 2)), so p +- digitsum(p, base b) is even for odd b. This sequence is for b = 10 (where "-" is also excluded, see comment in A243442), see A243441 for b = 2. - M. F. Hasler, Nov 06 2018
See subsequence A048523 for primes which only once give another prime under iteration of A062028, and A048524 .. A048527, A320878 .. A320880 for primes starting longer chains. See A090009 for their initial terms, starting the earliest chain of given length. - M. F. Hasler, Nov 09 2018
LINKS
FORMULA
Primes in A047791, i.e., intersection of A047791 and A000040. - M. F. Hasler, Nov 08 2018
EXAMPLE
a(9) = prime 97 because 97 + sum-of-digits(97) = 97 + 16 = 113 also a prime.
MAPLE
select(n -> isprime(n) and isprime(n + convert(convert(n, base, 10), `+`)), [$1..10^4]); # Robert Israel, Aug 10 2014
MATHEMATICA
Select[Prime[Range[500]], PrimeQ[#+Total[IntegerDigits[#]]]&] (* Harvey P. Dale, Oct 03 2011 *)
PROG
(PARI) select( is(p)=isprime(p+sumdigits(p))&&isprime(p), primes([0, 2000])) \\ M. F. Hasler, Aug 08 2014, edited Nov 09 2018
(Haskell)
a048519 n = a048519_list !! (n-1)
a048519_list = map a000040 $ filter ((== 1) . a010051' . a065073) [1..]
-- Reinhard Zumkeller, Sep 27 2014
(Magma) [p: p in PrimesUpTo(1200) | IsPrime(q) where q is p+&+Intseq(p)]; // Vincenzo Librandi, Jan 30 2018
CROSSREFS
Cf. A007953 (digit sum), A062028 (n + digit sum of n), A047791 (A062028(n) is prime), A048520.
Sequence in context: A145482 A167497 A068579 * A288878 A102907 A235479
KEYWORD
nonn,base
AUTHOR
Patrick De Geest, May 15 1999
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 18:04 EDT 2024. Contains 371254 sequences. (Running on oeis4.)