%I #17 Sep 08 2022 08:45:19
%S 0,2,8,12,24,42,44,48,50,54,68,78,80,84,92,108,110,120,134,140,162,
%T 164,168,174,182,188,194,210,218,222,230,240,264,272,278,282,300,308,
%U 314,320,330,338,342,344,350,362,378,392,402,408,414,420,422,440,444,458
%N Numbers k such that 29*k + 31 is prime.
%H G. C. Greubel, <a href="/A108978/b108978.txt">Table of n, a(n) for n = 1..10000</a>
%e If k=0 then 29*k + 31 = 31 (prime).
%e If k=92 then 29*k + 31 = 2699 (prime).
%p a:=proc(n) if isprime(29*n+31)=true then n else fi end: seq(a(n),n=0..530); # _Emeric Deutsch_, Aug 04 2005
%t Select[Range[0,500],PrimeQ[29#+31]&] (* _Harvey P. Dale_, Jan 13 2011 *)
%o (PARI) is(n)=isprime(29*n+31) \\ _Charles R Greathouse IV_, Jun 13 2017
%o (Magma) [n: n in [0..500] | IsPrime(29*n+31)]; // _G. C. Greubel_, May 19 2019
%o (Sage) [n for n in (0..500) if is_prime(29*n+31)] # _G. C. Greubel_, May 19 2019
%o (GAP) Filtered([0..500], k-> IsPrime(29*k+31) ) # _G. C. Greubel_, May 19 2019
%K nonn,easy
%O 1,2
%A _Parthasarathy Nambi_, Jul 27 2005
%E More terms from _Emeric Deutsch_, Aug 04 2005
|