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.
%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