A046121 Duplicate of A023271.
5, 11, 41, 61, 251, 601, 641, 1091, 1481, 1601, 1741, 1861, 2371, 2671, 3301, 3911
Offset: 1
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.
251, 257, 263, 269 are consecutive primes: 257 = 251 + 6, 263 = 251 + 12, 269 = 251 + 18.
N:=10^5: # to get all terms <= N. Primes:=select(isprime,[seq(i,i=3..N+18,2)]): Primes[select(t->[Primes[t+1]-Primes[t], Primes[t+2]-Primes[t+1], Primes[t+3]-Primes[t+2]]=[6,6,6], [$1..nops(Primes)-3])]; # Muniru A Asiru, Aug 04 2017
A033451 = Reap[ For[p = 2, p < 100000, p = NextPrime[p], p2 = NextPrime[p]; If[p2 - p == 6, p3 = NextPrime[p2]; If[p3 - p2 == 6, p4 = NextPrime[p3]; If[p4 - p3 == 6, Sow[p]]]]]][[2, 1]] (* Jean-François Alcover, Jun 28 2012 *) Transpose[Select[Partition[Prime[Range[16000]],4,1],Union[ Differences[ #]] == {6}&]][[1]] (* Harvey P. Dale, Jun 17 2014 *)
p=2;q=3;r=5;forprime(s=7,1e4,if(s-p==18 && s-q==12 && s-r==6, print1(p", ")); p=q;q=r;r=s) \\ Charles R Greathouse IV, Feb 14 2013
p = 409 then the AP-9 is {409, 619, 829, 1039, 1249, 1459, 1669, 1879, 2089} with the difference 9# = 2*3*5*7 = 210.
Clear[p]; d = 210; ap9p = {}; Do[If[PrimeQ[{p, p + d, p + 2*d, p + 3*d, p + 4*d, p + 5*d, p + 6*d, p + 7*d, p + 8*d}] == {True, True, True, True, True, True, True, True, True}, AppendTo[ap9p, p]], {p, 3, 10^9, 2}]; ap9p
v=[1..8]*210; forprime(p=1,,for(i=1,#v,isprime(p+v[i])||next(2));print1(p",")) \\ M. F. Hasler, Jan 02 2020
p = 11 then {11, 11 + 1*30, 11 + 2*30, 11 + 3*30, 11 + 4*30} = {11, 41, 71, 101, 131}, which is 5 primes in arithmetic progression with the difference 5# = 30.
Clear[p]; d = 30; ap5p = {}; Do[If[PrimeQ[{p, p + d, p + 2*d, p + 3*d, p + 4*d}] == {True, True, True, True, True}, AppendTo[ap5p, p]], {p, 3, 25000, 2}]; ap5p
lst={};Do[p=Prime[n];If[PrimeQ[p+6]&&PrimeQ[p+12]&&PrimeQ[p+18], AppendTo[lst, p+6]], {n, 8!}];lst (* Vladimir Joseph Stephan Orlovsky, Aug 29 2008 *)
lst={};Do[p=Prime[n];If[PrimeQ[p+6]&&PrimeQ[p+12]&&PrimeQ[p+18], AppendTo[lst, p+12]], {n, 8!}];lst (* Vladimir Joseph Stephan Orlovsky, Aug 29 2008 *) Select[Prime[Range[3000]],AllTrue[#+{-12,-6,6},PrimeQ]&] (* Harvey P. Dale, Jan 25 2023 *)
[p+18: p in PrimesUpTo(30000) | IsPrime(p+6) and IsPrime(p+12) and IsPrime(p+18)]; // Vincenzo Librandi, Jan 07 2015
lst={};Do[p=Prime[n];If[PrimeQ[p+6]&&PrimeQ[p+12]&&PrimeQ[p+18], AppendTo[lst, p+18]], {n, 8!}];lst (* Vladimir Joseph Stephan Orlovsky, Aug 29 2008 *)
p = 179 then the AP-5 is {179, 389, 599, 809, 1019, 1229, 1439} with the difference 7# = 210.
Clear[p]; d = 210; ap7p = {}; Do[If[PrimeQ[{p, p + d, p + 2*d, p + 3*d, p + 4*d, p + 5*d, p + 6*d}] == {True, True, True, True, True, True, True}, AppendTo[ap7p, p]], {p, 3, 10^9, 2}]; ap7p Select[Prime[Range[15000]],And@@PrimeQ[NestList[210+#&,#,6]]&] (* Harvey P. Dale, Nov 16 2013 *)
is(p)=forstep(k=p,p+1260,210,if(!isprime(k),return(0)));1 \\ Charles R Greathouse IV, Dec 19 2013
Clear[p]; d = 210; ap8p = {}; Do[If[PrimeQ[{p, p + d, p + 2*d, p + 3*d, p + 4*d, p + 5*d, p + 6*d, p + 7*d}] == {True, True, True, True, True, True, True, True}, AppendTo[ap8p, p]], {p, 3, 3000000, 2}]; ap8p Select[Prime[Range[260000]],AllTrue[NestList[#+210&,#,7],PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, May 03 2018 *)
p = 186874511 then the AP-11 is {186874511, 186876821, 186879131, 186881441, 186883751, 186886061, 186888371, 186890681, 186892991, 186895301, 186897611} with the difference 11# = 2*3*5*7*11 = 2310.
Clear[p]; d = 2310; ap11p = {}; Do[If[PrimeQ[{p, p + d, p + 2*d, p + 3*d, p + 4*d, p + 5*d, p + 6*d, p + 7*d, p + 8*d, p + 9*d, p + 10*d}] == {True, True, True, True, True, True, True, True, True, True, True}, AppendTo[ap11p, p]], {p, 3, 40*10^9, 2}]; ap11p ap11Q[n_]:=AllTrue[Rest[NestList[2310+#&,n,10]],PrimeQ]; Select[Prime[ Range[ 148*10^7]],ap11Q] (* The program uses the AllTrue function from Mathematica version 10 *) (* The program will take a long time to run *) (* Harvey P. Dale, Oct 27 2019 *)
Comments