OFFSET
1,1
COMMENTS
From first 10^5 primes, only 92 are terms. Indices of primes are 1252, 1451, 2032,..., 95460, 97950, 98973.
Note that p == q == r (mod 6), e.g., {10193, 10211, 10223} == 5 mod 6 and {17683, 17707, 17713} == 1 mod 6.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
a(1) = p = 10193; s = {p, q, r} = {10193, 10211, 10223};
{{1,2,2}.s, {2,1,2}.s, {2,2,1}.s} = {51061, 51043, 51031} all primes.
MATHEMATICA
pr=Partition[Prime[Range[40000]], 3, 1]; Reap[Do[s=pr[[k]]; If[Union[PrimeQ[{{1, 2, 2}.s, {2, 1, 2}.s, {2, 2, 1}.s}]]=={True}, Sow[Prime[k]]], {k, Length[pr]}]][[2, 1]]
tcpQ[{p_, q_, r_}]:=AllTrue[{p+2q+2r, 2p+q+2r, 2p+2q+r}, PrimeQ]; Select[ Partition[ Prime[Range[36000]], 3, 1], tcpQ][[All, 1]] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Feb 03 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Mar 15 2013
STATUS
approved