Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #9 Sep 08 2022 08:46:11
%S 2,3,7,19,37,43,79,113,151,163,193,229,251,281,317,373,397,433,463,
%T 503,577,757,827,911,953,997,1069,1123,1321,1399,1423,1481,1657,1693,
%U 1747,2029,2143,2267,2311,2473,2503,2551,2593,2687,2753,2791,2917,3313,3323
%N Prime(n) is included iff prime(n) + n^2 is also prime.
%F a(n) = prime(A064711(n)). - _Michel Marcus_, Feb 04 2015
%e 7 is in this sequence because 7+16=23.
%e 19 is in this sequence because 19+64=83.
%t Prime[Select[Range[500], PrimeQ[Prime[#] + #^2] &]]
%o (Magma) [NthPrime(n): n in [1..500] | IsPrime(NthPrime(n)+n^2)];
%o (PARI) lista(nn) = forprime (n=2, nn, if (isprime(n+primepi(n)^2), print1(n, ", "))); \\ _Michel Marcus_, Feb 04 2015
%Y Cf. A064711, A061067.
%K nonn
%O 1,1
%A _Vincenzo Librandi_, Feb 04 2015