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.

A137626 The largest prime in the first set of n consecutive primes for which p+4 is semiprime.

This page as a plain text file.
%I A137626 #38 Jun 22 2018 10:40:43
%S A137626 2,31,181,733,1777,8363,8369,19273,175333,175349,33952819,4377722977,
%T A137626 4377723013,1242030992717,1242030992723
%N A137626 The largest prime in the first set of n consecutive primes for which p+4 is semiprime.
%C A137626 a(n) = last prime in the first run of n primes such that p+4 is semiprime for each prime p in the run. - _Sean A. Irvine_, Feb 13 2012
%C A137626 a(n) > 5 * 10^9 for n > 13.
%e A137626 a(2)=31 is the largest in a set of 2 consecutive primes {29,31}, and 29 + 4 = 33 = 3*11 and 31 + 4 = 35 = 5*7 are both semiprime.  No smaller number has this property.
%e A137626 59 is not in the sequence because although 47 + 4 = 51 = 3*17 and 53 + 4 = 57 = 3*19 are both semiprime, 59 + 4 = 63 = 3*3*7 is not.
%t A137626 With[{prs=Table[If[PrimeOmega[n+4]==2,1,0],{n,Prime[Range[21*10^5]]}]}, Prime[ #]&/@Flatten[Table[SequencePosition[prs,PadRight[{},n,1],1],{n,11}],1]][[All,2]] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, May 10 2018 *)
%o A137626 (PARI) a(n) = {my(t = 0); forprime(p = 2, oo, if(bigomega(p + 4) == 2, t++; if(t==n, return(p)), t = 0))} \\ _David A. Corneth_, May 10 2018
%Y A137626 Cf. A001358, A137625, A137627, A137628.
%Y A137626 Subsequence of A289250.
%K A137626 more,nonn
%O A137626 1,1
%A A137626 _Enoch Haga_, Jan 30 2008
%E A137626 a(11) from _Sean A. Irvine_, Feb 12 2012
%E A137626 a(1) corrected by _Harvey P. Dale_, May 10 2018
%E A137626 a(12)-a(13) from _David A. Corneth_, May 10 2018
%E A137626 a(14)-a(15) from _Giovanni Resta_, Jun 22 2018