cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A066082 Prime octets: numbers k such that 210*k - 105 +- 2^j are prime for all 1 <= j <= 4.

This page as a plain text file.
%I A066082 #17 Jun 17 2018 09:46:41
%S A066082 242590,1175444,2416288,2583146,2596049,2796151,4953911,5574794,
%T A066082 6127655,6396209,6460877,6625438,8521234,11025856,11352491,15482298,
%U A066082 16228703,18861024,19048003,20043534,22835193,31519781,34399756
%N A066082 Prime octets: numbers k such that 210*k - 105 +- 2^j are prime for all 1 <= j <= 4.
%H A066082 Harry J. Smith, <a href="/A066082/b066082.txt">Table of n, a(n) for n = 1..100</a>
%e A066082 a(1)=242590 because 210*a(1) - 105 = 50943795 and 50943795 -2, - 4, - 8, - 16, + 2, + 4, + 8 and + 16 are all prime.
%t A066082 With[{c=2^Range[4]},Select[Range[35*10^6],AllTrue[Flatten[210#-105+ {c,-c}], PrimeQ]&]] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Jun 10 2015 *)
%o A066082 (PARI) { n=0; for (m=1, 10^12, b=210*m - 105; c=0; for (j=1, 4, if (isprime(b - 2^j) , c++, break)); if (c<4, next); for (j=1, 4, if (isprime(b + 2^j) , c++, break)); if (c == 8, write("b066082.txt", n++, " ", m); if (n==100, return)) ) } \\ _Harry J. Smith_, Nov 11 2009
%Y A066082 A proper subset of A061671 (prime sextets). Cf. prime quadruples A014561. See also A066081.
%K A066082 nonn
%O A066082 1,1
%A A066082 _Frank Ellermann_, Dec 03 2001
%E A066082 More terms from _Don Reble_, Dec 07 2001