|
|
A167414
|
|
Primes p such that sum of their squared decimal digits, minus one, is prime.
|
|
2
|
|
|
2, 127, 149, 211, 251, 271, 277, 347, 419, 457, 491, 521, 523, 541, 547, 587, 727, 743, 853, 857, 941, 1021, 1049, 1061, 1087, 1201, 1223, 1229, 1249, 1289, 1373, 1409, 1423, 1429, 1483, 1553, 1559, 1601, 1621, 1627, 1733, 1753, 1861, 1867, 1951, 1973
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
127 is a term because 1^2 + 2^2 + 7^2 - 1 = 53 (prime);
149 is a term because 1^2 + 4^2 + 9^2 - 1 = 97;
347 is a term because 3^2 + 4^2 + 7^2 - 1 = 73.
|
|
LINKS
|
|
|
MAPLE
|
A003132 := proc(n) local d; add(d^2, d=convert(n, base, 10)) ; end proc: A167414 := proc(n) local p; if n = 1 then 2; else p := nextprime(procname(n-1)) ; while not isprime(A003132(p)-1) do p := nextprime(p) ; end do ; return p end if ; end proc: seq(A167414(n), n=1..80) ; # R. J. Mathar, Nov 04 2009
|
|
MATHEMATICA
|
Select[Prime[Range[3000]], PrimeQ[Total[IntegerDigits[#]^2] - 1] &] (* Vincenzo Librandi, Sep 25 2014 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base
|
|
AUTHOR
|
|
|
EXTENSIONS
|
2 and terms larger than 941 added by R. J. Mathar, Nov 04 2009
|
|
STATUS
|
approved
|
|
|
|