|
%I #13 Jun 21 2018 01:59:37
%S 1,130,253,1222,1657,2437,2476,3589,5188,5395,6073,9763,10063,10225,
%T 10687,11176,11677,13297,13399,13522,13588,14008,14470,14614,16168,
%U 16276,16531,16591,20164,20620,21418,21445,21775,23734,25351,27049
%N Numbers k such that 100k+1, 100k+3, 100k+7, 100k+9 are all primes.
%H Harry J. Smith, <a href="/A064687/b064687.txt">Table of n, a(n) for n = 1..1000</a>
%t Select[Range[30000], PrimeQ[100# + 1] && PrimeQ[100# + 3] && PrimeQ[100# + 7] && PrimeQ[100# + 9] &]
%t Select[Range[30000],AllTrue[100#+{1,3,7,9},PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Sep 28 2015 *)
%o (PARI) for(n=1,100000, if(isprime(100*n+1) && isprime(100*n+3) && isprime(100*n+7) && isprime(100*n+9),print1(n,",")))
%o (PARI) { n=0; for (m=1, 10^9, if(isprime(100*m + 1) && isprime(100*m + 3) && isprime(100*m + 7) && isprime(100*m + 9), write("b064687.txt", n++, " ", m); if (n==1000, break)) ) } \\ _Harry J. Smith_, Sep 22 2009
%K nonn
%O 1,2
%A _Jason Earls_, Oct 10 2001
|
