login
A162250
Values of the form prime(prime(i)) with a prime digital sum.
0
3, 5, 11, 41, 67, 83, 157, 179, 191, 241, 283, 331, 353, 401, 461, 599, 739, 773, 797, 919, 991, 1031, 1217, 1297, 1433, 1471, 1499, 1523, 1723, 1741, 1787, 2027, 2063, 2081, 2221, 2269, 2351, 2609, 2647, 2683, 2719, 2803, 3019, 3109, 3169, 3259, 3299
OFFSET
1,1
FORMULA
{A006450(i) : A007953(A006450(i)) in A000040}. [From R. J. Mathar, Aug 03 2009]
EXAMPLE
Prime(prime(6)) = 41. 4+1=5, prime. So 41 is in the sequence.
PROG
(PARI) sodip(n) = {
local(s=0, a, x, y, j, p);
for(x=1, n, p=prime(prime(x)); a=eval(Vec(Str(p))); y=sum(j=1, length(a), a[j]); if(isprime(y), print1(p", ")); )
}
CROSSREFS
Intersection of A006450 and A028834.
Sequence in context: A035345 A254401 A174915 * A055511 A236457 A105236
KEYWORD
nonn,base
AUTHOR
Cino Hilliard, Jun 28 2009
EXTENSIONS
Definition rephrased by R. J. Mathar, Sep 11 2009
STATUS
approved