|
%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
|
