OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
The first three primes == 1 (mod 4) are 5, 13, 17, but 5+13+17=35 is not prime.
The second, third and fourth primes == 1 (mod 4) are 13, 17, 29, and 13+17+29=59 is prime, so a(1) = 59.
MAPLE
N:= 2000: # to use primes <= N that == 1 (mod 4)
P:= select(isprime, [seq(i, i=1..N, 4)]):
select(isprime, P[1..-3]+P[2..-2]+P[3..-1]);
MATHEMATICA
M = 2000;
P = Select[Range[1, M, 4], PrimeQ];
Select[P[[1;; -3]] + P[[2;; -2]] + P[[3;; -1]], PrimeQ] (* Jean-François Alcover, Apr 24 2019, from Maple *)
Select[Total/@Partition[Select[Prime[Range[250]], Mod[#, 4]==1&], 3, 1], PrimeQ] (* Harvey P. Dale, Aug 20 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, May 11 2018
STATUS
approved