A296507 - OEIS

%I #19 Mar 06 2023 09:55:23

%S 4,6,12,18,24,30,36,54,72,84,90,96,102,114,120,138,168,186,198,204,

%T 210,216,228,240,276,294,318,330,354,360,372,378,402,414,438,444,456,

%U 480,498,504,588,600,612,618,630,636,666,678,690,714,720,726,732,738,762

%N Numbers m such that m^2 - 13 is a prime.

%C All terms except 4 are divisible by 6. - _Robert Israel_, Dec 13 2017

%H Daniel Starodubtsev, <a href="/A296507/b296507.txt">Table of n, a(n) for n = 1..5000</a>

%p select(n -> isprime(n^2-13), 2*[$2..10^4]); # _Robert Israel_, Dec 13 2017

%t Reap[m=4;Do[If[PrimeQ[m^2-13],Sow[m]];m=m+2,{1000}]][[2,1]]

%t Select[Range[800],PrimeQ[#^2-13]&] (* _Harvey P. Dale_, Mar 06 2023 *)

%o (PARI) isok(n) = isprime(n^2-13); \\ _Michel Marcus_, Dec 14 2017

%Y Cf. A028870, A028873, A028876, A028882, A090696.

%K nonn

%O 1,1

%A _Zak Seidov_, Dec 13 2017