OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
100 is in the sequence because the 100th prime is 541, and 100^2 + 541^2 = 302681 is prime.
MAPLE
with(numtheory): nn:= 150: T:=array(1..nn):k:=1:for n from 1 to 1764 do:p:=ithprime(n):if type(p^2+n^2, prime)=true then T[k]:=n:k:=k+1: else fi:od:print(T):
MATHEMATICA
Select[Range[600], PrimeQ[#^2 + Prime[#]^2] &] (* Amiram Eldar, Aug 11 2019 *)
PROG
(Magma) [m:m in [1..600] | IsPrime(m^2+NthPrime(m)^2)]; // Marius A. Burtea, Aug 11 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, May 03 2010
STATUS
approved