OFFSET
1,1
COMMENTS
The first 3 terms were found by Rotkiewicz.
The generated tetrahedral pseudoprimes are 3776730328549, 4765143438329, 15680770945781, ...
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Andrzej Rotkiewicz, On some problems of W. Sierpinski, Acta Arithmetica, Vol. 21 (1972), pp. 251-259.
EXAMPLE
1179 is in the sequence since 8*1179+1 = 9433, 12*1179+1 = 14149 = 107^2 + 27*10^2 and 24*1179+1 = 28297 = 163^2 + 27*8^2 are primes.
MATHEMATICA
sqQ[n_] := n>0 && IntegerQ[Sqrt[n]]; sqsumQ[n_] := PrimeQ[n] && False =!= Reduce[ x^2 + 27 y^2 == n, {x, y}, Integers]; aQ[n_] := PrimeQ[8n+1] && sqsumQ[12n+1] && sqsumQ[24n+1]; Select[Range[100000], aQ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 20 2018
STATUS
approved