OFFSET
1,1
COMMENTS
This is to Pythagorean primes (A002144), primes of the form 4n+1, as primes with prime subscripts (A006450) is to primes (A000040). Hence this is one of four related sequences into which every prime with prime subscripts (A006450) may be classified: Pythagorean primes (A002144) with Pythagorean prime index; Pythagorean primes (A002144) whose indices are of the form 4n+3 (A002145); primes of the form 4n+3 with Pythagorean prime index; and primes of the form 4n+3 whose indices are primes of form 4n+3.
EXAMPLE
a(1) = 37 because the smallest prime of form 4n+1 is 4*1+1 = 5, and the fifth prime of the form 4n+1 is 4*9+1 = 37. a(2) = 109 because the second smallest prime of form 4n+1 is 4*3+1 = 13, and the thirteenth prime of the form 4n+1 is 4*27+1 = 109.
MATHEMATICA
A002144=Select[Range[5, 100000, 4], PrimeQ];
PROG
(PARI) c=0; forprime(p=2, 10^5, if(p%4==1, c++; if(isprime(c)&&c%4==1, print1(p, ", ")))) \\ Max Alekseyev, Feb 08 2010
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jonathan Vos Post, Dec 19 2009
EXTENSIONS
More terms from Max Alekseyev, Feb 08 2010
STATUS
approved