|
| |
|
|
A046704
|
|
Additive primes: sum of digits is a prime.
|
|
37
| |
|
|
2, 3, 5, 7, 11, 23, 29, 41, 43, 47, 61, 67, 83, 89, 101, 113, 131, 137, 139, 151, 157, 173, 179, 191, 193, 197, 199, 223, 227, 229, 241, 263, 269, 281, 283, 311, 313, 317, 331, 337, 353, 359, 373, 379, 397, 401, 409, 421, 443, 449, 461, 463, 467, 487, 557, 571, 577, 593
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,1
|
|
|
LINKS
| Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
|
|
|
MATHEMATICA
| Select[Prime[Range[100000]], PrimeQ[Apply[Plus, IntegerDigits[ # ]]]&]
|
|
|
PROG
| (PARI) isA046704(n)={local(s, m); s=0; m=n; while(m>0, s=s+m%10; m=floor(m/10)); isprime(n) & isprime(s)} [From Michael Porter (michael_b_porter(AT)yahoo.com), Oct 18 2009]
(MAGMA) [ p: p in PrimesUpTo(600) | IsPrime(&+Intseq(p)) ]; // Bruno Berselli, Jul 08 2011
(Haskell)
a046704 n = a046704_list !! (n-1)
a046704_list = filter ((== 1) . a010051 . a007953) a000040_list
-- Reinhard Zumkeller, Nov 13 2011
|
|
|
CROSSREFS
| Indices of additive primes are in A075177.
Cf. A046703.
Cf. A007953, A010051; intersection of A028834 and A000040.
Sequence in context: A087521 A078403 A129945 * A089392 A089695 A070027
Adjacent sequences: A046701 A046702 A046703 * A046705 A046706 A046707
|
|
|
KEYWORD
| base,nonn
|
|
|
AUTHOR
| Felice Russo (frusso(AT)micron.com)
|
| |
|
|
