%I #21 Sep 08 2022 08:46:08

%S 4,10,100,1480,16060,19420,21010,22270,43780,55330,144160,165700,

%T 166840,195730,201820,225340,247600,268810,326140,347980,361210,

%U 397750,465160,518800,536440,633460,633790,661090,768190,795790,829720,857950,876010,958540

%N Numbers m such that m+1, m+3, m+7, m+9 and m+13 are all primes.

%D W. Sierpiński, 250 Problems in Elementary Number Theory. New York: American Elsevier, 1970. Problem #82, variant.

%H Reinhard Zumkeller and Jens Kruse Andersen, <a href="/A245304/b245304.txt">Table of n, a(n) for n = 1..1000</a> (first 120 terms from Zumkeller)

%F a(n) = A022006(n)-1. - _Jens Kruse Andersen_, Jul 18 2014

%t Select[Range[10^6],AllTrue[#+{1,3,7,9,13},PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, May 07 2015 *)

%o (Haskell)

%o a245304 n = a245304_list !! (n-1)

%o a245304_list = map (pred . head) $ filter (all (== 1) . map a010051') $

%o iterate (zipWith (+) [1, 1, 1, 1, 1]) [1, 3, 7, 9, 13]

%o (PARI) forprime(p=2, 10^7, m=p-1; if(isprime(m+3)&&isprime(m+7)&&isprime(m+9)&&isprime(m+13), print1(m", "))) \\ _Jens Kruse Andersen_, Jul 18 2014

%o (Magma) [n: n in [0..10^6] | IsPrime(n+1) and IsPrime(n+3) and IsPrime(n+7) and IsPrime(n+9) and IsPrime(n+13)]; // _Vincenzo Librandi_, Jun 15 2015

%Y Cf. A010051, A022006, A245305, A007811, subsequence of A125855.

%K nonn

%O 1,1

%A _Reinhard Zumkeller_, Jul 18 2014